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STEADY STATE SOLUTIONS OF TRANSMISSION 
LINE EQUATIONS 


S. O. Rice 


Methods of obtaining the steady state voltages and currents in a 
uniform transmission line consisting of several parallel wires are described 
in Part I. This line may or may not be acted upon by an externally 
impressed field distributed along its length. A square matrix I’, whichisa 
generalization of the propagation constant y for a single circuit, is intro- 
duced. Matrix expressions obtained for the voltages and currents in- 
volve I in much the same way as the corresponding single circuit expres- 
sions involve y. In Part II similar methods are described for obtaining 
the voltages and currents in a transmission line composed of a number of 
multi-terminal symmetrical sections connected in tandem. Expressions 
for the voltages and currents in a line composed of unsymmetrical sections 
are also given. ‘These sections may or may not contain generators. 


HE transmission lines considered here are of two kinds, namely the 

uniform transmission line, and the transmission line consisting of a 
number of identical sections connected in tandem. The problem discussed 
is that of determining the steady state electrical behavior of these lines 
when the terminal conditions are given. Often there arises the problem 
of determining the currents induced in a uniform transmission line by an 
arbitrary impressed field of some fixed frequency or of determining the cur- 
rents produced by generators placed in the branches of the sections if the 
line is of the second kind. This is the type of problem with which we shall 
be particularly concerned. 

In dealing with the uniform transmission line it is found convenient to 
introduce a matrix I’, which is a generalization of the propagation constant ¥ 
for a single wire with ground return, or for a single circuit. This enables us 
to obtain matrix expressions for the currents and voltages which are similar 
in form to the single circuit expressions. 

A similar situation exists for the transmission line composed of a number 
of symmetrical sections. However, when the sections are unsymmetrical 
the corresponding procedure does not appear to yield a corresponding sim- 
plification and the formulas are considerably more complicated than in the 
symmetrical case. 

This paper is divided into two parts corresponding to the two kinds of 
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transmission lines. The first part discusses the uniform line. After q 
statement of the transmission equations in matrix form, expressions for the 
voltages and currents are given. Two methods of evaluating these ex. 
pressions are described. The first is based upon a property possessed by 
many transmission systems, namely that the various modes of propagation 
have nearly the same speed. The second method is based upon equations 
which may be obtained by the formal application of a theorem due to 
Sylvester. The first part concludes with the proof that these two methods 
lead to the correct results. 

After a short introduction the second part discusses the difference equa- 
tions which govern the transmission in a line composed of multi-termina! 
sections. The sections may contain generators. Expressions for the volt- 
ages and currents in a symmetrical section line, i.e. a line whose sections are 
symmetrical, are stated and proved in much the same order as the corre- 
sponding expressions for the uniform line. A discussion of the unsymmetri- 
cal section line concludes the second part. 

A sketch of the solution of the uniform transmission line equations by 
the classical method is given in Appendix I. In Appendices IT and III 
methods are described for solving the symmetrical section line difference 
equations. These methods are similar to the one of Appendix I. The 
method of Appendix III uses section constants which may be obtained from 
measurements made at one end of a typical section. 


Part I 


UNIFORM TRANSMISSION LINES 


1.1 Differential Equations 

For the sake of convenience in writing down equations we shall assume 
that the particular line under consideration consists of three parallel wires 
with ground return, or of three parallel circuits, denoted by the subscripts 
a, 6, and c respectively. The differential equations for this line in an arbi- 
trary impressed field are’ 





dq = —Zaata aes Za to ig + ee + 1,(x) 

dx 

dv, - . 7 ° 4 

— = —Z,i, - Lah — Ents + dle) (1.1) 
dx 

dv. = —Zeata —_ Zats ine Leote + 1.(x) 

dx 


1 These equations are given in substance by J. R. Carson and R. S. Hoyt, B.S.T.J/., 
Vol. 6, pp. 495-545 (1927). Equations (1.2) are equivalent to their equation (90) and 
equations (1.1) ray be obtained by combining their equations (83), (84), and (94). We 
shall use the term “impressed field” to mean a field distributed along the line. According 
to our convention there is no impressed field when the line is energized only at the terminals. 
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and 
di ‘ ; . 
-=-— y aaVq — } ab Up — Vacte + ta(x) 
dx 
diy : . . 
— J baVa — } bb Up ~~ } te Ve + ty(x) (1.2) 
dx 
di : ' A 
: == } cala — } Vy — } ee ve + t.(x) 
dx 
where Zan = Zora, Van = Vou, etc. If we are dealing with three parallel 


wires Ja(x), Us(x), 1-(x) are the longitudinal components of the electric force 
of the impressed field at the wire surfaces; ¢,(x), (x), ¢-(x) are specified 
by the admittance of the direct leakage paths and the values of the im- 
pressed potentials at the wires. If there are no direct leakage paths the 
t's are zero. 

In order to put these equations in matrix form? we introduce the column 


matrices 


Va ta 1,(x) t,(x) 
v={ 0%], i= iy , Ua) = 1,(x) P t(x) = ty(x) 2 (1.3) 
Ve i. 1.(x) t.(x) 


and the symmetrical square matrices 


Zea Zab 2. Fes Yoo Ves 
Z= |Z Zo Zr Y=] Vo. Vo Yee (1.4) 
Lea Zeb Zee Vcc Ye» Vec 


The equations (1.1) and (1.2) may now be written as 


dv 
dx 
di 


= —Zi + I(x) 


and these are the equations to be solved. 
When there is no impressed field equations (1.5) give 
d’v 
dx? 
d’i 
dx? 


2 Cf. L. A. Pipes, Phil. Mag., Vol. 24 (1937), p. 97. 


= ZYv 
(1.6) 
= YZ 
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and the analogy with the one circuit case leads us to put 
m=ZY, r=~</ZY 17 


where I is a square matrix representing a generalization of the propagation 
constant. Putting aside for the moment the question of interpreting the 
square root, we note that interchanging the rows and columns in I = Zy 
gives 


[2 = y’zZ' ae ga T’ wa VVZ 1.8 


where the primes denote transposition. Y’ and Z’ are equal to Y and Z 
respectively because of their symmetry. We thus expect I” to be associated 
with the propagation of i in the same way that T is associated with the 
propagation of 2. 
1.2 Statement of Results for an Infinite Line—No Impressed Field 

It is shown that when there is no impressed field the voltages and currents 
at any point x in a transmission line extending from x = 0 tox = &« are 
given by 


—2rT Ty: 
v(x) = e€ * v0) =e” Z,i(0) 
‘ —sI'’. 
i(x) = €” 1t(0) (1.9 
v(x) = Z,i(x) 


where ¢ * is the square matrix defined by the convergent series of matrices 


allied I xT 4 x I“ x" 4 4 im 
ee eee ae. t 
T . . r . 
I denotes the unit matrix. Z, is a 
square matrix and is called the characteristic impedance matrix: 


r’ 
and e* is the transposed of e * 


Z=V'ZeaQryY° (1.11 


Additional expressions of the same type for Z, are given by equations 
(1.45). The matrix e *'Z,, being of the nature of a transfer impedance, 
is symmetrical. 

The matrices e *' and Z, may be computed in several ways, the choice 
depending upon the circumstances. The first method to be described is 
useful when «x is not too large and when the propagation constants of the 
various modes of propagation are nearly equal to each other. In the case 
of open-wire lines these propagation constants are grouped around the 
value jw/v where v is of the order of 180,000 miles per second. The second 
method may be used for all cases, including those for which the series in 


3 Frazer, Duncan and Collar, “Elementary Matrices,” Cambridge University Press, 


§2.5. In the work which follows, this text will be referred to as “F.D.C.” 
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the first method converge too slowly to be of value. However, it requires 
the solution of an mth degree equation and the determination of the m 
modes of propagation where m is the number of circuits. 
this is no handicap and the method is quite convenient. In 


For m ) 


this case 
the method is closely related to one described by John Riordan in an un- 
published memorandum. 

First Method: Multiply the matrices Z and JV together to obtain ZY. 
Choose the number 7¥? in 


ZY = Iy* + R, (1.12) 


where J is the unit matrix, so that the elements of R are small in comparison 
For many transmission lines it is possible to do this. I may be 
obtained by using the binomial theorem to expand the square root in the 


with 7’. 


formula 
P= VZY = y(I + y°R)', 
where y is that square root of y? whose real and imaginary parts are non- 


negative. In carrying out the work it is convenient to introduce the matrix 
S whose elements are small in comparison with unity. 


(1.13) 


r= y+ S) (1.14) 


To compute S, first compute the matrix R/2y? and then use the power 


s= (4) -3(&) +3(%) -2(8) 
: ay?) 2 \2y? 27? 8 \2y? 
772 wvnizY 
+i(¥) : 3 (3) Mg 


This series will usually converge rapidly. 


series 


Nl 


(1.15) 


rT 


The matrix ¢ is given by 


or asa (1.16) 
where z is a number, z = yx, and e ™ is to be computed from 
f v\ 2 4 +\3 
28 ZS (2S) ( 2S) 
= _ — _ eee ( } 
€ I ii + 7 3) + (1.17 


rs ‘ r M - as 
e* is obtained from e * by interchanging the rows and columns. The 
characteristic impedance matrix may be obtained from (1.11), 


after computing I’ from S as in (1.14). 
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ef 


If only e “" is required the following series may be used. 


mer. Qa fReY 5,(s) 
euliens 2 (#) p! _- 


where R, y, and z have the same meaning as above and the coefficients ar 
computed from 


=e", be) = (1 + 4) 


2 1 
bpi2(z) = b,(z) — a bp41(z) 


: : Rx\° 
In the first term of the series (3 ) denotes J. 
<Y 


. ‘ r 2 . ‘ 

Second Method: TI’, e *" and Z, may be regarded as functions of the squar 
matrix ZY. In order to express these functions in a form suitable for calcu- 
. . ; 4 - . ° * cs : 
lation we apply Sylvester’s theorem’. The characteristic matrix of Z)’ is 

fq) = 771 — ZY (1.19) 
where now y° is regarded as a variable instead of a fixed number as in the 
first method. We shall suppose that ZY is a square matrix of order m 


2 2 2 or ° ° 
and that the roots yi, y2, --- Ym of the characteristic function, i.e. of the 
determinantal equation 


fo?) | = 0, (1.20 
are distinct. Let the matrix F(y?) be the adjoint of f(y?) and denote the 
derivative of the characteristic function by 

d 
d(y*)- 
Since vi, 7, cee 7, are all different f(y?) Pde unequal! to zero for r 


1, 2, --- m. Sylvester’s theorem says that if P(ZY) is any polynomial 
in ZY then 


| f(y’) |? = f(y’) (1.21 


P(ZY) = > N(y7?)P(7?) (1.22) 
r=] 
where P(y;) is a scalar (and thus deviates from our convention that capital 
letters denote square matrices). N(y?) is a square matrix: 
> ..2 
| F (y;) 
N(y;) = aD 
I (yr) | 
When m = 2, N(y73) is equal to I — N(y}). 


*F.D.C. §3.9. The « and d of the reference are the ZY and 7? of the present section 


wher 
are t 
nega 
mod 
of th 
prop 
asa 
are | 
colu 
corr 
1.3 
W 


exp! 





TRANSMISSION LINE EQUATIONS 137 
Applying (1.22) to I, e =F and Z, even though they are not polynomials 
in ZY gives results which may be verified to be true. 
r= V/ZY = > Ny, 
=¢*V2Z7¥ =} N(y2)e*” 
Z, = (ZY)'Y' = > N(y?)v,¥ 
e at 3 > N(y2)y-e rYr Y 1 


where the summations extend from r = 1 tor = mandy, y2, «++ Ym 
are the square roots of vi P v3 ree 7, respectively whose real parts are non- 
negative. 1, Y2, °°: Ym are also the propagation constants of the ‘‘normal 
modes” of propagation. Some light is thrown on the physical significance 
of the matrix N(y7) by supposing that only the rth normal mode is being 


propagated on the transmission line. V(y,) is such that it can be expressed 


as a column matrix times a row matrix. The voltages in circuits 1, 2, --- m 
are proportional to the first, second, --- mth elements, respectively of the 
column matrix. The currents in circuits 1, 2, --- m are proportional to the 
corresponding elements in the row matrix. 
1.3 Results for Any Uniform Line—No Impressed Field 
When the length of the line is finite the voltages and currents may be 
expressed as 
v(x) = cosh aI v(0) — sinh aI Z,7(0) 
(1.25) 
i(x) = — sinh aI” Z,'v(0) + cosh xI” i(o) 


where Z, and I have the same meaning as before. The matrices sinh aT'Z, 

e chien . ws ‘ ‘ ‘ ‘ 
and sinh xI’Z, are symmetrical. The square matrices cosh xT and sinh aI 
are defined by the series 


cosh xl = 


cT x TY’ 

sinh «I i! + 31 + 
cosh xI” is obtained by interchanging the rows and columns of cosh «I 
and sinh xI” is obtained similarly from sinh xf. Solving (1.25) for v(o) 
and i(0) gives 

v(o) = cosh aI v(x) + sinh aT Z,i(x) 

‘ ; — 3 

i(o) = sinh xI’ Z, v(x) + cosh aI” i(x) 

As in the case of the infinite line, we have two ways of computing the 

coefficients of v(o) and i(o) in the expressions (1.25) for v(x) and i(x). 
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First Method: Choose a number y? and compute the matrices R, S, I’, 7 
as described in the first method for the infinite line. The matrix ¢’' 
given by 


zT s 8S 
c =€-¢ 
where z = yx and e** is computed from the series 
28 25 , 2S 
é = I — _ 
ba 1! ™ 2! ™ 


If the elements of 2S are so large that the series converges slowly it may lx 


; i ' 2S 
worthwhile to divide zS by 16, say, compute exp (;*) from the series, an 
) 


then obtain e*” by four matrix multiplications. When e*” is known. its 


: z8 
inverse e ~ can be computed and e * 


r ‘ ° rar 
obtained from (1.16). The hyper- 
bolic functions are given by 
. r r 
cosh aI = 3 (e&* +e”) 
(1.27 


s ¢2r zT 


sinh al = §$(€ —e”~ ) 


which follow from the series definitions of the various matrices. 
If only the coefficients in (1.25) are required we may choose y? and con 
pute R and powers of the matrix Rx/2y. Then the coefficients in (1.25 


2 p 
cosh «I = D3 (**) a,y(z) 


are given by 


p=0 27 p! 
’ : — ( Rx\" apyilz) , 
sinh al Z, = £(F) oes! Z (1.28 
p=0 \2y ply 
, -  ( R’x\" apiilz) y, 
sinh xI’Z,' = a (% PH ) } 
p=0 \ 2 py 
where R’ is the transposed of R, and the scalar coefficient a,(z) is a function 
of s = yx given by 
a,(z) = cosh z a;(z) = sinh s 
ES ree sinh z 
a2(z) = cosh z — (1.29) 
2p+1 
Gp+2(3) = a,(z) — i — Gp41(z), 


n 


and it is understood that (Rx/2y)° = I. 
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Second Method: Compute the propagation constants y:, Ye, ++: Ym 
and the square matrices V(y,) given by (1.23) as in the second method for 
the infinite line. Then 

. , 2 
cosh xl = J N(y;) cosh xy, 
° “ar ° ymnr!l 
sinh al Z = sinh al T} 
a Se a , 
> N(y;) sinh xy, yr] (1.30 
° sp wl . wo-ly, 
sinh xI’ Z, = sinh aI’ I’ '} 
+7. 2 Sinh xy, ,, 
=D N'(y) } 
Yr 
ew oe oe 2 : . : 
where N’(y,) is the transposed of V(y,), N(y;) being defined by (1.23), 
and the summations extend from r = 1ltor = m. 
When the transmission line consists of perfectly conducting wires strung 
on perfect insulators over a perfectly conducting earth the magnetic and 
- ¢ ‘ 1 
electrostatic fields are related so as to make Z equal toy,)' ~~ where 
Yo a qw/ C, 
w being 24 times the frequency and c¢ the speed of light. 

It is interesting to apply the first method of solution to this line. Even 
though the proof of the first method, which is given in §1.10, does not cover 
this case there seems to be little doubt that the correct answer is obtained. 

We have 

our 2 
ZY = y,1 
Choosing y = Yo gives R = 0 and therefore S = 0. It follows that 
' e wl « aaa 
r=7/,2Z.=T" Z2=%7 2 
cosh xT = cosh (xy.J) = cosh xy, / 
‘ — P = 
sinh aI Z, = sinh xy, 7, Z 
P —— F a | 
sinh xI” Z, = sinh x70 YoZ 


When these are put into equations (1.25) the expressions for v(x) and i(x) 
in a perfect transmission line are obtained: 


sinh x75 ,-, 
sinh 2y Zi(o) 
Yo (1.31) 


v(x) = cosh xy,v(0) — 


i(x) = —-y, sinh xy,Z ‘v(0) + cosh xy, i(o) 




































—4 (°F — 1) (I — OL) (r + ZA) 


, . i ‘ po . r’ 6 7 
provided that the inverse matrices exist. The matrix (e* — e *J)(I’ + 


5 F.D.C. §2.10. 
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1.4 Results for Any Uniform Line—Impressed Field 
The differential equations to be satisfied in this case are given by (1.5 
é A solution which reduces to v(0) and i(o) at x = Ois 
r. v(x) = cosh aI v(0) — sinh aT Z,t(0) 
‘ z z 
&s + cosh (x — &)Pl(é) dé — | sinh (x — §)I Z, ¢(€) dé 
: 0 0 
(1.32 
i(x) = —sinh xI” Z,'v(0) + cosh aI” i(o) 
ax z 
— | sinh (x — &)I'’Z7'I(t) dé + / cosh (x — &)I't(&) dé 
“0 0 
The matrices cosh aI’, sinh aT and Z, are the same as the ones discussed i; 
$1.2 and §1.3. The elements of the integral’ of a matrix U (U is not neces. 
4 sarily a square matrix) are given by the integrals of the corresponding 
ga elements of U. 
In many cases of practical interest the impressed field varies exponential], 
Es with respect to x. The column matrices /(x) and ¢(x) may then be ex 
he pressed as 
e Na Ta 
I 6 2 
Wx) =e "|r i(x) =e"! ty (1.33 
Ne Te 
k: where the \’s and 7’s are constants and @ is the propagation constant of 
's the impressed field in the direction of the line. The integrations in the 
E expressions (1.32) may be performed with the result 
v(x) = cosh aT v(0) — sinh aT Z,i(0) 
. + 3 (e*" — & 7) (CT + 61) "(A — Zor) 
: — 3 (e7" — 6 *71) (TF — 61) "(A + Zor) 
(1.34 
' i(x) = —sinh xI” Z;'v(0) + cosh xI” i(o) 
+4 (er — 67) (I + 61) (r — Z,'d) 


° ° : O@r\/n 1 
61) ' is the transposed of (e*" — e 7)(T + 07)", etc. If one of these 
matrices, say I' — @/, has no inverse then it is necessary to evaluate the 
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corresponding integral in some other way. Thus it may be advantageous 
to use the formula 


_(e rT ws “rr _ 61) 1 =e rT [ ef 1 de 
2 ‘ (1.35) 
. ar a ) 
=¢* E Vale a Se Sy + | 
as Ji 
[wo special cases of (1.34) are of interest. When the line is shorted at 
both ends, v(0) = v(x) = 0, where x is the line length, and 
i(o) = 4 Z,' (sinh aT) ‘[(e"" — e “/)(T + 61) (A — Zor) 


-zT 


— (e —~¢ "I(T — 6!) (x + Zor)] 
i(x) = Z,' (sinh x1) ‘[(e"* — e? 1)(" — 61) (X + Zor) 


— (e*" — 7)(T + Of) (A — Zo7)) 


When the line is terminated in its characteristic impedance at both ends, 
vo) = —Z,i(o), v(x) = Z,1(x), and 


io.) = 3 (1 — &%e*"(r’ + 61) "(Z.'d — 7) 
i(x) = —¥ (07 — PNY’ — OF) (ZX + 17) 


The matrices occurring in the expressions (1.34) for o(«) and i(x) may 
be computed by the first or second method described for the uniform line 
in the absence of an impressed field. The second method involves the use 
of expansions similar to 


(f? — &"1)(T + Of) (A — Zor) 


ry x8 a 

; r— Z 
N r(° : Y(r - ) 
X Ye) Yr + 0 ‘: 


7T’ 6 1 - (1.36) 
—e° It’ + ol) (ZA — 2) 


=> vo (— at a A- r) 
; fe Yr + 0 Yr 


where the summations run from r = 1 tor = mand N'(y?) is the transposed 
of the square matrix N(y;) given by (1.23). In obtaining these expansions 
by Sylvester’s theorem, Z, in the first is replaced by ’Z and Z," in the 
second by IY" Y. 

If we assume that an impressed field acts upon the perfect transmission 
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line of equations (1.31), we see that f(x) = 0 because there are no direc: 
leakage paths. We may also write 


(7 —e @I(r +61)! = ("Il —e*1)(y.1 + 61)” 


er er x6 
= I 
Yo t+ 6 


From this and similar equations it follows that 


sinh xy, 


v(x) = cosh xy,v(0) — Zi(o) 


re 1 = a" —_ TYo ee “]a 
2L Yt 8 Yo — 9 (1.37 


i(x) = —y sinh xy, Z ‘v(0) + cosh XYo1(0) 
ro x6 TYo x6 
sie °° =e e —e 
2 Yo + 6 Ye™ 0 


1.5 Results for Infinite Uniform Line—Impressed Field 
When the line extends from x = 0 tox = » and the impressed field is suc! 


Yo 


that the voltages and currents remain finite at x = 2%, the appropriate solu- 
tions may be obtained from the results of §1.4 by a limiting process. The 
condition that v(x) remain finite suggests that the coefficient of e7" be zero 
in the expression (1.32) for o(x). This gives a relation between v(0) and 
i(0) which must be satisfied: 


v(0) = Z,i(o) — e (l(t) — Z,4(€)| dé (1.38 
0 


If the impressed field varies exponentially with x expression (1.34) gives 
v(0) = Z,i(o) — (1 + 61) '(\ — Zr) (1.39 


Expressions for (x) and i(*) may be obtained by using relations (1.38 
and (1.39) in (1.32) and (1.34) respectively. As these are somewhat 
lengthy we shall state only two which follow from (1.39). 


v(x) = e *" v(0) 
+ 3e7* —e *7)[(T + 61) (A — Zor) 
— (T — 61) “(A + Z,7)] 
(x) = ze as v(o) 
+ 4Zo'(e*" +e MIT + 61) (A — Z7) 
— 3Z0'(e*" —e "I(r — ol) (A + Zor) 
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[wo similar expressions may be obtained in which the initial current i(o) 
' instead of o(0) appears on the right. If the line is terminated in its char- 
acteristic impedance at « = 0, v(0) = —Z,i(0), and the voltages and currents 
produced by the impressed field are 

o(o) = —4 (1 + OL) (A — Zor) 

i(o) = 3 Z,'( + 61) (A — Zr) 

As in §1.4 these expressions may be computed by the first and second 
methods described in §1.3. For example, the application of the second 
method to the relation (1.39) which must exist between v(0) and i(0) in an 
infinite line gives 

m oe Z ; : 1 Z 
v(o) = 2y N(y7) |2 i(o) — h-—7 (1.42) 
r=1 Yr Yr + 6 r 
where (73) is the square matrix (1.23). 
1.6 Outline of Proofs 

The proof of the results which have been stated is divided into three parts. 

In the first part it is shown that if T’ is a matrix such that (a) its square is 


(1.41) 


ZY and (b) every element in the matrix e~* approaches zero as x > ~, 
then the expressions for o(«) and 7(x) involving I’ and Z, satisfy the trans- 
mission line equations. In the second part of the proof it is shown that if 
certain requirements are met I’ as obtained by the first method satisfies the 
conditions (a) and (b) and hence the expressions for v(x) and i(x) given by 
the first method are correct. ‘The third part of the proof discusses a general 
procedure which may be used to prove the equations which constitute the 
second method. 

Both the second and the third parts of the proof are based upon the solu- 
tion of the transmission line equations which is sketched in Appendix I. 
This solution assumes that the propagation constants of the normal modes 
of propagation are unequal, and our proofs are limited accordingly. How- 
ever, considerations of continuity seem to show that the first method is 
valid even when two or more propagation constants are equal. Under the 
same circumstances the second method suggests the use of the confluent 
form of Sylvester’s theorem.° 
1.7 Relations Obtained by Considering An Infinite Line 

We suppose that we are going to deal with transmission lines possessing 
the non-singular, symmetrical impedance and admittance matrices Z and Y. 
We further suppose that, by some means or other, we have determined a 
matrix I which satisfies the two conditions; (a) the square of I’ is 


lr = ZY, (1.43) 
. . —2zT 
and (b) every element in the matrix e ~ approaches zero as x > ~. 


°F.D.C. §3.10. 
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Consider a line extending from x = 0 tox = ©, there being no impressed 


field. Viewing the line at x = O as an m terminal network shows tha; 
there is a symmetrical matrix Z, such that v(0) = Z,i(o). Let this be 


taken as the definition of the characteristic impedance matrix Z,._ We shal 
show from the differential equations of the line that 
1. The voltages and currents in the infinite line are given by 


—eF_, 


Wx) =e § (0) 
? a (1.44 
it) = € 1(0) 
2. The matrix Z, satisfies the relations 
Z=V'Z=29"'=Try'=sy'r 
Z = 2 = PZ" = YT oe YY a 
v(x) = Z,i(x) (1.46 
3. The matrices Z,, Z, and }’ obey the commutation rules 
@(1')Z, = Z,o(T"’) 
@(T)Z = Z&(T’) (1.47 


Yo(r) = &(I’) VY 


where #(I°) is any square matrix, such as e*', representable as a 
convergent power series in I with scalar coefficients. Furthermore, 
the matrices #([)Z,, ®([)Z, and Y#(T) are symmetrical. 
The differential equations of the transmission line are 
. 2 
Cat, Sant, Da dre (1.48) 
dx dx dx? 
the third following from the first two when 7 is eliminated. That o(x) = 
e *"y(o) is a solution of the third equation may be verified by direct substi- 
tution and differentiation’. Since this expression for v(x) approaches zero 
as x —> © and reduces to 2(0) at x = 0, it represents the voltages in an 
infinite transmission line. Hence the first equation in (1.44) is true. Set- 
ting it in the first differential equation of (1.48), putting x = 0, replacing 
v(0) by Z,i(0), and noting that i(0) may be regarded as an arbitrary column 
gives 


TZ,=Z (1.49 

- 7 . oe s ol é -* 
Since I was assumed to be non-singular, Z, is equal tol’ Z. Z is sym- 
metrical and the reciprocity theorem for electrical networks requires that 7, 


7 The differentiation of the exponential function is discussed in F.D.C. §2.7. 
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be symmetrical, hence 
Zo= 1 'Z = ZI" 


[he first group of equations in (1.45) follow from this together with the 
expression IY ' for Z obtained from (1.43). The second group in (1.45) 
is obtained from the first group. 

The commutation rule for Z, is obtained from (1.49) together with the 


equation obtained from (1.49) by transposition. Since Z is symmetrical 


TS = Z.J’, Zz, =1Z,I" = ZI". 
rz. «a £y* 


and the first of equations (1.47) follow from this. The second and third 
of equations (1.47) may be obtained similarly from the relations (1.45). 
The matrix &([')Z, is symmetrical since its transposed is Z, [@(I)]’ and 
this is equal to Z,#(I’) = #(1)Z,. A similar argument applies to the 
other matrices in (1.47). 

The expression for i(x) in (1.44) may be obtained by Maclaurin’s ex- 


pansion. Setting x = 0 in the second differential equation of (1.48), 
di , a i 
(¢ = —Yor(0o) = —YZ,i(o) = —I'i(o) 
dx/}o 


where we have used the equality between the first and last members of the 
first equation of (1.45) and where the subscript 0 denotes the value of the 


derivative at x = 0. Repeated differentiation gives 


@i__y#. ae ie 
dx* dx 


a 1 492 di 93° 
(33), = I (t) = I i(o) 
Hence 


i(x) = [ _ + wi < .-- ilo 


and so on. 


Equation (1.46) may now be obtained by using the commutation rule 
for Z,: 


rT 


w(x) =e€e° oe) =e #0 7Z.i(0) 


= Ze 7" i(o) = Z,i(x) 


This completes the proof of equations (1.44) to (1.47). 
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1.8 Proof of Relations for Any Uniform Line—Impressed Field 

Here it is shown that if a matrix T satisfies the two conditions of $1.7 
and if Z, is the characteristic impedance matrix defined there, then th 
voltages and currents in any uniform line are given by the expressions (1.32 
If suitable conditions are fulfilled the relation (1.38) between v(0) and : 
for an infinite line may be obtained from (1.32). 

First of all, v(«) and i(x) reduce to the required values of v(0) and i( 
at x = 0. All that remains to be shown is that o(x) and i(x) as given by 
(1.32) are solutions of the transmission line equations (1.5). By substitut 
ing (1.32) in (1.5) and using the formulas 


: cosh sf = r'sinh aI = sinhalT 


dx 


qd. 
—- sinh af = [cosh af = cosh al T 
dx 
which follow immediately from the series definitions (1.26) of the hyper- 
bolic functions, we obtain two matrix equations corresponding to the two 
differential equations. The terms in these equations involving 2(0) may 
be canceled out provided 

I sinh af = Z sinh al” Z, 

(1.50 

, Vv yl , ’ 

I’ cosh aI” Z, = Y cosh al 
and these are seen to be true from (1.45) and (1.47). The terms involving 
i(o) may be canceled by a similar argument. The terms involving / 
may be canceled provided 


| sinh (x — £)P P(é) dé = | Z sinh (x — &)I’ Z,'U(&) dé 
0 0 


: : 
[ I’ cosh (x — &)I'’Z7" Ue) dé = [ Y cosh (x — €)I 1(E) dé 
0 J 
and these are seen to be true when «x in (1.50) is replaced by (x — &). The 
terms involving f(x) may be similarly canceled. Thus we have verified 
that 2(x) and i(x) as given by (1.32) are solutions of the transmission lin« 
equation provided that the commutation rules (1.47) and the relations 
(1.45) involving Z, of §1.7 are satisfied. This is the case when T is such 
that (a) I” is equal to ZY and also (b) every element in ¢ '* approaches 
zero aS ¥ > &, 
In order to establish equation (1.38) for the I of §1.7 several assumptions 
regarding the impressed field are required. Writing the hyperbolic fun: 
tions in the first of equations (1.32) in exponential form and premultiplying 
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both sides by 2e FS gives 

2e *' o(x) = (0) — Z,1(0) + | e (We) — Z(@) 4 
0 


+e *"[x(o) + Z,i(o)] 
at 
+e" | e 7 ©F ie) + Z,2(€)| dé 
Jo 
When x — ® equation (1.38) is obtained provided that the impressed field 
and the terminal conditions at the far end are such that (a) v(x) remains 
finite, (b) the integral in (1.38) converges, and (c), the last expression on 
the right in the equation above approaches zero as x > <. 
1.9 Derivation of Equations (1.25) 

Although equations (1.25) may be obtained by setting /(x) = tx) = 0 
in $1.8, it is of some interest to derive them directly. By repeated differ- 
entiation of the equations 

° ~#i, 2 an —Yo (1.48) 
dx dx 


the second, third and higher order derivatives may be obtained. Using 
these in Maclaurin’s expansion about x = 0 gives 


v(x) 


x . x* inetd 
[r+ 3,21 + 7 (ZY) + ++ |a( ‘ 


3 5 
~ E I+ 3, 2¥ +35 (ZY) + .-- |zito 
x ro 2, 5 iia 
i(x) = -|* I+ 31 YZ + A (YZ) + vo] Yolo) 
2 - 
+4 [ +572 +5 (vz) + .--Jivo 
These series converge for all values of x and could be used for computation 
were it not for the unfortunate fact that in most problems a great many 
terms would be required for a satisfactory answer. For the time being, 
let T be any matrix whose square is ZY. The definitions (1.26) of the 
hyperbolic functions enable us to write (1.51) as 
o(x) = cosh aT v(0) — sinh al T 'Zi(o) 
(1.52) 
i(x) = —sinh xT’ I” 'Yo(0) + cosh aT” i(o) 
If in addition to being a matrix whose square is ZY, I’ is also such that 
every element in e *" approaches zero as x — ~«, then we may use the 
relations (1.45) for Z, and obtain (1.25). 
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Incidentally, when we put ZY = Jy’ + R in (1.51) and rearrange thy 
terms so as to get a power series in R we get the series (1.28). 
1.10 Proof of the First Method 

The first method consists essentially of determining [ from the serie, 
expansion of (1.13): 


(1.53 


> —$)n (—R)” 
ao «8! >" 


b 


where (—3)n = (—3)(3)(#) --+ (w — 3) when > O and (—3), = 1, and 
then computing Z, and the required exponential and hyperbolic functions 
of xf. From §1.7 and §1.8 it follows that the first method gives the correct 
result provided that I as determined by (1.53) satisfies the conditions: 


z 


(a) its square is equal to ZY and (b) every element of e 7" approaches zero 
asx—> @. 
These two conditions are satisfied by the matrix 
Tr = PGP” (1.54 
where P and G are matrices defined by equations (A1.1) and (A1.3) of Ap 
pendix I, G being a diagonal matrix whose rth element is y,. For fron 
(A1.9) the square of T is 


r’ = PGP = ZY 
Furthermore, 


cae —o (PGP ')" 
=P> Se GP (1.95 
0 nN: 


= PM(x)P* 


where M(x) is diagonal matrix (A1.5) whose rth element is e . Since 


the real part of y, is positive and the elements of P are independent of . 
it follows that the second condition is satisfied. 

. ‘ 1 . : x? 

It will now be shown that PGP may be expanded in the series (1.53 


provided that y may be chosen so as to make all of the points ¢, = 


Y 
r = 1,2, --+ m, in the complex ¢ plane lie within that loop of the lemniscate 
|¢° — 1| = 1 which contains the point ¢ = 1. For then we may write 
the rth element in G as a convergent series: 
2 2\4 
ae. © 
> r(1 oe 
: . 
(1.56 


= (—4), (y’ — v7)" 
n=() n! yy?" 
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and PGP™ may be written as a convergent infinite series, the nth term 
of which contains the matrix (assuming only three circuits for the sake of 
simplicity) 


vi-7v 0 0 7 
| a I de P* = &, (1.57) 
0 0 ¥; ~ y 


where the equality follows from the definition (1.12) of R and equation 
A1.9) of Appendix I. This series for PGP is exactly the same as the 
series (1.53), and this completes the proof of the first method. 

The equations (1.18) and (1.28) which are incidental to the first method, 
will now be established for the case in which the matrix I occurring in 
them is equal to PGP". For then we have equation (1.55) and the 
equations 


cosh xy; 0 0 
cosh xl = P 0 cosh x72 0 Pp (1.58) 
0 0 cosh xy3 


sinh xT Z, = sinhal 'Z 


where sinh x! I~ may be expressed in the same fashion as cosh al’, the 


sinh xy, 


rth element of the diagonal matrix being The elements in the 


Yr 


diagonal matrices occurring in these expressions may be expanded in series 


, P , 2 eu fy p 
by replacing y, by its representation (1.56), assuming -"— 1 < 1, and 
7 


. oo P (__\P a 
geen = 2 (*) ( Z 4/* K,-\(z) 
0 Pp: 7 
y, < "Zz 4 1 14 
cosh 2/1 +r = »» (3) ry i T p—4(2) 
( . yA 


sinhavV/1 +r G(rz\’ 1 12 . 
eee 7G) pF ime 


where [,4(z) and K,_;(z) are Bessel functions of the first and second kinds, 
respectively, for imaginary argument. Equations (1.18) and (1.28) are 
obtained when equation (1.57) and the Bessel function recurrence relations 
are used. 


. & 
using 


a 


bo 


§ These are special cases of formulas given in “Theory of Bessel Functions,” by G. N. 
Watson, page 141. 
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1.11 Proof of the Second Method 
To establish the second method we must prove the various formulas 
which are used. These formulas all involve the square matrix N(y;) de 
fined by (1.23). 
Since N(y?) is proportional to F(y;) it follows that N(y?) may be ex 
pressed as 
N(y?) = Pip +> (1.59 
where p, is the column matrix defined in Appendix I and p, is a row matrix 
specified by p, and N(y;). Applying Sylvester’s theorem to the unit 
matrix gives 
Pi 
IT = D) N(vr) = D Pror = (pr, po, pal] pe 
_P3 
where the two matrices on the extreme right are partitioned square matrices 
From the definition of P in Appendix J it follows that 


[pr » pr, ps] = Pf, p2| =P ' (1.60 
_P3 
These relations enable us to verify the equations (1.24) when I is equal 
to PGP. Thus for the first of equations (1.24) 


Pi 
tT = PGP . = [pr ) pe; pslG p2 
P3 
Y1p1 | 
= [pi, po, pall yep | = Zz. Pr¥rPr = > N(y3)¥ 
Y3 P3 


The second equation of (1.24) follows likewise from the expression (1.55 
for e 7" 
The third equation of (1.24) follows at once from the first when we use 
(1.45), Z, = TY". The fourth equation is obtained by writing 
e*"Z, = PM(x)P'PGP'Y" 
PM(x)GP'Y" 


II 


and proceeding as in the case of the first equation. 

All of the other equations connected with the second method may be 
proved in a similar way. Incidentally, the formulas obtained by the second 
method are closely related to the “special form of solution” described in 


§6.5 of F.D.C. 
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Part II 
TRANSMISSION LINES COMPOSED OF MULTI-TERMINAL SECTIONS 


21 Introductory 

Some transmission systems may be regarded as consisting of a number 
of identical sections connected in tandem. The question of determining 
the steady state electrical behavior of such a system from a knowledge 
of the properties of a single section will be considered here. 

Each section will have a certain number, say m + 1, terminals on its 
left end and an equal number on the right. The case in which there are 
only two terminals (m = 1) has been completely worked out, and some 
studies of more general cases have been made. The ones which most nearly 
approach the point of view of the present paper are those due to S. Koizumi’. 

In the present work difference equations are used to solve the general 
case in much the same manner as they have been used in studying the 
two-terminal case. This approach differs from that used by Koizumi and 
throws additional light on the problem. 

In several lists of formulas, particularly in Appendix IV, I have included 
a number of results due to Koizumi for the sake of completeness. 

2.2 Transmission Equations for a Typical Section 

We consider the equations connecting the input and output currents and 
voltages for the mth section which is shown in Fig. 1. The directions which 
are assumed for positive current flow are indicated by arrows. The leads 
marked 0 play a special role in that all the voltages are 
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measured with respect to them, and the currents which they carry are the 
sum of the currents flowing into or out of the remaining terminals at the 
end under consideration. In applications to transmission lines the terminals 
() would correspond to the ground or the cable sheath. 

The currents and voltages shown in Fig. 1 are related by a number of 


® Archiv fiir Electrotechnik, Vol. 33, pp. 171-188, 609-622 (1939). See also a paper 
by M. G. Malti and S. E. Warschawski, Trans. A.I.E.E., Vol. 56, pp. 153-158 (1937). 
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sets of 2m linear equations which may be conveniently written in matrix 
form. One such set is 


v(m) = Zyi(n) — Zyoi(n + 1) + v°(n) 
21 
vin + 1) = Zai(n) — Zoot(n + 1) + u°(n) 


Zi, Z12, Zu, Ze are square matrices of order m whose elements are in 


I 


pedances. v(m) and i(m) are the column matrices 


V4(2) i (n) 

vo(n) io(m) 
v(n) = ° in) = ° 

Um(n) bm (1) 


The column matrices v°(m) and u°(m) arise from generators which may 
be acting within the mth section. If both ends of the section are open 
circuited so that i(m) = i(n + 1) = O the equations show that o(n 
v°(n), vm + 1) = u°(n). Consequently, v°() and u°(n) give the open 
circuit voltages produced on the left and right ends of the mth section by 
the internal generators. If the section is a passive network then v°(1) 
u°(n) = 0 and they do not appear in the equations. The subscripts on 
the square matrices, the 7’s, are chosen so as to preserve the analogy for 
the simple case m = 1, where the left and right ends of the section are 
denoted by the subscripts 1 and 2, respectively. 

Solving the equation (1.1) for i(m) and i(m + 1) gives 


i(n) = Vyo(n) + Vigo(n + 1) + 1°(n) 
—i(n + 1) = Vov(m) + Veooo(n + 1) — 7°(n) 


pet 


where the elements of the Y’s are admittances and i°(m), j°(m) are the 
currents produced by the internal sources when the terminals on the right 
and left are short-circuited so that v(m) = v(n + 1) = 0. 

A third set of equations is 


v(n) = Av(n + 1) + Bi(n + 1) — Bj*(n) 


) 3) 
i(n) = Co(n + 1) + Di(n + 1) — Cu?(n) 
Solving these equations for v(m + 1) and i(m + 1) gives 
v(m + 1) = D’v(n) — B’i(n) + B'i?(n) 
(2.4) 


i(n + 1) = —C’o(n) + A’i(n) + C’v°(n) 


There are a great many relations between the square matrices appearing 
in the equations (2.1) to (2.4). These are discussed in Appendix IV. 
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For symmetrical sections Vo = Vi2, Ye = Vu, Za = Zi2, Ze } = 
and equations (2.1) and (2.2) become 


v(m) = Zyi(n) — Zyoi(m + 1) + v°(n) 


Nm 
uw 


vn + 1) = Zyi(n) — Zyi(n + 1) + u°(n) 
i(n) = Vyv(n) + Vyoo(n + 1) + 79°(n) 
) 
2.0) 


—i(n + 1) = Vyov(n) + Vyyv(n + 1) — 7?(n) 


Eliminating i(#) from (2.5) and v(z) from (2.6) and using, from (A4.4), 


aa tay ; ; 
{ = Louden = —V\2 Vy leads to the difference equations 
v(m + 1) + v(m — 1) — 2Av(m) = Bli?(n) — 7°(n — 1)] (2.7 


i(n + 1) + i(m — 1) — 2A’i(m) = Clv°(n) — u?(n — 1)] 2.8) 
Since we also have B’ = B, C’ = C, D’ = A for symmetrical sections 
equations (2.4) become 
v(m + 1) = Av(m) — Bi(n) + Bi?(n) 
i(n + 1) = —Co(n) + A’i(n) + Cv°(n) 
We assume that the distribution of the sources in the branches of a sym- 
metrical network need not be symmetrical with respect to the two ends, 
even though the impedances of the branches are. 
2.3 Statement of Results for Infinite Symmetrical Section Line—Passive 
When the sections are passive the equations to be solved are, from (2.9), 


vin + 1) = Av(n) — Bi(n) 


(2.10) 
i(n + 1) = —Cov(n) + A’i(n) ” 
If the line extends from x = 0 ton = = the solution is 
o(n) = e "' v(0) 
i(n) = e "" i(o) (2.11) 


v(n) = Zoi(n) 


where the matrix e ' is such that (a) the equation 

e€ dhe = 2A (2.12) 
is satisfied, e” being the inverse of e ', and (b) all the elements of the matrix 
e"" approach zero as n — ©. In dealing with sections we shall never 
have occasion to consider IT itself but only its exponential and associated 
functions. The characteristic impedance matrix Z, is defined by the rela- 
tion between the initial currents and voltages in an infinite line 


v(o) = Z,t(0) 
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A formal solution of (2.12) may be obtained by writing it as 
cosh r = A (2.14 
Then 
e ' = coshT — sinhT’ 
= A — (A?— 1)! = A — (BC)! 


where the square root is to be chosen so that condition (b) for e ~ is satis- 
fied. The characteristic impedance matrix Z, is given by equations (2.34 
of which the following two are representative. 


Zo = (sinh T) 'B = sinh rc! (2.15 


: r r 
where sinh I is given by 2 sinh. =e —e 


The wide variety of sections makes it appear unlikely that there is a 
general method of determining e * analogous to the first method discussed 
for the uniform line. However, in some cases rapidly convergent series 
fore ‘ ande’ may be obtained. For example, suppose that the elements 
of (2A)"! are small compared to those of 2A. Then, from (2.12), 

e| = 2A — (2A) — (2A)-* — 2(2A) — 

e = (2A) + (2A)-* + 2(2A)-§ + - 
Again, if A* — I = BC is expressible as Iy? + R where the elements of 
R are at in comparison with y*, we have (cf. equations (1.14), (1.15)) 


r R 1/R*\ 
e -4to[1+%-3(BY +.) | 
R eT eo 
ne ee ] : (2 :) | 
| | T a aha) Tt 


Finally, it follows from a comparison of equations (2.11) and (A2.12) that 
7 “a - 
a suitable e ~ is given by 


e' = PAP", e* =QAQ"' (2.16) 


where P, Q and A are the matrices designated by the same symbols in 
Appendices II and III. 

The formal application of Sylvester’s theorem leads to a method of solving 
the symmetrical section line which is analogous to the second method 
discussed for the uniform line. Thus, if P(A) is any polynomial in A, then 


m 


P(A) = >> NG) PE) (2.17) 


r=] 
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S where P(¢,) is not a square matrix but a scalar and \V(¢,) is the square 


matrix 


N{{,) = . (2.18) 


| F(¢) is the adjoint of the characteristic matrix 


I(©) = It —A (2.19) 


and ¢1, f2, +++ &m are the roots, assumed to be unequal, of the character- 
istic equation 
f() | = 0. 


The denominator in the expression for \V(¢,) is the derivative of the char- 


i... 
f(¢,) |" — | f(r) 
si E " | t 


The formal application of Sylvester’s theorem then gives 


acteristic function: 


cosh Tr = A = IN (f,)f; 


e' = A — (A? — J)’ = SN(E,)A, 


; (2.20) 
ni TAY n yA A 
e = DN{(E,)A, 
; , 2 N(¢,) 2 
Z, = (sinh T)'B = = — B 
(A, — A;) 
where V({,) is given by (2.18), the summations run from r = 1 tor = m, 


and X, is related to ¢, through 


a 


2, =rA+X, , w= h— Ve 1 (2.21) 


where the sign of the square root is chosen so that A,| <1. A, is related 
- iy 
toe in the same way that ¢, is related to cosh I. 
24 Resulis for Any Symmetrical Section Line— Passive 
The solutions of equations (2.10) which reduce to the given values o(0), 
io) atn = Oare 
v(n) = cosh nT v(0) — sinh nT Z,7(0) 
at ‘ at unl ih se 
i(n) = —sinh mI’ Z, v(o) + cosh nI” io) 
r - : 2 am eee ; ‘ 
where e and Z, are the matrices of §2.3. These may be put in slightly 
different form by using the relations 


sinh nT Z, = Z, sinh nI” 


‘ os jd a. oe . 
sinh mI” Z, = Z, sinh nl 
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When (2.22) are interpreted by Sylvester’s theorem we obtain 


o(n) = EN (t;) | 305" + d*)o(0) — = = = BiCo) | 


r 


tm 
~m 


pe, 
i(n) = IN’(f,) |- = " Colo) + 407" + ait) | 
z ; 
where N’(¢,) is the transposed of V(¢,) and N(¢,) is given by (2.18 
the summations run from r = 1 tor = m., 
2.5 Results for Any Symmetrical Section Line—Active 

When the sections contain generators the equations to be solved ar 
those of (2.9). The solutions corresponding to the initial values o(0 
i(o) are, for xm = 1, 


v(m) = cosh nT v(0) — sinh nT Z,i(0) 
n 


+ >> {cosh (n — p)P Bi°(p — 1) — sinh (nm — p)V Z,Cv°(p — 1)} 


p=l 


= 


. ‘ oo a ® 
i(m) = —sinh nI” Z, v(0) + cosh nI” i(0) 


+ > {cosh (n — pt’ Cv°(p — 1) — sinh (n — pt’ Z,;" Bi°(p — 1)} 


p=l 


These may be simplified somewhat by replacing Z,C and Z;'B by sinh | 
and sinh I’, respectively. 

The series in the above expressions may be summed when the generators 
are such that 


v°(n) = e-"%°, i°(n) = ey? (2.25 


where @ is a scalar and 7° and v° are column matrices whose elements ar 
independent of x. Thus 


v(m) = cosh nV v(0) — sinh nT Z,i(0) 
iT Tr l/ * r ‘ 
s(e"” — e "1 )(e — e*1) (Bi? — Z,Cv°) 


Lie "* — ¢-nof)(e * — ¢*7) (Bi? + Z.Cv°) 


+ + 


‘) Ti 
(2.4U 
‘ . re ee er 
i(n) = — sinh mI” Z, v(0) + cosh nI” (0) 


+ 4(e""’ — e-*7)(e"’ — eI) "(Cv® — Z,'Bi°) 
+36" — e-nbqy(e?’ — e897) "(Cv° + Z7'Bi°) 


provided that the inverse matrices exist. 
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We may interpret these expressions by Sylvester’s theorem. For 


example, 


m - a z i = »* 
o(n) = > ne) | 40. + A*)v(0) — = Bi(o) 
r=] a vs cee 
nd > Ri 
+ td, 5 F 4 (sis -- et ) (2.27) 
2 Ay = ¢ Ar a A, 


n oa nO ) “s 
4 1A; —« (ni 4 2B )| 
2h, —e i" 


where V(¢,) is given by (2.18). 
When the line extends to » = & and the sources and end conditions 
satisfy suitable conditions we have the relation 
ie) 
oe. Trp: 7 , 
v(0) = Z,i(o) — Die” [Bi°(p — 1) — Z.Cv%(p — 1)] (2.28) 
p=l 
When the impressed field is of the form (2.25) this becomes 
r . r 1 . r ’ 
v(o) = Z,i(0) — (e — e 81) (Bi? — Z,Cv°) (2.29) 
provided that the inverse matrix exists. Expressions for v(m) and i(m) 
in such an infinite line may be obtained by using (2.28) or (2.29) in (2.24) 
or (2.26). 


Applying Sylvester’s theorem to (2.29) gives 


= 2Bi(o) Bi® Me — de 
v(o) = b> N(¢,) ( —- + = v° (2.30) 
ral My hr Ap ee”) OA, —e*) 
The last term within the braces may be replaced by 
2BCv° 
(Ar — Ar — e 


6 


) 
2.6 Derivation of the Properties of an Infinite Line 
We shall consider a symmetrical section line which is specified by the 
equations 
v(m + 1) = Av(n) — Bi(n) 
(2.10) 
i(n + 1) = —Co(n) + A’i(n) 


From these equations and the relations A? — BC = J, AB = BA’, A’C = 
CA of (A4.6) it follows that 


vin + 1) + v(m — 1) = 2Ad(n) 
i(n + 1) + a(n — 1) 


(2.31) 
2A‘i(n) 

























a Se” 


sme ge 82 












158 





BELL SYSTEM TECHNICAL JOURNAL 
= ‘ a te ile , 
Ife * isa matrix satisfying the conditions of §2.3, namely, (a) e * satisties 
the equation 
. r 
2coshl =e +e * 2A 2.14 
and (b), every element in e”' approaches zero as n—> &, and if Z, js 


defined by v(0) and i(o) for an infinite line as in (2.13), then 
1. In an infinite line 


iT 
vin) = e " v0), 
(9 29 
. nT’. o.WJe 
(in) = e 1(0), 
vin) = Z,1(n) (2.33 
2. The characteristic impedance matrix Z, is given by 
+ ° . 1 . ) l vl , . ee | 
Z, = (sinh) B= B(sinh I’) = C~ sinhI” = sinhT ¢ 
‘iP 
| - . ° i l P . la~ ws 8 ‘ am 
Z, = B sinhT = sinhTI’ B’ = (sinh I’) C = C(sinh Lr) 
3. The matrices Z,, B and C obey the commutation rules 
Py; F r’ 
®(e )\Z, = Z,P(e ) 
r ” he 
(e')B = Be&(e' ) (2.35 


C&(e') = &(e"')C 


where #(e') is a square matrix representable as a sum of powers of ¢*! 
The matrices (e") Zo, &(e") B, and C&(e") are symmetrical. 

To prove these statements we proceed as follows: By direct substitution 
into (2.31) it is seen that o(z) = e ""v(0) is a solution by virtue of condi 
tion (a) satisfied by e '. Since, by condition (b), o(2) > 0as n—> @ it fol- 
lows that v() is the voltage in an infinite line. Similarly, i(m) = e al ilo 
is the current in such a line. Substituting the expressions (2.32) for o(n 
and i(m) into the difference equations (2.10), setting » = 0, using the 
definition of Z,, and regarding v(0) and i(0) as arbitrary columns gives 

e* = A — BZ," 


(2.36 
I 


e = —CZ,+ A’ 
Applying condition (a) in the form of (2.14) to these equations gives 
re | . . _ . ’ aie 
BZ, = sinhT, CZ, = sinh I’ (2.33 


Since the sections are symmetrical, B and C are symmetrical matrices, and 
from the reciprocal theorem for networks it follows that Z, is also sym 
metrical. These remarks and (2.37) lead to (2.34). Setting the expres 
sions (2.32) for o(#) and i() in the second of the difference equations (2.10 





uSil 


y 


and 
tions 


Add 


fron 
(2.3 


the 


low 


lin 
ex] 


2e 


) 


14 


“a 
wn 











TRANSMISSION LINE EQUATIONS 159 


using the definition of Z,, and regarding i(0) as an arbitrary column gives 


(n4+1)T’ + —nI'- nT’ 
: = —Ce" Z,+ A’e 

ah nT’ ’ nT - 
(A’—e- Je =Ce" &4 


Replacing A’ —e by CZ,, as follows from the case » = 0, and pre- 


i : oe 
multiplying by C ~ gives 
, nT’ aT - 
Zoe = ¢€é be 


and this leads to the first of equations (2.35). From (2.34) and the rela- 
tions AB = BA’, CA = A'C we have 


sinh T B = B sinh I” cosh I B = B cosh I’ 
C sinh f = sinh I’ C C cosh T = cosh I’ C 
Addition and subtraction leads to 
e*'B = Be" Or ae’ ¢ 


from which the last two of equations (2.35) follow. Since each of equations 
(2.35) expresses the equality of a matrix and its transposed, it follows that 
the matrices are symmetricai. 

Equation (2.33), which is almost self-evident on physical grounds, fol- 
lows from 


/ nv al = -« 
vm) =e vlo)=e Zpt(a) 


; lr’: ae 
= Zoe" i(o) = Z,i(n). 


2.7 Proof of Relations for Any Symmetrical Section Line 

The expressions (2.24) for v(m) and 7() in a line whose sections contain 
generators may be verified to satisfy the difference equations (2.9). The 
expressions (2.34) for Z, and the commutation rules (2.35) for B and C 
are used in the verification. Setting » = 1 in the expressions for v(m) and 
i(n) gives the difference equations (2.9) and hence v(m) and i(m) are the 
solutions which correspond to the initial values 2(o) and i(o). 

In order to derive the relation (2.28) between v(0) and i(0) for an infinite 
line we put the hyperbolic functions in the expression (2.24) for v(m) in 
exponential form and multiply through by 2e "* 


de" o(m) = v(0) — Zoi(o) + Die [Bip — 1) — Z.Co%(p — 1)] 
p=l 
+ €**""[v(0) + Zoilo)] 
+t Doe Bip — 1) + Z.C0%X(p — 1)] 


p=1 












i) oe ae 


| 
















160 





BELL SYSTEM TECHNICAL JOURNAL 


Hence, letting  —>+ ~ and using condition (b) satisfied by e ‘, equatio 
(2.28) is obtained provided that (i) the terminal conditions at the far end 
are such that v(m) remains finite, (ii) the sum in (2.28) converges, and (jj 
the expression in the last line in the equation just above approaches zero 
as m—> ©, 

The results obtained by the formal application of Sylvester’s theorem 
may be verified by using the results of Appendix II and writing N(¢,) as 
the product of a column matrix and a row matrix. They may also |y 
verified more directly. For example, setting » = 0 in the expressions (2.23 
for v(m) and i(z) in any passive symmetrical section line and using 


> NV, = 7, (2.38) 


which follows from Sylvester’s theorem, we see that v(m) and i(m) reduc 
to the appropriate values v(0) and i(o) at m = 0. Substituting o() an 


ad 
i 


i() into the difference equations (2.10) and using 
BC = A? — I] 
(Ig, — A)N(v-) = N(e,)Ut, — A) = 0 
BN'(¢) = N(&®)B 
CN(E) = NSC, 
shows that they are solutions. The second of the relations (2.39) follows 


from the fact that N(f,) is proportional to the adjoint F(¢,) of f(¢,). In 
the third and fourth relations 


PO) 
IMS) | 


which is in agreement with the definition (2.18) of N(¢,). To establish 
the third relation we start from,”® 


NO) = 


($7 — A) F(g) = T\ f(g) | 
(«I — A) NY) = TOW SO |” 
Postmultiplication by B gives 
(¢1 — A) N(g) B = BI f(g) \/\ f(e) | 
We also have 
($7 — A’)F'(S) = I| f(g) | 


(I — AYN") = TSO) ASO |” 
10 F.D.C. §3.5. 
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Premultiplication by B and use of BA’ = AB gives 
(cl << A)BN (¢) = B w AC ) IMG ) 


Hence, the third equation in (2.39) holds except possibly for ¢ = ¢,, and 
from the concept of continuity it holds there also. The fourth equation in 
2.39) may be proved in the same manner. 

[he expression (2.20) for Z, may be obtained by letting m become very 
large in the expression (2.23) for v(m). (0) and i(o) must be related so that 
o(m) remains finite. Since | A,| < 1 and the X,’s are unequal the coefficients 
of A," must vanish. This requires 


IV ° 
N({,)v(0) = oe 


Ay — Ar 


Summing r from 1 to m and using (2.38) gives the required expression for Z, . 
28 The Unsymmetrical Section Line 

The method used here is analogous to those described in Appendices I and 
II for the uniform line and the symmetrical section line. The other methods 
apparently do not lead to the simplification which occurs in the symmetrical 
case. 


Equations (2.2) and (2.1) lead to the difference equations 
Vyo0( + 2) + [Var + Vos}o(m + 1) + Vor(n) = —i?(n + 1) +7°(n) (2.40 
Zi9i(n + 2) — [Zia + Zosli(n + 1) + Zaii(n) = v°(n +1) — u(m) (2.41 
Both of these equations are of the form 
Gx(n + 2) + Hx(n + 1) + G’x(n) = g(n) (2.42) 


in which G and H are square matrices of order m, H being symmetrical and 
G’ being the transposed of G. When the sections are passive equations 
(2.40) and (2.41) become 


V yov(n + 2) + [Var + Voo]lo(an + 1) + Voyo(n) = 0 (2.43 
Zit(n + 2) - [Zi1 + Z22\i(n + 1) + Zo1(n) =: ( (2.44 


In the passive, unsymmetrical case the expressions for o(#) and i(m) are 
of the form 


o(n) = PA"a + PA-"G 


tN 
rp 
un 


i(n) = OA"a — QA-G 


Comparison with (A2.8) shows that in the symmetrical case P = P and 
@ = (Q. The minus signs over P, Q, and d indicate that they are associated 


with propagation in the negative direction. The propagation constants of 
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the m modes of propagation are the same in the positive as in the negative 
direction, as indicated by the appearance of A" and A in (2.45). Cor. 
responding to any given propagation constant say A, , there are two modes of 
propagation, one in a positive direction and the other in the negative 
direction. The distribution of the voltages corresponding to these tw: 
modes are given by the rth columns in P and P, respectively. The fact 
that P and P differ shows that the distributions differ according to the 
direction of propagation even though the propagation constant is the same. 
A is still the diagonal matrix defined in (A2.3) but now the computation oj 
the elements A, is more difficult than when the section is symmetrical 
They are defined as the roots of the equation 


GM + Hx 4 G’ = 0) (2.46 


which are less than unity in absolute value. The second of the equations 
(A4.5) shows that the roots of (2.46) are the same whether the Z’s or the 
Y’s are used in place of G and H. Of course, this is to be expected on 
physical grounds. The third of the equations (A4.5) may be used to show 
that the roots of (2.46) are also the roots of 
1A — 1 AB 
= 0 (2.47 
AC AD — I 
From the form of (2.46) it follows that if \, is a root so is \;'. This fact 
may be used to simplify the determination of \,.. When the substitution 
%=AtrA', AHE—-Ve=-1 (A2.4) 
is made equation (2.46) may be written as 
0=|(G+G)%°+H+ GC -GQVe-1| 
0=|G+G)+H 


‘ : m(m — 1) : 
+ (¢ — 1) times the sum of 3 determinants each ob- 


tained by replacing two columns of |(G + G’)¢ + H | by the cor- 


responding columns of (G — G’) 
% a? . mim — 1)(m — 2)(m — 3) : 
+ (¢? — 1)* times the sum of 4 determi- 


nants each obtained by replacing four columns of |(G + G’) ¢ + 
H \by the corresponding columns of (G — G’) 


eax: 


The last equation is a polynomial of degree m in ¢ which is to be solved 


for its roots (1, £2, +++ ¢m. For simplicity we assume that these roots are 


distinct. A, is then determined from ¢, by the relations (A2.4), the sign 








of th 
his § 
alter 


W 


whe 


Let 
adj 
The 


give 


whi 
Oa 


wh 


an 
wh 


Sa 


ar 





tive 


‘or- 


Live 
tw 


fact 


me. 


1 of 








TRANSMISSION LINE EQUATIONS 163 


of the radical being chosen so that | A,| < 1 asin the symmetrical case. In 
his second paper Koizumi has given a procedure which amounts to an 
alternative method of determining A. 
We shall first assume that the Y’s are known and that our equations are 
i(n) = Vy0(n) + Yy0(n + 1) 


(2.48) 
—i(n + 1) = Vow(n) + VYoov(n + 1) 


As described above A may be computed from the determinantal equation 
| f(A) | = 0 
where f(A) represents the matrix 
f(A) = Vd? + (Vir + VYo2)dA + Va (2.49) 


Let p, be proportional to any non-zero column in F(A,) where F(A) is the 
adjoint of f(A) and let p, be proportional to any non-zero row of F(A,). 
Then the matrices P and P in the expressions (2.45) for v(m) and i(m) are 
given by 


P= [pr, po, +++ Pm 


“ (2.50) 
P = [pi, pr, +++ Dm 
where p, is the column obtained by transposing the row p.. The matrices 
Q and Q are obtained from P and P by means of the equations 
O= YuP+Y¥ePA = —YnP — YnPA™ 
(2.51) 


Q = —VyP —* VPA < => VoP + VoPA 


which are derived from (2.45) and (2.48). 
The properties of the individual columns of P and P lead to the relations 


Vi »PA? + (¥in+ Ya)PA + Y¥nP =0 
YpPA? + (Yiu + Yo2)PA* + YuP = 0 


and these guarantee that the difference equations (2.48) will be satisfied 
when the expressions (2.45) for v(m) and i(m) are used. 

When the Z2’s are known instead of the Y’s the procedure is much the 
same. The difference equations are 


v(m) = Zyi(n) — Zy2i(n + 1) 
v(n + 1) = Z21(n) — Zooil(n + 1) 
and the equation to determine the \,’s is 


|f) | =90 
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where now f(A) represents the matrix 
SO) = Zw — (Zi + Z22)A + Zar 54 


Let g, be proportional to any non-zero column in F(A,) where F(A) is the 
Oi -8 ~ of . Pa 

adjoint of f(A) and let @, be proportional to any non-zero row of F() 

The matrices Q and @ in the expressions (2.45) for v(m) and i(m) are giver 


V= Iq 25°** Omi 
Q = (91, G2, -++ Imi 


where @, is the column obtained by transposing the row q, . From (2.45 
and (2.53) equations for P and P in terms of Q and @ are obtained: 
P= Zy0—ZpQA = —Zn0 + ZnQA”" 
(2.56 
P = —-2::Q0+Z2:QA' = Zo — ZnQA 


The difference equations (2.53) are satisfied by our expressions for v() and 
i(n) by virtue of the relations 


ZV — (Zi + Z2)QA + ZnO = 0 


Nm 
va) 


ZiwQA * — (Zi + Zo) QA =o ZnQ = 0 


which are a consequence of the properties of the individual columns of 
O and Q. 

If the system extends tom = + & and if the voltages and currents are to 
remain finite at nm = the elements of d@ must be zero and the expressions 
(2.45) for v(m) and i(m) become 


o(n) = PA"a = PA"P™'v(o) 


i(n) = OA"a = OA" “i(o) (2.58 
36 ° 1 
vn) = PO “i(n), in) = OP “v(n) 


where we have assumed that P “ and Q" exist. We accordingly introduc 
the characteristic impedance and admittance matrices Z, and JY, associated 
with propagation in the positive direction, i.e., in the direction of » in 
creasing. 


vn) = Zot(n), i(m) = Y,v(n), Z.= Ys" 
Ai = PO : —— Zi — Z120AQ , = — Zo + ZoaOA ‘0 , (2.59 
Y, = OP" = Yn + VuPAP™ = —VYo — YnPA'P™ 


Incidentally, since Z, must be a symmetrical matrix the above equations 
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show that Z:0AQ™* and Zn0A 0" aresymmetrical. Z,and Y, satisfy the 
relations 

ZCZ.+ 2D —- AZ,—- B=0 Y.BY,+ Y.A — DY, -—C=0 
Zon + Zo)Zia (Zi — Zo) = Zax, = (Va2 + Vo) Vie (Yi — Yo) = Yn (2.60) 
Z0AO' = PAP"Z, Y,PAP"' = QAQ’Y, 


[he characteristic and admittance matrices Z, and Y, associated with 
propagation in the negative direction are introduced in a similar way. 


Suppose the system extends ton = — ©. Thena = oand 
vn) = PA*G@= PAP v(0) 
i(n) = QA "4 = —QA "Q“i(0) (2.61) 


o(n) = —PQi(n), i(n) = —QPv(n) 


Hence we write 


v(n) = —Z,i(n), i(n) = —Y,v(n) 
Z, = PQ? = —Zun + ZwQA'Q" = Z2 — ZunQaQ' (2.62) 
Y, = QP° = —¥u — YePA'P" = Yo + YnPAP" 


Z, and Y, satisfy the relations 
ZZ, -ZD+AZ,-B=0 Y,BY,—Y,A+ DY,—C=0 
de: an ry - (2.63) 

(Zu+Z.)Za (Zee — Z,) = Zi. (Vut ¥o)Va(V¥22 — Y.) = Vin 
The fact that Q’(Z, + Z.)Q = P'(Y, + Y.)P is a diagonal matrix may 
be used as a check on computations. 

When the expressions (2.45) for v(#) and i() are placed in (2.3), 7°(m) and 
u°(n) being zero, we obtain the relations 

PA'=AP+BO PA= AP-— BOQ 


(2.64) 
1 


QA*=CP+DO QA= —CP + DQ 


When the typical section contains generators the difference equation to 
be solved is of the form (2.42) 
Ga(n + 2) + Hx(n + 1) + G’x(n) = gin) (2.42) 
This is true for symmetrical as well as unsymmetrical sections, G being a 
symmetrical matrix in the former case so that G’ = G. The expressions for 
v(m) and i(m) are those of (2.45) with the particular solutions added: 


v(m) = PA"a + PA "G4 + u(n) 


i(m) = QA"a — QA “a + j(n) 
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where P, P, Q, Q are determined as before and u(m) and j(m) depend upon 
the generators. 
Here we shall consider only the physically important case in which the 
° ° : —né 
voltages of the generaters in the mth section are proportional to e " where 6 
is a constant. In this case g(m) may be expressed as 


nO 


g(n) = ge (2.66 


where g is a column matrix whose elements are independent of ». A 
particular solution is obtained by assuming 
-n8 
a(n) = ye" 
Setting this in (2.42) gives 
. 26 —6 7A 
(Ge + He’ + G')jy = g 
Hence a particular solution is 
\ + —20 6 . 1 n6 ~ 
x(n) = (Ge ~ + He™~ + G’) ge (2.63 
rw . ° . . —] 1 
rhis method fails when @ is equal to one of the roots Ay, -»+ Am,A1 , «++ An. 
In this case, a particular integral may be obtained by a method similar to 
one described in $5.11 of F.D.C. 


APPENDIX I 
CLASSICAL SOLUTION OF UNIFORM TRANSMISSION LINE EQUATIONS 


For the sake of convenience we again assume that there are three circuits 
in the transmission line. The equations to be solved are: 
dv di - 
-= —Zi, — = = Y§ (1.48) 
dx 


We adopt here the notation associated with equations (1.19) and (1.20), 


f(y’) being the characteristic matrix of ZY, F(7’) the adjoint of fly’), and 


vi, ¥2, ¥3(m = 3) being the roots, supposed distinct, of | f(y’) |= 0. The 
propagation constants 71 , Y2 , Ys are those square roots of “ ¥2 » 3 Which in 
physical systems have a positive real part. 

The solution may be constructed” as follows: Let the column p, be 
proportional (with any convenient constant of proportionality) to any non- 
zero column of F(y?). The non-zero columns of F (73) are proportional to 
each other according to a theorem in matrix algebra.’? Construct the 
square matrix P from the columns f1, p2, ps: 


P = (pi, p2, pal (Al. 


"The method is that described in F.D.C. §5.7(a) and §5.10 
2F.D.C. §3.5 Theorem D. 
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and obtain the square matrix Q from P: 


QO = Z'PG = YPG" (A1.2) 
where G is the diagonal matrix 
1 0 O 
G=]0 yy O (A1.3) 


The voltages and currents at any point x are 


v(x) = PM(x)a + PM(—x)a 
(A1.4) 
i(x) = OM(x)a - OM(—x)a 


where a and 4 are arbitrary column matrices associated with propagation 
in the positive and negative directions of x and M(x) is the diagonal matrix 


e™ 0 0 
M(x) = 0 eo 0 (A1.5) 
0 e 


The values of a and 4 are to be determined from the boundary conditions. 
When the line extends tox = « 


v(x) = PM(x)P™'v(0) = Zoi(x) 


R (A1.6) 
i(x) = QM(x)QO i(o) 
where the characteristic impedance matrix Zo is given by 
Zo = PO' = PG'P'Z = PGP'Y 
(A1.7) 


= Z0G'Q"' = Y'QGQ" 


Since v = p,e””* and i = q,e’”, where q, is the rth column of Q, are solutions 
the differential equations give 

(yi — ZY)p,=0, (ly — ¥Z)q, = 0 (A1.8) 
and from these it follows that 


P'zYP=0'YZ0=G@ (A1.9) 


The relations (A1.8) may be used to prove the following: 

1. The elements in the rth column of (Q are proportional to those in the 
non-zero rows of F (2). 

2. The matrix P’Q is a diagonal matrix and from this it follows that if 
is any diagonal matrix 


(PyP") = QvO" (A1.10) 


is 
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3. The characteristic impedance matrix Zp satisfies the relation 





Z=LVLo (ALM 





ran ° ° 1 1 : op , 
4. The inverse matrices P and Q “always exist if y1, v2, 3 are distinc; 





APPENDIX II 





CLASSICAL SOLUTION OF SYMMETRICAL SECTION LINE Equations —] 


The method of this section is very similar to that of Appendix I. Th 
equations to be solved are (2.10) or one of the sets 





Vie) = Z1(n) _ Zi01(n + 1) 











\2.1 


vit + 1) = Zyot(n) — Zyi(n + 1) 











i(n) = Vyyv(n) + Vyov(n + 1) 
—1(n + 1) = Vyov(n) + Viyo(n + 1) 


(A2.2 









which are obtained from (2.5) and (2.6). We shall use the notation asso- 
ciated with equation (2.19), /(¢) being the characteristic matrix of A, F(¢ 
the adjoint of f(¢), and (1, f2, +++ &m the roots, assumed unequal, of the 
characteristic equation | /(¢)| = 0. The diagonal matrices A and © are 
defined by 


where 





ze ae dy + Sa Ay = ie mis Ve —_ 1 (A2.4 


fe + V3 - | 

In general, electrical energy will be dissipated in the typical section and 
from the physical significance of \, , as seen from equations (A2.8) below, it 
follows that the sign of the radical in (A2.4) may be chosen so that | A, | < 1 
Since V/¢? — 1 = ¢ — Ay = 3A; — X,) it follows that 


> = }(A'— A) (A2.5 





Let the column matrix s, be proportional to any non-zero column in 
F(¢,) where F(¢) is the adjoint of f(¢). (It follows from the theory of 
matrices that the non-zero columns of F(¢,) differ from each other only by a 
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multiplying factor.) The matrix S is then formed by taking s; to be the 
first column, S2 the second and so on. 


S = [s,, S52, +++ Sm} (A2.6) 


. , . . . - 
Similarly let the row matrix ¢, be proportional to any nonvanishing row of 


Fic,) and form the matrix 7 where 
T = [ti , te, +++ tm] (A2.7) 


in which ¢, is the column matrix obtained by transposing ¢, .” 
Solving our difference equations for the passive case by the customary 
method gives the expressions 
vn) = PAta+ PA "a 
(A2.8) 
i(n) = OA"a — OA “4 


for the voltages and the currents. P and Q are square matrices and a and d 
are column matrices whose elements are determined by the boundary con- 
ditions. a and 4 are of the same nature as constants of integration. The 
minus sign over @ indicates that it is associated with propagation in the 
negative direction, i.e., in the direction of m decreasing. 

P and Q may be chosen in a number of ways, each choice requiring 
different values of a and d to represent the same system. In all cases, 
however, the rth column of P may be expressed as a@,s, where a, is a scalar 
multiplier which may depend upon r. Similarly the rth column of 0 may be 
expressed as §,t,. When either P or Q has been chosen the other one is 
fixed since equations (A2.2) and (A2.1) require 

O = VYuP + VPA = —Y¥unP — VPA * 


1 


P= Z0 — ZWVA = —Zin0 + ZwVA 


Some useful choices are, 


1. = § QO = —¥,,St = B’ 


SY 


2. 0 = Z25 =CS 

(A2.10) 
3. T P = ZyxT= = CTY 
4, =T P= —YusT = BT 


The particular choice to be made depends upon the system of difference 
equations which is being used. In choices 1 and 2, T is not required and in 
3 and 4, S is not required. However, if both S and T are known some of 


13 Methods of determining s, and ¢/ are available. A description will be found in 


F.D.C. §4.12, 
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the matrix multiplication may be avoided. Taking choice 1 as an exampk 
we may determine the rth row of Q from the expression 8,t,.. To determi, 
8, only one element in the rth column of — l’y2SE need be known, for 3, 
the quotient obtained by dividing this element by the corresponding elemen; 
int,. The product P’Q must be a diagonal matrix, and the same is true o; 
S’T. This may serve to check computations. 

That the expressions for v(m) and i(m) given by (A2.8) and (A2.10) satisfy 
the transmission equations (A2.1), (A2.6) and (2.10) may be verified bj 
direct substitution and use of 


S(A+ A") = 24S  T(A+A") = 24’'T (A211 
These relations follow from the properties of the individual columns | 
Sand 7. 

When the system extends to nm = © d must be zero in order that the 
voltages and currents may remain finite. This is true because X, is chose 
so that |A,| <1. From equations (A2.8) it follows that 

v(m) = PA"a = PA"P 'v(0) 
i(n) = QA"a = QA"0 “i(o) (A2.12 
o(n) = PO *i(n) i(n) = OP 'd(n) 
the reciprocal matrices Q-' and P' always exist when the sections are 
symmetrical and the roots £1, f2, +++ &m distinct. The last equations in 


(A2.12) suggest the introduction of the characteristic impedance and 
admittance matrices Z, and Y,: 


v(n) = Z,i(n), i(n) = Y,v(n), Z,= Y,'. 
Z, = PO'= Zn — Z:0A0" = —Zn + Z20A 0" 
= SES 'B = SES ‘Zu 


= Z.TIT = BTS 'T" (A2.13 
Y, =QOP¢ = ¥un+ Y¥yPAP* = —¥u — YPA'P™ 

= —V,StS = CS>'s" 

e T°TC = —TIT YY » 





aaa ra Nnca 





Not all of the expressions for Z, and VY’, obtainable from (A2.10) have been 
included in (A2.13).  Z, and Y, are symmetrical matrices. Although P and 
Q are arbitrary to some extent the same is not true of Z, and Y,. Com- 
puted values of Z, and VY, may be checked by use of the relations 


A? —] = (Z2R) = (YrY.) 
Z.CZ. = B, Y,.BY,=C (A2.14 
YZ: = —V¥ ZZ, 
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Sometimes it is desirable to terminate a line consisting of a finite number 
of sections by a network which simulates an infinite line. As is known, the 
elements in one such network may be obtained from },. Every terminal 
is joined to every other terminal, including the return terminal (denoted 
by the subscript 0), by the branches of this network. The admittance of 
the branch connecting terminal 7 to terminal 7, i # 0,7 # 0, is —y,;; where 
y,;; is the element in the ith row and jth column of V,. The admittance 
of the branch connecting terminal 7 to terminal o is via + Yio + --- + Yi 
+...» + Yim, i.€., it is the sum of all the elements whose first subscript is 7. 


APPENDIX III 
CLASSICAL SOLUTION OF SYMMETRICAL SECTION LINE EQuations—II 


When the electrical properties of a typical symmetrical section are to be 
determined by measurement, equations (A2.1) and (A2.2) show that 21 
and ¥'; may Le obtained by measurements at one end. In order to obtain 
VY 2 and Z,. measurements have fo be made at both ends. Expressions for 
v(m) and i(m) will now be given which depend only upon Z,; and J}; and 
hence are useful in case the measurements are restricted to one end. 

The method is based upon the equations 


v(n + 2) + v(m) = Zyli(n) — i(n + 2)] 
i(n + 2) + i(n) = Vulvo(n) — v(m + 2)] 


which may be derived from (A2.1) and (A2.2). Combining these equations 
leads to 


[T — ZV ujlo(m + 2) + v(m — 2)) + 2 [7 + ZyVujo(n) = 0 
[J — Vu2Zu)[i(n + 2) + ifm — 2)) + 2[7 + VuZulji(n) = 0 


The first step in the solution is to solve the equation 


| pl = ZuV |= 0 (A3.2) 
for its roots 41, “2, --+ &m Which we shall suppose are distinct. The diag- 
onal matrices M and M? are defined by 
i Bac 8 a Bee ® 
0 Me j 0 Me ; 
M =j - ; M* = (A3.3) 
De shuwas bm eer ut, 


and A is defined as in (A2.3) where X, is given by 


fj 1+] 
— M; po r 
1 = far Mr = l} san 4 (A3.4) 
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rr . ° A. e ° 

rhe sign of yw; is chosen so that | A, | < 1, and this is the value to be us 
in M*. However, there is an ambiguity in the sign of \, which is inher 
in this method. A relation between A and M? is 


M = (I+ A*)\(I — A’) (35 


. , 
Let u, be proportional to any non-zero column and w, be proportio: 
to any non-zero row of the matrix adjoint to [u,J — ZV] and form th 
matrices 


U [ui , M2, +++ Um 
W = [wy ye, tes Wm] 


(cf. equations (A2.6) and (A2.7) for S and 7) where w, is the colun 
. a , 
obtained from w,. 
The voltages and currents are given, as before, by 
vm) = PA"a + PA "a 
A28 
° n nm. 
(n) = OA a— OA a 
and there is again a number of ways in which P and Q may be chosen. | 
all cases the rth column of P may be expressed as a,u, and the rth column of 
QOas 8,w,. The equations fixing 0 when P is chosen and vice versa ar 
from equations (A3.1) 


O = Y,PM™ 


(A3.6 
P = Z,OM 


Ds ° - i 7 . > 
where M * is the inverse of M°. Equations (A3.6) may also be obtaine 
from (A2.10). 
Suitable choices for P and QO are 


1. P= U, O= Y,UM™? = ZUM 
2.P=UM', Q=Y,U = Z,;'UM 


we 
~ 


W, P= Z,WM* = Y;'wM! 
4.0=WM', P= ZW = Y;;WM 


P’O and U'W must be diagonal matrices. That the expressions for o(” 
and i(#) just derived satisfy the difference equations (A3.1) may be verified 
by making use of 


UM Zul uu, WM = VyZiW (A3.8 


Equations (A3.8) follow from the properties of the individual columns of 


U and W. The characteristic impedance and admittance matrices ar 
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sven te 
Z, = PO' = Z,0M "0" = PM'P'Y;) 
= UM*U"'Zy, = UM'UY;7 
= Z,WM ‘WwW = Y,'wM'w" 


A3.9 


Y,= OP" = YuPM'P" = OM'0"'Z;' 
= YyUM*U" = Zy'UM'U" 
= WM*W'Yy = WW 2; 
The matrices Z, and Y, may be checked by means of the relations 
Z.Y¥u = Zu¥,, ZY unZ. = Zu, VoZuVo = Vu (A3.10) 
Another set of solutions may be based upon the equations 
20(n) = —Zyli(n + 1) — i(n — 1)] 
(A3.11) 
2i(n) = Vyuolo(mn + 1) — v(m — 1)] 


which are derivable from (A2.1) and (A2.2). Combining these equations 
e . r—1 ~—1 ’ 
gives, upon using }j2Z;2 = — BC, 


v(n + 2) — 2v(m) + v(m — 2 


i(n + 2) — 2i(n) + i(n — 2) = 4CBi(n) 


= 4BCo(n) 


(A3.12) 


However, we shall not consider these equations here beyond pointing out 
that they lead to 


P = Z,A0%, OQ = —VpP> 
Px’ = BCP, Ox’ = CBO 


which may also be derived from (A2.10). 


APPENDIX IV 
RELATIONS BETWEEN THE SQUARE MATRICES OF 
A MULTI-TERMINAL SECTION 
When the reciprocal theorems of network theory are applied to equations 


(2.1) and (2.2) it is found that Zy, , Z22 


22 » 


Zn = Zi:; Vou = Vis (A4.1) 


Yui, V22 are symmetrical and 


i.e., Zo: and Vo: are the transposed matrices of Z;. and Vy», respectively. 
Solving equations (2.1) for the currents and comparing the result with 
(2.2) shows that 


Zu Zi in o Yu Yi2 i°(n) ae oS Yu Yie v°(n) 
Zu Zn ‘7 Yo Vx ‘ —j°(n) = Ya Yx u°(n) 
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These are partitioned matrices Ihe square matrices have 2m rows ay 


columns and the column matrices have 2m elements Phe first of the 
relations may be written as 


Sa Batti Vu Vu I ol 


Zn Zn\\¥n Yn) {0 1| \ 


L 
where / denotes the unit diagonal] matrix of order m. Multiplying the 
matrices on the left together and equating the elements of the pr. 


to the elements on the right gives 
; Met JL Ui Tipit PIy 


SN NN 


lransposing the matrices in these equations leads to other relations 


from the first we obtain V4Z1; + Vi2Ze) J These equations also vic 


f ' F y 
eynress< y forthe?) s mn termes of the ZAsand Vice vers: 
4 somewhat sin r treatment NVOlViIng e« ations Z ] ne Zo ¢ 
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Zi: NZ 4 MZa + ZadlYa + Vu 
(AZ i2 MPa + Aba + Al t Pay A4.5 
Vn 0 \t—] \i Vag + Vu \l 
0 221 AC \p { \/ AZ fas | 
Sometimes it is of interest to obtain the elements of |, ay. when 
Ary, Za, } il» } yy are known Relations helpful ih studying this problems 
ure 
TV ATO by } bay alia tii ! i} i 
Vi9)'aa } i bia ‘W ig iit ag . 
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ion is symmetrical some Siinipolite dlioi takes jp bet 


Vis Vyy AY 42 ! YY 1h bA 
Vis Vai Zi 421 hi fs’ 1 (4 
( ( 1° — he / A4.6 
Aiviait ZisVis / (‘RBA ( Is 
AiuViart Zab 0 
APPENDIX \ 
Provertivs or THE Marrix GA? + I) (,' 
The matrix 
fOr GX + Hr + G’ | 
which entered the discussion of the case of unsymmetrical sections has 
imber of interesting properties which are given belo G al Hf are 
square matrices with m rows each, and // is required to be symmet: 
4s before, we shall denote by A; Ae, ri A, the 2 wm roots ol ‘ 
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These are partitioned matrices. The square matrices have 2m rows and 
columns and the column matrices have 2m elements. The first of thes 
relations may be written as 


Zi Zi Yu Yio = I 0 (A ; 
Zn Ze\| Vn Yo} |0 I M42 


where / denotes the unit diagonal matrix of order m. Multiplying the tw 
matrices on the left together and equating the elements of the product 
to the elements on the right gives 


ZuVu t+ ZV 
ZuVi2+ Zi2Vo = 0 


II 
— 


(A4.3 
ZaVir + ZeVo = 0 


ZuVi2 + Z2V22 = 1 
Transposing the matrices in these equations leads to other relations. Thus, 
from the first we obtain 1121, + Yi2Z21 = J. These equations also yield 
expressions for the ]’’s in terms of the Z2’s and vice versa. 
A somewhat similar treatment involving equations (2.1) and (2.3) leads 
to expressions for the Z’s in terms of A, B,C and D. The Y’s may be like- 
wise expressed. These relations are given in the following table. 


VY; = DB" io = 2a - Zaks. 
Vie= C — DB'A = —B’ Vit = Zn — ZeZeZu 
Yn= —B" Va = Z2 — ZuZnZe 
Vo = B ‘A V2 = LZ — ZnZu Zi 
Zu = AC" Zn =Yu- V2 22 Vou 


Zax AC'D—BzC" Ze w= Yu — Ya¥ul¥ un 


Zn = € Zu = VYy- Vuln Yo 
Z2=C'D Zu = Vo — Yalu V2 
(A4.4 
A=Zy2Zn = —YuaYn 
B= Zy2ZnZy— Zp. = —VYn 
C = Zn = Ve —- Vu¥un¥ 2 
D = ZaZy= —Vu¥un 
A B/'_[D' -—-B’ AD’ — BC' = 1 AB’ = BA’ 
C i: Lt CD’ = DC’ DA’ — CB' =1 


“On A < hn 
—i(n)} LC DJL—j\(n) 
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3 


| The following equations in which ) is an arbitrary scalar multiplier may 
be verified by equating coefficients of powers of \ and using the relations 


* just given. 


Zu — A“%u)(Viur + AV a2) = (AZi2 — Zo2)(AVo2 + Yo) 
(Zig — AZ — AZo2 + Zar)( Yur + AV w) 
= (AZ. — Zo2)(MV a2 + AV + AV 22 + Vos) (A4.5) 


—V¥, O }[A4A —I \B _ [AVa2 + Yor MI 
0 Zu AC AD —-I| MI NZ — Zu 


Sometimes it is of interest to obtain the elements of Vi2, say, when 
Zi, Z2, Yi, Ye2 are known. Relations helpful in studying this problem 
are 

VuZuVie = Vi2Ze2V 2, PuZulu — Yu = VileVo 
YeVe¥n=V¥n-Zn Ze = —ZuVaVe 
YalunVie= VYa—- Ze Zn = —Y¥a(VnZn— 7) 
When the typical section is symmetrical some simplification takes place 
and we have 
Vu a Yo Zi) = Zo A = D’ AB = BA’ 
V1 = Yo Zi2 = Za B = B’ AC = CA 
C=C’ A?—BC = I (A4.6) 
Zu¥u t+ Z2¥i2n = I A'B"A—C=B" 


ZuViet+ ZeV¥u = 0 


APPENDIX V 
PROPERTIES OF THE MATRIX GA? + HA + GC’ 
The matrix 


f+) = GV + A+ G' (A5.1) 


which entered the discussion of the case of unsymmetrical sections has a 
number of interesting properties which are. given below. G and H are 
square matrices with m rows each, and H is required to be symmetrical. 
As before, we shall denote by \1, --- Am, Al, +--+ Am_ the 2 m roots of the 
determinantal equation 
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and we shall suppose these roots to be distinct. Let the column &, and the 
row /, be such that 


k,l, = F(A,) (A5.2) 


where F(A) is the matrix adjoint to f(A), and let the square matrices K and 
L be defined by 


h 
ly 


K = [ki, ko, +++ Rul, L= (AS.3 


Im 
Comparison of (A5.3) and (2.50) suggests that when G and H are expressed 
in terms of the Y’s we have the relations 


K=P, L=FP' (A5.4 


The method of choosing the column p, and the row p, shows that they are 
related by 


bed, _ yrF (Xr) 


instead of (A5.2) where y, may turn out to be any non-zero constant, and 
consequently equations (A5.4) are not satisfied in general. Nevertheless A 
and L may be regarded as particular choices for Pand P’. In the same way 
K and L may be regarded as particular choices for Q and Q’ when G and // 
are expressed in terms of the Z’s. There is still some arbitrariness connected 
with K and L since the product ,/, is unchanged when the &, is multiplied 
by some number and /, is divided by the same number. 
The relations which correspond to (2.52) and (2.57) are 


GK’ + HKA + G’K = 0 


(A5.5 
GUA* + BUA ' +G'L' =0 


where A is the diagonal matrix whose elements are \;, \2, +++ Am. These 


relations are consequences of the properties of k, and /,. Two more rela- 


tions may be obtained by transposition. From the first of (A5.5) and the 


transposed of the second it follows that 


GKAK' + H+GKA'K™ = 0 
: (A5.6 
L*sALlG + 822s" = 0 


. 2 ° P a 1 . . 
where it is assumed that the reciprocal matrices K and L~ exist. Com- 


. . — . | ye eee ‘ » 
binations similar to KAK ', KA K , etc. enter the expressions (2.59 
for Z, and Y,. 
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By differentiating the equation 


fF) = AA), 


‘where A(A) is the determinant 


A(A) = | f(A) | = |G + AAFC 
t may be proved that 


Gear 40 + CALG = LRK 


ee 2 (AS.7 
G’KA'K '+H+L°A'LG@ = —L"'EK™ 
in which £& is the diagonal matrix 
r (1) 
| AQa) 0 a at 0 
} (1) 
0 Arg) 
E= 
0 1) 
A(Am) 
and 
. d 
AQ,) = A(\) 
dX = 
1 
Since the roots A, are assumed to be distinct, A(A,) + 0. 
We also have the equations 
KE'L = L'E'K' 
A5.8 


GKAE'L — GLU’A ‘EK’ = I 


rhe first equation of (A5.8) shows that AE L is a symmetrical matrix. 
From this and the second equation it follows that 


‘E pil ‘yr lyr-l Lypyp-l = 

GKAK © — GL'A  L LEK A5.9 
From the first of equations (A5.7) and the second of (A5.6 

. l 1 v lp 1 . 

GKAK~ — L'A'LG' = LEK A5.10 
. ° ~ . , rl - 

and the comparison with (A5.9) shows that the matrix GL’A ‘L’™ is sym- 
metrical. The other matrices of this type are also symmetrical as may now 
be seen from equations (A5.6) and (A5.7). Results of this sort may be 


obtained from physical principles by noting that Z, and VY, must be sym- 
metrical matrices. 
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As an example of the application of these formulas we assume that ( and 

H are expressed in terms of the Y’s. Then we may take K and L’ to be 

I ) € 
particular choices of P and P and equation (A5.9) becomes 


ViePAP' — VpPA'P* = (P’)"EP™, 
From equations (2.59) and (2.62) 
¥o+ ¥, = YuPAP — YyPa'P' 
and hence 
P'(Y,+ ¥,)P = E. 


For the more general choice of P and P allowed in §2.8 the diagonal matrix 


é 


£ is replaced by a general diagonal matrix. Similarly it follows that 


Q'(Z. + Z.)O 


is a diagonal matrix. 
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Engineering Problems in Dimensions and Tolerances 
By W. W. WERRING 


DIMENSIONAL UNITS 


The basic unit in most considerations of dimensions in the United States 
is the inch. The value of the inch is so important that many companies 
including the Bell System maintain in their measurement laboratories a 
standard yard bar calibrated against the standard at the National Bureau 
of Standards. In spite of this it is an interesting and curious fact that 
though all have been much concerned over the legal value of the dollar 
there has been little interest even among engineers in the exact legal value 
of the inch. Actually there is no single answer to so simple a question as 
“What isan inch?” In fact, we have changed from a British inch and our 
own legal meter, to our inch and the International meter and now through 
action of the American Standards Association we are actually using an inch 
based on conversion from the International meter which is neither our 
own legal inch or the British legal inch— and the British are using it too. 
Table I shows this history of the legal inch in the United States. 

It will be seen that under the present status there exists a difference of 
two parts in a million between the legal inch and the inch used in the di- 
mensional work of industry. This difference is more theoretical than real 
in small dimensions and industrial use. The bill before Congress, sponsored 
by the Bureau of Standards is intended to eliminate this as well as any 
possible ambiguity in the U. S. inch. 


DECIMAL DIMENSIONING 


In subdividing the inch the modern trend in industry is toward the use of 
decimals instead of the older common fractions although fractions continue 
to be used, especially for dimensions of certain materials such as iron pipe, 
lumber, phenol fiber. In fact even a special decimal system based on using 
only the tenths and fiftieths of an inch is being considerably discussed by 
general industry. This system would use a scale on. which the smallest 
division is sy” or .020” instead of gy” = .0156”. It is in use by the Ford 
Motor Company and the values shown in Table II are some of those used 
in place of common fractions. The decimal equivalents of these common 
fractions are also shown rounded to 3 decimal places in accordance with 
American Standard Rules for Rounding off Numerical Values Z25.1-1940. 

In the Ford system one and two-digit decimals carry the general toler- 
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ance of + .010”. When greater accuracy is required three-place decimals 
are used to express a minimum and a maximum value. 

The adoption of decimal dimensioning for all drawings prepared at Bell 
Telephone Laboratories is being actively considered. However, adoptio 


TABLE I 
History OF UNITED STATES DIMENSIONAL STANDARDS 


Year Action | Resulting Dimensional Relationships 
| _ 

1830-36 Adoption for Customs Service and for distribution to individual states oj 
standards intended to be the English yard based on a certain portion of 
an 82 inch bar imported in 1813. The portion selected was supposed 
to be identical with the English yard. 


1856 Official copy of new British Imperial International Meter = 39.370147 
Yard accepted as standard British Inch 

1866 Congress declared metric units law- Legal Meter in = 39.37 
ful and established legal equiva- u, 8. British Inch 
lents 

1893 Mendenhall Order set up Inter- | International Meter = 39.37 
national meter as the fundamental U. S. Inch 
standard 

1933 American Standards Association | International Meter = 39.370078 
(Representing Industry) adopts U.S. Inch 


1 inch = 2.54 centimeters 
1937-41 Bill before Congress but held in International Meter = 39. 370078 


committee for amendments U.S. Inch 


TABLE II 
EXAMPLES OF ForD DECIMALS COMPARED TO COMMON FRACTIONS 


: _o i Ste d 
Decimal of Existing American Standar 


Ford Decimal Common Fraction Commnon Fraction an — its 
( ace) 
02 1/64 -015625 .016 
03 1/32 -03125 .031 
05 3/64 -046875 | .047 
.06 1/16 .0625 .062 
.08 5/64 -078125 078 
.3 7/3 .21875 219 
46 | 15/32 46875 | -469 


of decimal dimensioning would not of itself result in any changes in our 
system for establishing tolerance values. 


Raw MATERIAL SIZES 


In contrast to this continued trend toward simplification and rationaliza- 
tion of our systems of dimensional units raw material supply is still com- 
plicated by a multitude of obsolete systems of gauge sizes in every day use. 
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Many in industry have probably grown used to the standard gauges in 
particular fields but though gauge numbers were undoubtedly initiated as 
a simplified identification the variety of gauges and the variety of names 
for the same gauge now merely increases confusion. Sheet metals are 
handled in terms of a number of gauges such as B&S gauge, U. S. standard 


gauge and BWG gauge; and sheet soft rubber is even designated in decimals 
of dy such as $?”. It has become good practice to specify sizes by decimal 
dimension values and not by gauge numbers and holes by actual decimal 
size rather than by drill numbers. The actual sizes used, however, are 
determined in many cases by the values corresponding to old gauge num- 
bers long used commercially, though in large running items mills will and 
do manufacture to any specified decimal size. For some time it has been 
the practice of material manufacturers and other large industries thus to 
discontinue the use of gauge numbers though still using the decimal values 
of gauge sizes. 

There is now under way an effort, organized under committee B32 of the 
American Standards Association, to eliminate the old wire and sheet metal 
gauge systems entirely and set up a rational series of American standard 
thicknesses for all metal sheets and preferred diameters for wire, and insure 
availability in these sizes. The basic conception of a rational series of 
sizes is that a uniform degree of choice should be presented between suc- 
cessive sizes. Therefore each size should differ from the next by a fixed 
percentage. The series should therefore be geometric. A variety of geo- 
metric series could be used but in order to permit extending the series in- 
definitely by shifting the decimal point, the particular series based on the 
root of 10 has been established internationally as the Preferred Numbers 
Series for standard sizes. The 5 series is one having 5 numbers between 1 
and 10 (or between 10 and 100) and is produced by using as the multiplier 
the fifth root of 10; the 10 series is produced by multiplying by the 10th root 
of 10; the 20 series by multiplying by the 20th root of 10 etc. The complete 
Preferred Numbers Series is explained and listed in various forms in Ameri- 
can Standard Z17.1-1936. 

The subcommittee working on the sheet metal sizes has recently issued a 
proposed American Standard of preferred thicknesses for all uncoated flat 
metals thinner than .250”. These thicknesses are all decimals based on the 
20 series of preferred numbers rounded in the standard manner to 3 decimal 
places. The Preferred Numbers and the proposed thicknesses are shown by 
Table III. It happens that this series closely approximates the Brown and 
Sharp gauge used in the nonferrous metals which simplifies that portion of 
the changeover. If this proposed American Standard is generally approved, 
as now appears most promising, we will be able to choose thicknesses of any 
metal interchangeably without the restrictions of ancient gauge sizes es- 





ez 


yrange soe 
pc Teed 


bax. 


















BELL SYSTEM TECHNICAL JOURNAL 


TABLE III 


Decimal Series of Preferred Numbers 10-100 Proposed Preferred American Standard Thick 
5 Series 10 Series 20 Series . 
5 10 . 20 Under .010 -010 to .100 11201 
V 10 = 1.6 V 10 = 1,25 V/10 = 1.12 

10 10 10 .010 

it .2 O11 112 

12.5 12.5 .012 125 

14 .014 140 

16 16 16 .016 . 160 

18 018 . 180 

20 20 .020 200 

22.4 .022 224 
25 25 25 .025 
28 .028 
i 31.5 .032 
35.5 .036 
40 40 40 .004 .040 
45 045 
50 50 .005 .050 
56 -056 
63 63 63 .006 .063 
71 .007 071 
80 80 .008 -080 
90 .009 .090 


100 100 100 . 100 


tablished for reasons which were possibly good and sufficient but which 
certainly have long been forgotten. Meanwhile, another subcommittee is 
investigating the possibility of applying a similar series to the diameters of 
wire. Probably diameters to 4 decimal places will be required. 


DIMENSIONAL TOLERANCES 


Part Tolerances 


Regardless of the dimension decided upon in a design it is obvious that it 
cannot be regularly manufactured to the exact size. Certain manufac- 
turing variations or tolerances must be expected and these introduce a large 
share of our dimensional problems. 

The usual statement on tolerances is that the larger the tolerance 
allowed the cheaper the part is to manufacture and, therefore, the tolerance 
specified should be the widest that will permit functioning. However, 
this is generally true only of overall tolerances which define the manufac 
turing methods that may be used. It is true in the sense that apparatus is 
inexpensive to manufacture if it can be so designed that its functioning i 
largely independent of variations in dimensions. However, such design 
not usually achieved and in much apparatus fairly good overall accuracy ‘ 
dimensions and fit is necessary for uniform functioning. 
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PROBLEMS IN DIMENSIONS AND TOLERANCES 
setting tolerances then becomes one of distributing certain tolerances over 
various dimensions and different parts. This is a very difficult problem 
and in the case of any individual tolerance a larger value does not neces- 
sarily mean lower apparatus cost and may even mean the reverse. 
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STOCK DIAMETER IN INCHES 


Fig. 1—Total diameter tolerances of commercial round stock 


This is easily demonstrated in the case of part tolerances on dimensions 
which correspond to the dimensions of raw materials. Figure 1 shows the 
tolerances of commercial grades of round stock. If, for example, engineer- 
ing requirements dictate the use of a particular material there is no gain in 
specifying larger tolerances than those to which it is regularly furnished and 
doing so may require greater accuracy in the mating part. There may even 
be economy in the use of higher priced material produced to closer toler- 
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ances, as for example, drill rod instead of cold drawn steel through economy 
in the manufacture of associated parts. Similarly manufacture of cantileve; 
springs from sheet stock produced to closer tolerances may reduce the cos; 
of subsequent adjustments. Therefore, when individual part tolerances are 
involved consideration must always be given to the size tolerances of ray 
materials. 

The same situation exists in the case of tolerances on dimensions pro- 
duced by a manufacturer’s own tools. While close overall limits will re. 
quire greater overall accuracy of the tools provided and greater frequency 
of set-ups the most economical distribution of tolerances will be that based 
upon the normal tolerances that can be expected from various manufactur- 
ing operations. Certain degrees of accuracy are inherent in certain types 
of machines and tools and allowing variations not in proportion to these 
values serves little if any purpose. Also there are types of combination 
tools and automatic machines, familiar in mass production practice with 
which wide tolerances are not an economy because accuracy is required for 
locating or nesting the part for subsequent operations. Since the dis 
tribution of tolerances involves such complex factors of manufacturing 
method and cost as these, it is desirable for the designing engineer to de- 
termine and to indicate unmistakably the effect of tolerances upon func- 
tioning and, where interchangeability of individual parts in service is not 
involved, to allow manufacturing considerations to determine the distribu 
tion of tolerances in an assembly. 

It is apparent that considerable study of the requirements for functioning 
of the design, of available materials and the limitations of manufacturing 
process are required to establish the most economic balance between per- 
formance of the apparatus and the required tolerances. Consideration 
should be given to these tolerance factors in cooperation with manufac- 
turing engineers in an early stage of a design problem so that they may 
influence the trend of design. This step may avoid the necessity for slow 
and costly manufacturing developments and delays in starting production. 
However, completely rigid adherence to the status quo of tolerances is not 
necessary in long range planning of major design projects. In such cases 
the trend of progress in materials and manufacture should be determined 
and anticipated. For example, some cantilever spring design requiring 
narrow control has been based on sheet material produced to tolerances 
not commercially available at the time but made so by the time it was 
needed for production. The extent of progress in this direction is shown 
by Fig. 2. 

Similar progress in manufacturing technique can also be expected. For 
example, the development of broaching from a comparatively crude opera- 
tion to the precision method it is today is recent and outstanding. 
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Cumulative Assembly Tolerances 


Another problem in choice of tolerances is in those cases where a \ 
siderable number of parts are additively assembled into a unit as in the cas 
of “spring pileups’’ used on electrical contacting apparatus such as relays 
and switches. These consist of considerable numbers of sheet metal spr ngs 
and insulators alternating and clamped by screws. If the overall tole: 
ance on such an assembly must be taken as the sum of the tolerances of th, 
individual parts various courses of action are presented, the extremes «| 
which are: 

1. Very small tolerances must be maintained on the individual parts or 

2. Adequate space must be provided in the apparatus for extreme 

large variations in the assembly. 

Small tolerances on the individual parts may be extremely expensive a1 
large space allowances and provisions in associated parts for variations ii 
the assembly may be a serious design handicap. 

However, it is recognized that there is obviously small probability that 
all minimum or all maximum parts will appear in any one assembly. I 
has been found satisfactory in certain types of such pileups to assume that 
the maximum dimensional variation that will actually be encountered in an 
assembly will not be greater than 70% of the sum of the part tolerances. A 
similar situation exists in many kinds of assemblies or associations of toler 
ances. 

The statistical relationships involved in this problem are indicated by 
Fig. 3. The curves show the percentage of the cumulative part tolerances 
within which 99.7% of the assemblies may be expected to be found wit! 
different numbers of similar units in the assembly. The solid line is de 
duced from theoretical relationships. It assumes that the parts are all of 
one kind, that the parts going to assembly are controlled, of normal dis 
tribution and the limits are rationally set to represent the actual conditions 
The dotted curves have been deduced from relationships which have been 
proposed as representing rectangular and triangular distributions of in 
dividual part tolerances. The curves may not be truly representative o! 
specific cases because of inconsistent selection of limits or erratic distribu- 
tions. However, they indicate that the 70% rule on pileups is probably on 
the safe side in most cases and that closer design of assembly or less re 
strictive tolerances and cheaper manufacture of piece-parts might be readily 
possible either (1) by better control, (2) by actual mixing of lots of piece 
parts or (3) even merely by knowledge of the actual statistical distribution 
of part dimensions. 

The three points indicated in Fig. 3 show the results of a limited experi 
ment’ in which pileups were assembled from 2083 individual insulators 0! 


37” phenol fiber taken from factory stock. The establishment of curves 
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by this type of experiment using a sufficiently large and representative 
sample would be practicable and would permit considerable condensation 
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NUMBER OF UNITS IN ASSEMBLY (ALL ONE KIND) 


Fig. 3—Statistical relationship of overall tolerance on an assembly and the sum of the 
individual part tolerances 
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THICKNESS AT CENTER OF CELL IN THOUSANDTHS OF AN INCH 
Fig. 4—Distribution of thickness in 2083 pileup insulators 
of design on a sound basis. In this particular case the parts used apparently 


came from only two different sheets of fiber as indicated by the distribution 
of thickness of the individual parts shown by Fig. 4. 
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Further statistical analysis of this type of situation is needed together 
with experimental determination of the distribution of dimensions actually 
encountered in specific cases 

The distribution of dimensions in a product arises from a variety of causes, 
One type of cause is the variations such as those in the dimensions and 
physical properties of raw material which may produce different product 
dimensions even from a particular tool. A different and more systemati 
type of cause is the change in the dimensions of tools as a result of wear 
The practice followed in establishing tool wear allowances will therefore 
affect the limits and statistical distribution of part dimensions during the 
life of the tool. Some designers and some tool makers consider that the 
specification of a nominal value with plus and minus variations requires a 
different handling of initial tool dimensions and wear allowances than does 
the specification only of minimum and maximum limits for a part dimen- 
sion. Equally good authority maintains that a manufacturer recognizes 
no difference. Establishment of standard practices in such matters is a 
needed step in determining the distribution of dimensions to be expected in 
machined parts. In the present absence of standards or of any consistent 
attitude on the subject it is necessary for designing and manufacturing 
engineers to reach an agreement in specific cases where this factor is im- 
portant. 

Such are the factors which determine the tolerances which can be obtained 
economically or which perhaps will be unavoidably encountered. It is 
necessary for a designer to keep informed of the interaction of these factors 
as his design crystallizes and he must also determine the effect of such toler- 
ances upon functioning in order to complete a design which will function 
properly when assembled in quantity production. 


FUNCTIONAL DIMENSIONING 
Effect of Tolerances 


If apparatus parts are minute or have complicated relative motions it is 
recognized that manufacturing drawings to the usual scale have serious 
limitations to their usefulness in the analysis of the effects of combinations of 
tolerances. In such cases designers frequently make layouts to larger 
scales or large scale adjustable models to investigate the effect of variations 
on functioning. Illustrations of this practice are numerous in the expe- 
rience of most designers of small apparatus. 

Even in large parts which are stationary in use the application of toler 
ances, in effect, establishes several possible positions for each element and 
may present problems similar to those involving motion. These are not 
easily recognized because of a curious limitation inherent in small scale 
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drawings. This limitation is probably well known to most engineers but 
it is worthwhile to analyze it because it is important to be always aware 
of it. 

This limitation is the fact that in drawings the shape of the part and the 
effect of all nominal dimensions are actually shown graphically whereas, 
it is possible to indicate tolerances numerically but not graphically. We 
are therefore apt to visualize the part as it is graphically shown, that is, 
without tolerances and to think of the numerical tolerances one at a time 
rather than in combinations as they affect each other and the shape of the 
part. 

If any dimension, significantly affecting the design of a part, is changed 
the drawing is immediately corrected so that its meaning will be clear and 
the functioning of the part can be checked. This obviously facilitates de- 
sign and manufacture. Yet because they cannot be shown directly by 
regular drawing methods, we have grown accustomed to not being shown the 
effect of tolerances or changes in tolerances upon the shape of the part. 
Nevertheless it is obvious that these effects are critical in the functioning 
of the part or tolerances would not be set. The fact that these critical 
features of the design are not actually graphically shown and therefore are 
not easily seen and understood on the drafting board is a serious detriment 
in working out a design and in all later analysis of it. The full effect of 
interrelated variations particularly if in three dimensional space may ap- 
pear only after tools are in process or the first parts produced and this may 
be rather late for economy. 

Originally this difficult analysis of the effect of tolerances upon function- 
ing probably involved only the designer. The manufacturer tried to make 
the part as nearly as possible to the nominal values shown and variations 
from them were accidental. Tolerances were looked upon as an indication 
of the care required and as a means of inspection for acceptance or rejec- 
tion. With increasingly complex manufacturing tools the permitted toler- 
ances are utilized more and more in the design of tools to allow the greatest 
possible wear before defective parts are produced and the tools must be 
replaced. For mass production parts progressive step type tools are used 
in which a continuous strip of stock advances by various stages from blank 
sheet tofinished part. Tools ofthis typeare extremely expensive and in order 
to obtain maximum life full use of allowed variations is made in their de- 
sign. Design of such tools and the gauges required to maintain quality in 
mass production therefore also requires analysis of the effect of combina- 
tions of variables upon the desired part. As the designer has presumably 
already made this analysis, and incidentally is best qualified to do it, econ- 
omy and accuracy dictate that his analysis be transmitted to the manu- 


facturing engineer. The problem is to find means by which he can indicate 
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unmistakably on the drawing his analysis of the required functioning of thy 
part and the manner in which he intends the tolerances to apply, in the 
event that there is any possibility of misunderstanding. 

The essence of this problem and some of the possibilities of solution ca; 
best be seen by reference to drawings which illustrate the major points. 

Figure 5 shows the drawing of a flat plate dimensioned from center lines 
but without any tolerances whatever. Some minor dimensions not in 
volved in this discussion are omitted in the interest of simplification but 
the part shown is in every way a normal one. The meaning of the drawing 
is completely clear and can be interpreted in but one way no matter fron 
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Fig. 5—Flat plate dimensioned without tolerances 


what standpoint the analysis is made. The reason for this is obviously 
that but one value is shown for every dimension. 

Figure 6 shows this same drawing dimensioned in exactly the same way 
with the exception that tolerances are shown for most of the dimensions. 
To the uninitiated it might appear to present no more problem than the 
previous drawing without tolerances because of the tendency to visualize 
the drawing in terms of the nominal dimensions only. 

When the engineer analyzes the effect of the combinations of the various 
tolerances shown, interesting questions immediately arise. In the first 
place the combination of holes dimensioned 1.25” + .002” from the center 
line appears to be definitely located because on the drawing the center 
line is shown in a definite position. Yet when the tolerances are considered 
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the center line of this drawing could actually be shown in several different 
places as, for example: 
1. It may bea line through the centers of the two large holes. 
2. It may be a line anywhere from 2.992” to 3.008” from the outside 
edges. 
3. It may be 2.247” to 2.253” from the small holes in the center of the 
plate. 
4. It may be 2.615” to 2.635” from the holes numbered 2 and 4. 
In brief, the center line which appears so definitely located on the drawing 
may actually be rather an indefinite location on the part when the various 
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Fig. 6—Flat plate of Fig. 5 with the addition of tolerances 


tolerances are considered. While the differences in the possible interpreta- 
tions are in the order of thousandths of an inch nevertheless this order of 
magnitude is critical in this part or the indicated tolerances would not have 
been used. The interpretation of the center line which should be adopted 
will depend entirely upon the manner in which the part is intended to func- 
tion and therefore should be indicated by the designing engineer. Obvi- 
ously, not all designs or all dimensioning will present this difficulty but all 
should be studied from this viewpoint to determine whether or not they do. 


Functional Datum Positions 


When the type of uncertainty illustrated exists, it is necessary to 
indicate clearly the effect of tolerances on functioning by establishing 
the functional positions to which dimensions should refer. It may be 
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difficult to do this graphically, in which case it is necessary to indicate by 
notes the particular interpretation which the designer intends. As ay 
example, if the part of Figs. 5 and 6 functions by being located in position 
by means of the four holes numbered 1, 2, 3 and 4, the intentions of the 


designer are readily indicated by the following notes: 
1. Functional datum line I is midway between the centers of holes 1 and 
2 and the centers of holes 3 and 4. 
2. Functional datum line II is perpendicular to datum line I at a point 
midway between the centers of holes 2 and 4. 
These notes establish both horizontal and vertical center lines specifically 
in terms of the center of the one set of dimensions between the holes marked 
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Fig. 7—Flat plate of Fig. 5 functionally dimensioned from outside edges with tolerances 


1, 2,3 and 4. The term functional datum line is suggested as completely 
descriptive but other equivalent terms might be used. This information 
could be indicated on the drawing without the use of notes by the adoption 
and use of some standard convention or symbol to indicate the particular 
dimension bisected by the center line. 

If the functioning of this part were determined by location against the 
outside edges, this could be readily indicated by dimensioning the part as 
shown by Fig. 7 and using notes establishing the line A~B as one datum 
line and the perpendicular to it through A as the other. 

In either of these cases the drawing becomes completely definite and sub- 
ject to only one interpretation. In drawings of this type no change in the 
method of dimensioning may be required and the problem is solved simply 


by the addition of suitable notes or symbols indicating the intention of the 


designer as to functional datum lines. 
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It is sufficient to establish datum lines in the case of parts which are 
practically flat pieces with little depth but when a part has substantial 
depth it will be noted that center lines or other datum lines on a drawing 
really represent planes in space. In such parts it becomes necessary to 
establish datum planes rather than lines and three planes at right angles to 
each other are required. 

Figure 8 illustrates such a part which might be an armature such as is 
used in many pieces of electrical contacting apparatus. In the typical 
operation of such a part its functioning is determined by the relation of its 
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Fig. 8—Functional datum plane dimensioning of magnetic armature type of part 
Functional datum plane I passes through the common axis of the two .1285 in. diameter 
holes and .265 in. above the pole face gauge position. 
Functional datum plane II is perpendicular to plane I and passes through the common 
axis of the two .1285 in. diameter holes. 
Functional datum plane III is perpendicular to planes I and IJ and passes midway 
between the finished surfaces which are 1.578 + .003 in. apart 


various dimensions to the position of the pole face and the axis support. 
In order to indicate this on the drawing it is necessary to establish dimen- 
sioning as shown and add to the drawing the notes shown. 

These notes establish three functional datum planes, the first through 
the axis at the point of support and a distance .265” from the pole face 
area; the second at right angles to the first through the axis of support and 
the third at right angles to both the first and second and halfway between 
the finished surfaces 1.578” apart. With these planes established the 
application of all the limits and tolerances shown is based on the operating 
position and analysis of the design is simplified. The drawing and the 
intentions of the design engineer cannot be misunderstood. 
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The clear expression of the designer’s intentions by datum plane dimep. 
sioning will be appreciated by all concerned with the drawing or the r 
sulting part. Inspection of the part is expedited no less than productio 
The inspector can usually by means of gauge blocks or simple fixtures se: 
the part up on a surface plate as indicated by the drawings datum planes 
and positions. He can then establish the conformance of the part with th 
drawing by simple measurements to the indicated horizontal and vertica 
planes. When production quantities justify special gauges the required 
design of the gauge is established clearly by the datum planes. 


Invariable or Gauge Dimensions 


The drawing of Fig. 8 just described illustrates the use of gauge dimen- 
sions. The dimensions .265” and .718” and the indicated half-inch d 
ameter for the pole face are all gauge dimensions without tolerances and 
some statement must be made or understanding reached that they ar 
considered invariable and tolerances not permitted. They represent, 
might be said, theoretical dimensions, on the drawing, or in practice they 
represent tools or gauging apparatus made to the highest standards of 
accuracy. ‘These invariable dimensions are necessary in order to establis! 
a starting point for the dimensioning of the part. It may appear at first 
that stating that a dimension has no manufacturing tolerance or variatior 
is a hardship upon the manufacturer but this is not really so because th 
dimensions are not ones which are actually manufactured in the part. They 
represent usually dimensions built into tools or gauging equipment whic! 
are made toa precision greatly superior to that represented by part tol 
erances. 

Invariable dimensions, or better, gauging dimensions or whatever it is 
proposed to call them are really not a new invention and it is possible to 
cite easily recognized examples. For instance, the dimension 2.473” on 
Fig. 8 isan invariable gauging dimension not associated with the setting 
up of datum planes but typical of long standing use of invariable dimensions 
We all can recall also the use of the term “theoretically correct position” 
and it is present practice in the case of vacuum tube bases and similar ap- 
paratus to designate the location of the contact studs in terms of a gauge 
having holes located on “‘true centers.’”’ Last but not least a minimum or 
maximum limit in its application is itself an invariable dimension. 

In effect, datum lines or planes established when necessary by use ol! 
invariable or gauging dimensions remove the uncertainty as to the de- 
signer’s intentions and prevent misunderstandings between design, pro 
duction and inspecting engineers. Admittedly they do not completely 
solve all problems of dimensions as probably nothing will. They do, how- 
ever, transfer whatever problems remain from the field of tolerances 01 
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finished product to the realm of tool making tolerances and gauging toler- 
ances. The problem of how invariable is “invariable’’ remains but we are 
obviously then considering differences of an order of magnitude not usually 
vitally significant in the functioning of product parts. Theoretically, all 
“invariable” dimensions should be taken to the best accuracy of good 
gauging methods which means that any differences of opinion will be re- 
duced at least to one-fifth and probably to one-tenth of the order of magni- 
tude of those where tolerances themselves are involved. 

It will be necessary to specially identify gauging dimensions on drawings 
to distinguish them from ordinary unlimited dimensions and to indicate 
that they are dimensions for gauges to which only gauge tolerances apply. 


Practical Use of Datum Lines and Planes 


It is not usual to establish datum lines on all d‘awings but if their use is 
necessary in the layout and design of the part they need to be permanently 
identified. This use of datum lines and planes on drawings, where neces- 
sary, may require somewhat greater drafting effort in the actual production 
of the drawing but their use results in a simplification of design and of the 
work of those subsequently using the drawings. It reduces the effort ex- 
pended in analysis of drawings preparatory to the construction of tools 
and minimizes the possibility of misunderstandings or errors in tools. In 
products manufactured only intermittently it is particularly valuable as 
it minimizes the need for understandings and instructions supplementary 
to the drawings which may be forgotten between production periods or lost 
through shifts in personnel. 

The overall economy in engineering effort and the reduction of the numer- 
ous possibilities of error more than compensate for the increase in the actual 
work of indicating datum positions, lines or planes upon drawings. In 
addition the choice of design of punches and dies and similar tools by pro- 
duction engineers is better guided by the designer’s requirements if func- 
tional datum lines are clearly identified. An obvious example is the use 
of either the inside or outside of a punched and formed part as the starting 
point. In brief datum plane dimensioning is a more explicit expression on 
the drawing, of the designers ‘end point requirements”. 

When establishing datum planes, it is important to consider them in 
terms of the actual physical part rather than in terms of the drawing. 
Lines which appear as definite points on a drawing may not be actually part 
of the product when it is completed or may be on surfaces shown as a line 
on the drawing but rough or unfinished in the part. It is difficult to es- 
tablish any set of rules covering what shall or shall not be done because 
each drawing and each part must be considered practically as an individual 
case. That this is so will be amply demonstrated by a serious study of even 
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one part. However, there are obvious generalities which can be established 
and Fig. 9 shows some of them. 

An example of functional datum plane analysis and dimensioning in three 
dimensions of a complicated part is shown by Fig. 10. This is the die cas; 
frame for a special selector switch. It is the base upon which many in. 
terrelated parts and subassemblies are mounted. The proper functioning 
of the completely assembled switch depends in large measure on proper 
manufacture of this casting. In effect, the switch is designed around a 
vertical shaft passing through points P and Q and planes 1 and 2 ar 
therefore, established through the axis of this shaft. The production 
planning engineers intend to design the die and withdraw die plugs from 
such directions that the mounting surfaces will be smooth, flat and without 
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Fig. 9—Types of datum positions 


any taper and they intend to use these surfaces as guiding points for their 
jigs and fixtures. It is for this reason that Plane 1 is established paralle! 
to these mounting surfaces and an indicated distance from them. The 
other planes are established as shown on this drawing and described by the 
notes. With this arrangement of planes the designer’s analysis in terms of 
Plane 1 is easily worked out and the reference of Plane 1 to the mounting 
surfaces permits the preduction or tool engineer to translate the design of 
the part into the design of his tools without necessity for further analysis 
and without the possibility of different interpretations. It will be noted 
that invariable or gauge dimensions are again used. The complete draw- 
ing of this part is very complicated and occupies a drawing practically 4 
ft. x 6 ft. The perspective sketch shown and the accompanying notes are 
incorporated in the drawing as a separate view. 
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Fig. 10—Functional datum planes of complicated switch frame 

Dimensions to datum line “X” or “Y” of the drawing of the frame refer to functional 
datum planes 1, 2 or 3 described below. Points “P” and “Q” are gauge points used in 
establishing these datum lines and planes. Points “‘P” and ““Q” shall be half-way between 
the surfaces “A” and “B” and “C” and “D” respectively and 4.358 in. from the plane of 
surfaces ‘‘M” and ‘“‘N”’ on the mounting lugs. 

Datum line ““X” shall pass through the points “P” and “Q”. 

Plane 1 shall be parallel to surfaces ‘““M” and ‘‘N” and shall include datum line “X” 

Plane 2 shall be perpendicular to plane 1 and shall also include datum line “X” 

Plane 3 shall be perpendicular to plane 1 and to plane 2 at the point “P” 
Datum line “Y” passes through point “P” and is the intersection of planes 2 and 3 
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REQUIRED STANDARDIZATION 


It is not suggested that the drawings shown and the notes referred | 
represent a final practice on datum planes. A standard practice in desiy 
nation of planes and standard terminology and understanding on gaug 
points and gauge dimensions is required. It will probably be desirable | 
adopt some symbol or designation for use on drawings to distinguish gaug: 
dimensions which are invariable from ordinary unlimited dimensions 1 
which manufacturing engineers for their own purposes usually add sho; 
tolerances. One thing is certain and that is that datum planes, dimensions 
and tolerances when established should be primarily in terms of the 
quired functioning of the apparatus. When that is done no one using th: 
drawing in any capacity will have any doubts as to the designer’s intentio: 
and this results in a great reduction in the discussions and analysis whic! 
might otherwise be necessary. 


SUMMARY 


In summary it may be said that the whole approach to these problems i: 
dimensions and tolerances should be on the basis of functioning. However, 
good engineering of dimensions and tolerances requires knowledge of what 
can reasonably be produced and the sources of reasonable tolerance values 
are: 

1. Raw material limits including some knowledge of future trends anc 

developments. 

2. The normal accuracy of manufacture, also including anticipation o/ 

future improvement. 

3. Discussion of trend of design with manufacturing engineers. 

Solution of tolerance problems in the final design may involve all of th: 
following steps: 

1. Study of the effect of combinations of tolerances on functioning, 

allowing for statistical effects in accumulations of tolerances. 


me 


Discussion of this analysis with the production planning engineer 
because the analysis of tolerance combinations is important in the 
design of long life tools. 

3. Indication of the results of such an analysis by the method of dimen- 

sioning drawings. 

4. Indication on drawings of functional datum positions, lines or planes 
established on geometrically correct principles to permanently and 
unmistakably record the intentions of the designer regarding com- 
binations of variations wherever this is necessary. 
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Time Division Multiplex Systems 
By W. R. BENNETT 


INTRODUCTION 


HE idea of transmitting and receiving independent signals over a 
common line by means of synchronized switches at the terminals is 
quite old and has been used in multiplex telegraphy for many years. In 
general if V signal channels are to be provided over one line, the switching 
cycle includes N equal time intervals, one of which is allotted to each 
channel. Each channel is connected to the line throughout a part of its 
particular time interval and is disconnected throughout the remainder of 
the cycle. Absence of interference between the channels depends upon 
the fact that the channels are connected to the line throughout mutually 
exclusive time intervals. It is thus possible to avoid the use of channel 
band filters such as are necessary in carrier systems employing frequency 
as the basis of separation. 
Application of time division multiplex methods to telephone channels 
has been proposed from time to time and some experiments have been 


1,2,3,4,5,6 
made. 


It is fairly evident that the concept of simple on-and-off 
switching giving alternately transmission and complete suppression for the 
signal from a particular channel on the line is inadequate for speech waves 
in actual telephone circuits. Imperfections in the transmission properties 
of the line tend to distort the wave form of the successive signal components 
and prolong the contribution of one signal into the time allotted for a differ- 
ent channel. It is the object of this paper to present a general quantitative 
discussion of the factors which enter into the transmission of any type of 
signal by a system of this kind. It has been found possible to arrive at 
definite criteria for such matters as the required switching frequency, the 
conditions to be imposed on contact time for good crosstalk suppression 
with economy of frequency band, and the transmission requirements which 
must be met by the intervening circuit to hold the interference between 
channels to tolerable values. The analysis leads directly to a physical 
viewpoint of the whole process which, to those familiar with the carrier and 


' Patten and Minor, U’. S. Patent 745,734, 1903 

* Electrical World, Dec. 5, 1903. 

3 Goldschmidt, U. S. Patent 1,087,113, Feb. 17, 1914. 

* Poirson, Soc. Fr. El., Apr. 1920. 

5 Marro, L’Onde Electrique, Oct. 1938. 

°M. Cornilleau, Revue de Telephones, Telegraphes et T. S. F., 13 (1935), pp. 625-643 
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sideband philosophy of signal transmission, illuminates the manner 
which departures from ideal amplitude and phase characteristics cause 
crosstalk between the several message channels. It further leads direct]; 
to other physical methods for producing and detecting a transmitted signa! 
identical with the essential components derived in time division or switching 
processes. 

A first step in the theoretical solution of the problem was taken by 
Dr. J. R. Carson, who, in an unpublished memorandum of May 25, 1920, 
derived quantitative relations between band width and interchannel inter 
ference in time division multiplex transmission. Applying Fourier series 
analysis to on-and-off switching, he showed that if the transmission mediun 
had constant attenuation and linear phase shift for all frequencies below 
cutoff and no transmission of frequencies above the cutoff, the band width 
required for satisfactory multichannel telephony would be much wider than 
needed in conventional carrier methods. A further step was taken by 
Dr. H. Nyquist, who, in unpublished memoranda of August 24, 1936 and 
November 12, 1936,’ showed that the width of band necessary may be 
reduced by providing a specially devised type of non-uniform transmissio1 
characteristic. In the following discussion, we shall see that a similar 
result can be obtained by control of the switching, and specific switching 
processes will be described which allow a flat transmission band of minimum 
width to be used. 

In order to arrive at requirements which must be imposed on the various 
components of the system, we shall first give a theory of time division 
multiplex transmission in which both the switching processes and_ the 
transmission characteristic are completely general. Specific forms which 
give crosstalk suppression will then be discussed and effects of small de 
partures estimated. 


GENERAL THEORY 


We shall assume an \-channel system with a sinusoidal signal impressed 
on the j® channel. An illustrative arrangement is shown in Fig. 1. Since 
the system is linear, we may represent currents and voltages by complex 
quantities with the understanding that the actual currents and voltages 
are the real components of the expressions used. Accordingly, let the 
signal voltage impressed on the j** channel be 


E,(t) = E;e"*'' (1 


7 Basic concepts used in Nyquist’s analysis were included in his paper, ‘“Certai: 
Topics in Telegraph Transmission Theory,” A. J. E. E., Trans., April, 1928, pp. 617-644 
Mention is also there made of the equivalence of signal shaping and equalizing in effect 
on reception of telegraph signals. 
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TIME DIVISION MULTIPLEX SYSTEMS 


and let the switching between the j"® channel and the line at the sending 
end be represented by: 


I 


where J,;(t) is the current flowing into the line from the j"® channel. The 


sit) = F,(HE;,(0), (2) 
function F(t) has the dimensions of an admittance and, in the arrangement 
shown in Fig. 1, is periodic in time with fundamental frequency g = 2%/T 
radians per second, where T is the time eccupied by one cycle of the switch- 
ing operation. In the interests of economy of analysis, it is preferable for 
our purposes to assume for F';(/) a somewhat more general function of time 
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Fig. 1—Elementary arrangement for time division multip!ex system 


than is directly obtainable with the elementary arrangement of Fig. 1. 
We shall let 


oo 


F(t) = ps Am; cos [(v + mq)t — On;). (3) 


m=0 
To make the results applicable to Fig. 1, we merely let v = 0; then by the 
usual Fourier series analysis, 


Aoj = 40/2, Anji = Onj + Oni 
7m > 0 
60; = 0, tan O0nj = On;/anj |} 
9 T+t, } F 
amj = =| F(t) cos mgt dt i \4) 
1 


T+t, 
: F(t) sin mgt dt, t; arbitrary 


‘1 











a 


a 


Bia PP MB ACI 





202 BELL SYSTEM TECHNICAL JOURNAL 


The wave (3) consists of the output of the circuit of Fig. 1 with all freque; 
cies shifted by a constant amount v radians per second; various means 0; 
accomplishing this result in the switching process will be discussed |ater 
It is sufficient to point out here that such a shift in frequency is often ¢ 


sirable for optimum utilization of the transmission medium. Combining 
and (3), we then have: 
E, 2 ! 
4} i(v+mq+wj;)t—i0 i(v+mq—w;)t+16 
I ,;(t) oa ; ; Amjle mqtw; mij + Pe q-@; me 


= m=0 


a 


It is clear from (5) that the result of the switching process is the produ 
tion of upper and lower side frequencies from the signal on each harmoni 
of the switching frequency. It is also evident that if more than one signa 
component is superimposed, the resulting side frequencies constitute sic 
bands of the same nature as used in amplitude modulation systems. 
significant difference between time division and amplitude modulatio: 
appears in that in the latter only one sideband or at most one pair of side 
bands is transmitted, while the essential character of time division depends 
on the transmission of a plurality of sidebands. Thus if one pair of side 
bands were selected from the output (5) by filtering, the time divisior 
process would merely be a particular way of generating the sidebands 
required in an amplitude modulation system. 

The next step in a time division system is the transmission of the wave (5 
over a line. The properties of the line in general may be specified by a 
complex transfer impedance, which we may express here by the ratio oi 
open-circuit output voltage to input current: 


E,/I, = Z(w) 6 


The result of applying the wave (5) to the line is then the open-circuit 
voltage: 


. E “ ets i(wtmg tw) t—6 
E,(t) = , pe AmjZli(v + mg + werner mt 
& m= 


> 
+ 7 p> AnjZ*[i(v + mg — u;)Je 0m 

In the above we have adopted the notation Z*(iw) to represent the con 
jugate of Z(iw) and have made use of the fact that the response of a network 
to the applied wave e “‘ is the conjugate of the response to e 

At the receiving end another switching process takes place synchronous!) 
with that at the transmitting end. We shall assume that the switching 
process between the k* channel and the line is represented by the relatior 
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In(t) = G.(HE,; (0), te) 
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where J,,(t) is the current received in the &* channel and G,(¢) is a periodic 
function of time with fundamental frequency g. It is understood that j 
and k may be any two of the N channels. We shall express G;(¢) in a 
manner analogous to the corresponding switching function at the trans- 


mitter, i.e., 


G,(t) = > Bix cos [(v + ng)t — ®,,) (9) 


n=( 


Combining (7), (8), and (9), we find 


T x(t) = 7 bs ps Amp Bux Z\ilv + mq os w;)] 
n=0 


m=() 


ai ac an 
if2v+(mtn)qtwjlt—j(Omj+Pnk il(m—n)q+wj]t—i(m j—-Pnk 
(e +e ) 


+ _ D dD Amp BueZ*[i(v + mg — w;)] 


m=0 n=0 


(e jl2et (m+n) q—wj)t+i(Omjt+Onk as elm m)qtwj]t+ilOm j—Pnk ) (10) 


The received wave thus consists of a doubly infinite set of side frequencies 
involving harmonics of g. It is, however, possible to set up conditions 
under which the original signal may be selected and the frequencies involving 
the switching rate may be suppressed by filtering. If v = 0, such separation 
is possible provided 


w; < 9/2, (11) 


for it then follows that a low-pass filter with cutoff frequency at g/2 will 
not pass any of the components with frequencies dependent on g. The 
condition (11) follows from the fact that the lowest frequency of (10) 
dependent on q is g — w;, and hence we must make g — w; > w; in order to 
separate w; from g — w;. In other words the sidebands on adjacent har- 
monics must not overlap. If v > 0, the condition (11) also suffices as far 
as suppression of terms dependent on g are concerned, but an additional 
condition is required to suppress frequencies dependent on » in the special 
case in which v < q/2, i.e., the case of v less than the maximum allowable 
value of w;. For in the latter case the frequency 2v + (m + n)q — w; is 
less than g — w; in the special case of m = n = 0. The additional condition 
needed is evidently either Aj; = 0 or By = 0. If v = Oor if vy > g/2, 


this condition is unnecessary. 

Assuming then that (11) is fulfilled, and that a low-pass filter with cutoff 
at g/2 is inserted in the output of each channel, we calculate for the typical 
channel output current: 


Tea(t) = V jE je", (12) 
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where the value of Y ;, is as follows: 
Case 1, vy = 0 


V 54 me 4 { 0) Box Z| 14;) 
oo 


a : Anj Bmx Z[i(mg + w,) le Omi Pm 
1 


m 


i< wine ae 
+ ; pe Aaj Bm Z*[i(mg = w,)le! Omi -Pne 13 
Case 2, v > 0; Ao; or Box = 0 if v < g/2 


Y, = ; a Anj Bn Ziv + mg + w;) le 1m j— @mb) 


1 = a ° (A >,,. 
rs : > Ani BmiZ*[i(v + mg — w)JeOmi*™ (44 
m 0 

The combination of an \-channel time division multiplex system wit! 
low-pass filters in the receiving branches is thus found to be equivalent t 
a linear network having N pairs of input and output terminals with th 
transfer admittance from the j‘" pair of input terminals to the &* pair oi 
output terminals given by J’, in (13) or (14). The transfer admittance is 
calculated by summing the contributions of upper and lower sidebands o1 
harmonics of the switching frequency and is affected directly by the trans- 
mitting properties of the medium at the side band frequencies. The result 
we have obtained is of sufficient generality to include all cases we shall 
treat in this paper. We shall now proceed to specific examples. 

ON-AND-OFF SWITCHING WITH COMMUTATOR 

When an ideal commutator is used as a switching means, the switching 
functions for the NV channels are identical except for a time displacement 
which is the same between all pairs of consecutive channels. This condition 
is expressed by: 

F(t) = Ft — (j — 1)T/N] (15 
Thus F,(¢), the switching function for the first channel becomes a reference 
function, /2(¢) is the same except for a time delay of T/N, F(t) is delayed 
by 27/N, etc. Substitution of (15) in (3) gives the relations: 
A mj = Ami ) 
Omj = Omi + (FG — 1)2mx/w | 


If we further suppose that the commutator makes contact between the 
typical channel and the line throughout a fraction x of the time interval 
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TN allotted to that channel and breaks contact throughout the remainder 
of the switching cycle, we may write the reference switching function as: 
A, —aT/€n <t < xT/2N 
F(t) = eae ; ” " (17) 
0, «7T/2N <t < (2QN — x)T/2N 


+ Hence from (4) 


la = Ax/N, 
2A . mxn 

Am = sin —— , m > 0 > (18) 
mar N 

Oni = 0 


In the receiving device, the corresponding switching process should be 
delayed with respect to the transmitter by a time interval f) equal to the 
time of transmission of the line. Hence we write 


Bu = Bx/N 


2Be . mxr 
Bas om —- 
mr 


tre? (19) 


Pri = mato 


Bn; and 6, ; are related to B,. and ,,; in a manner analogous to (16). 

The time of transmission of a distorting line is not precisely definable, 
but may be represented for our purpose by a linear phase component of 
Z(iw). That is, we write 


Z(iw) = Zoliw)e ““*, (20) 


where f) is the slope of a straight line giving the best linear approximation 
to the phase vs. frequency curve, and Zo(iw) is the impedance function 
remaining after the subtraction of fow from the actual phase shift ordinates. 
Substituting (15)—(19) in (12), we find 


Vix = ABe row; (= Z(1w;) 


=. sin’ man /N 


r . 2m N 
oa be —_— Zoli(mg + w;)le ii 
m=1 m* 1° 
sin’ mxx/N _*,. i( j—k)2ma/N 
——. Zo[i(mg — w;)]e"” , (21 
m=] m* 1° 





\ 
} 


















remneecs 


Cr RE TREN AE Fs 











206 BELL SYSTEM TECHNICAL JOURNAL 
If the attenuation of the line is constant throughout the range w; to Mg ~ 
w; and all frequencies above the latter value are suppressed, (21) becomes 


’ 2 e itqw; M Fe IN 2 “= 
Yu = ABs =e E +2 pm (= ett ) cos (7 — k)2mr | 


m=1 mx /N 


The crosstalk ratio or ratio of amplitude of signal received in the k* chan 


to that received in the 7" channel when signal is transmitted in the 
channel is, therefore, 


M _s fe 2 
: 1+ 2 b he (= _ if) cos 2mr(k — 7)/N 
Vix m= \ ma x/N : 


Vij “ (sinmarx/NY ; 
1+2>( 


m=1 mr x/ N 


Results of calculations made for a 10-channel system from (23 
x = land x = .5, corresponding to no lost time and half lost time respe 
tively in switching are shown in Fig. 2. It may be noted that adjacen 
channel crosstalk with half lost time is equivalent to alternate chann 
crosstalk with no lost time. Examination of the curves reveals a number 
of significant facts, among which are: 

1. Crosstalk is quite imperfectly suppressed when the band width of +! 
line is smaller than the theoretical minimum—the width of one sideban 
multiplied by the number of channels. 


2. As the band width of the line is increased above the theoretical mi: 
mum, improvement in crosstalk suppression increases slowly, so that ii 
general the use of frequency range on the line is uneconomical compare: 
with other systems. For example, with no lost time in switching, the ban: 
width of the line must be increased tenfold to suppress adjacent channe 
crosstalk by 40 db. This conclusion is, however, to be qualified as follows 

3. When the duration of contact is decreased (less of the available channe 
time used) definite optimum transmission band widths appear which giv: 
better crosstalk suppression than bands somewhat wider or narrower. ‘This 
suggests the possibility of critical phase relations existing between the con 
tributions from the various sidebands such that if the right number having 
proper amplitudes and phases can be combined, complete suppression o! 
crosstalk may occur even when the transmitted band width is finite. 

When x, the fraction of contact time used, is made to approach zero, 
the limit of the amplitude factor (18) for the typical harmonic of the switch 
ing function is Am; = 2Ax/N, which is independent of m. This is con 
sistent with the known fact that a wave consisting of periodically repeated 
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amplitude. If we use very short contact durations in time division, we 
should accordingly expect a large number of sidebands of nearly equal 
amplitude. The combination of proper numbers and phases of these side- 
bands offers a key to the realization of a time division multiplex system 
giving good crosstalk suppression with economy of frequency band. 
Suppose that the duration of contact time is made sufficiently small to 
| realize approximately the limiting values Am, = 24x/N, Bny = 2Bx/N in 
transmitting and receiving respectively for the first 2M + 1 of the side- 
bands and that by means of a low-pass filter with linear phase shift and 
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NUMBER OF SIDEBANOS TRANSMITTED 


Fig. 2—Crosstalk between adjacent channels of ten-channel time division multiplex sys- 
tem with on-and-off switching. No attenuation or phase distortion in pass band of line 


uniform attenuation in its pass-band all other sidebands are removed from 
the line. The expression (23) then becomes: 


M 
1+2 > cos 2mmn(k — 7)/N 


Vix m=1 
Y;; 1+2M 
| 1, k= jj 


}sin (2M + 1)4(k — j)/N 
\(2M + 1) sin x(k — j)/N’ 


(24) 
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In particular, if 
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M = (N — 1)/2, ) 
it follows that f Ys 
V x/V;; = O,k ¥ j) 


Thus there exists in theory a system employing sidebands on zero fre. 
quency and the first (V — 1)/2 harmonics of the switching frequency, j; 
which multichannel transmission is possible without interchannel inter 
ference. Since the required condition (25) may also be written \ 
2M + 1, an odd number of channels is obtained. Since N sidebands ar 
transmitted, the band width used is the same as the minimum required for \ 
single sideband amplitude modulation channels on a frequency discrimina 
tion basis. Sidebands produced on higher harmonics in the time divisi 
process must be removed by filtering. 

It is to be noted that since it is equality of the V sideband contributions 
which is important and the amount of each contribution is determined }y 
the transmission characteristic of the line as well as the transmitting and 
receiving switching processes, it would be theoretically possible to make wy 
for sideband irregularities by equalizing the line. However, the equaliza- 
tion required in the line would be of “stairstep” type rather than smooth 
varying with frequency since an error in the value of one harmonic of the 
switching function produces the same error throughout the entire range 
occupied by the pair of sidebands associated with that harmenic. 


GENERAL SWITCHING FUNCTIONS WITH CROSSTALK SUPPRESSION AND 
Mixtmum BAND WIDTH 
The above discussion based on the properties of a commutator has led 
us to an ideal switching function which is, except for an unimportant 
proportionality factor, 


(N—1)/2 
F(t)=1+2 )) cosml[qt — (j — 1)2x/N], N odd (26 
m=1 

This type of switching is approximately realizable with synchronized com 
mutators having contact widths very narrow in comparison with the spacing 
between contacts. For a 3000-cycle speech band, the minimum switching 
rate would be 6000 cycles per second. Such a speed would be difficult to 
obtain with ordinary mechanical means but would be feasible with rotating 
electron beams. 

The concept of combining detected contributions from a number ol 
sidebands in proper phase to give in-phase addition of desired components 
and cancellation of unwanted ones leads to a generalization of the switching 
processes over those possible with synchronized commutators. We note 
that the switching functions F(t) of (2) and G,(t) of (9) are analogous to 
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carrier waves applied to a product modulator, and an electrical analogue 
of time division may be realized therefore by applying signal and a suitable 


carrier to a product modulator. Phase shifts in the carrier supply circuit 
mav be made to serve the same purpose as the angular displacements 
between commutator segments. It is thus of interest to examine various 
other possible forms of the function F(t) which are suitable for multiplex 
transmission and investigate methods by which they can be realized. 

We note that (26) is suitable for an odd number of channels because it 
makes use of the direct signal component (or sideband on zero frequency) in 
addition to the paired sidebands on harmonics of the switching frequency. 
It seems reasonable to expect that systems for even numbers of channels 
can be devised using only upper and lower sidebands on harmonics and 
omitting the signal itself. Complete information for the separation of .V 
channels should be contained in any set of V sidebands; hence we should 
not be forced to start with the sidebands of lowest frequency, but be able 
to use other sets with a more suitable place in the spectrum or with better 
equalization of amplitudes. 

We shall derive an expression for a quite general switching function 
meeting the desired conditions of freedom from crosstalk and economy of 
band width for an even number of channels by assuming the following 
forms for A»; and @,,; in (3), 

A,nimsn+N/2-1 

Anj = (27) 
O,m <n,orm>n+ N/2—1 

Omi = (FJ -1) (m+hyrt+a (28) 


The switching function assumed contains V/2 harmonics and hence will 
produce N sidebands. The values of n, /, a and y are first assumed to be 
arbitrary. At the receiving end, a switching function similar except for a 
time displacement fo will be assumed. That is, in (9), we take 
Bhnsmsn+N/2-1 


Bink _ (29) 
0,m<norm>n-+N/2—-—1 


Pre = (k — 1) (m+ hy + (mg+r)r + (30) 


Transmissicn over the line is assumed to be of the distortionless form ob- 
tained by setting Zo(iw) = Zo, a constant, in (20). Substituting (27)-(30) 
in (14), we then calculate 


j= 


Y it 


1, 
_NABZye J2 sin N(j — k)y/4 cos (4n + 4h + N — 2)(j7 —k)y/4 
4 | N sin (7 — k)y/2 


j#k (31) 
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Freedom from interchannel interference is obtained if the numerat,; 


vanishes and the denominator does not vanish for all unequal integer valyes 
of j and k less than V. We note that this can be accomplished by setting 


Ny = (4n + 44+ N — 2)p = 2r, 3) 


in which case the numerator becomes 2 sin (7 — k)x/2 cos (j — k)x 
sin (7 — k)r = 0. Solving (32) for h and y, we find 


h = —(2n — 1)/2| 
> (33 

y = 2n/N } 
The denominator thus becomes N sin (7 — k)w/N, and since the larges 
possible value of 7 — k is N — 1, the denominator cannot vanish for j = | 


The value of m remains arbitrary; hence we may start with any harmoni 

we please. The value of a is also immaterial. The general form of switc! 

ing function thus derived is, omitting the constant of proportionality: 
n+N/2—1 


F(t) = >> cos [(v + mg)t — (j — 1)(2m — 2n + 1)e/N + al, 


m=n 
n = (0,1, 2,---, N even 


Dependence of the phase angle upon the initial harmonic may |i 
avoided in certain special cases; for example, if we set n = rN/2, where» 


may be zero or any positive integer, we obtain the result that a particular 


switching function which satisfies the required conditions is: 


(r+1)N/2—1 


F(t)= >) cos [(v+ mg)t — Gj — 1)(2m + 1)x/N + al, 


m=rN /2 


w 
‘at 


ry = 0,1,2,---, N even 


Specific methods of realizing satisfactory switching functions will now be 


examined. 


PLUS-AND-MINUS SWITCHING WITH COMMUTATOR 


Instead of opening and closing the circuit between each individual 
channel and the line, the commutator may be designed to reverse the 
polarity of the connections on alternate contacts with a given channel. An 


appropriate switching function for alternate polarity reversal is: 
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0, -—(2N—x)T/2N <t<—xT/2N 


F(t) = 7 2 . AGES 
0,2x7T/2N <t< (2N —x)T/2N 


— A, (2N —x)T/2N <t<(2N+2)T/2N 





A, —2T/2N <t<xT/2N | 
| (36) 
| 


[he relation of F(t) to F(t) is as specified by (15). By the usual Fourier 
series expansion, we then find 


co 
, . (Qm + 1)xxr 
F(t) sin : cos (2m + 1)gt/2 (37) 
Mf) = 25 2m * 4)s 2N q 
where g = 2n/T. In the limiting case in which the duration of contact 


time is made small compared with the interval between contacts, the 
expression (37) combined with (15) gives the general result: 


2Ax < (2m + 1)gt (2m+1)(j —1)4 
teat Le (38) 
4 m=() 2 N 
By comparing this function with (35), setting v = g/2, we note that plus- 


and-minus switching with small contact time can be made to give freedom 
from crosstalk and economy of band width in a time division multiplex 
system provided that we select the sidebands on V/2 adjacent odd harmonics 
of q/2, beginning with m = rN/2, i.e., the (rN + 1)* harmonic of q/2, where r 
may be sero or any positive integer. It is necessary that even harmonics of 
q/2 be suppressed, which will be the case if the positive and negative 
switching contacts are identical except for change in sign. 


VESTIGIAL SIDEBAND SYSTEMS 


It is seen that the use of a commutator results in the production of a 
large number of sidebands, and that for economy of bandwidth in the trans- 
mission medium it is desirable to select only a definite group by filtering. 
The requirements on the filter are not easily met since there should be linear 
phase shift and uniform attenuation throughout the range of frequencies 
occupied by the desired sidebands and a high degree of suppression of all 
other frequencies. The requirements may be simplified somewhat by using 
the principle of vestigial sideband transmission which allows the attenuation 
to increase gradually within the desired band and compensates for energy 
removed from wanted sidebands by allowing a corresponding amount of 
energy to be transmitted from adjacent sidebands which would normally be 
suppressed. Linear phase shift is required for all frequencies at which 
appreciable transmission exists. 
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Figure 3 shows the curve of admittance vs. frequency required for on-a 

off and plus-and-minus switching. Referred to the mid-band admittance 4: 
unity, the admittance is reduced to one-half (six db loss) at the frequencie: 
rNq/2 and (r + 2)Nq/2 which are the nominal upper and lower cutoffs 
N is an even integer and r is zero or any positive integer. The admitta: 
curve has odd-symmetry about the cutoff frequencies—that is, if 
frequency « cycles below a cutoff frequency, the admittance has the valu 
it must be 1 — a at a frequency x cycles above the cutoff. The noming 
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Fig. 3—Vestigial sideband transmission in time division multiplex systems 


band transmitted in the case of on-and-off switching consists of the upper 
sideband on the harmonic rNq/2, the lower sideband on the harmoni: 
(ry + 2)Nq/2, and all intervening sidebands. In the case of plus-and-minus 
switching the upper and lower sidebands on frequencies (rN + 1)q/2 to 
[(r + 2)N — 1]q/2 inclusive are transmitted. Impairment of the nomina 
band by the filter is made up by transmitting the appropriate parts oi 
sidebands outside the nominal range. It is easily verified that either of the 
systems depicted in Fig. 3 satisfies the required conditions for multiplex 
transmission without interchannel interference when the sidebands produced 
by a given signal have equal amplitudes over the range utilized. 
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[ne vestigial method is required only when strong sidebands very near 
the desired ones must be removed. Modifications of the time division 
process exist in which vestigial filters are unnecessary because very little 
energy appears at unwanted side frequencies. It is clear that if we regard 
the problem as one of producing certain sidebands on carriers of definite 
phases, we are not restricted to commutating devices only but may make 
use of general modulator technique. Further details concerning specific 
circuit arrangements are described in LU’. S. Patent 2,213,938, W. R. Bennett; 
and U. S. Patent 2,213,941, E. Peterson. As a general guide the following 
table of carrier phases for an .V-channel system (.V even) is furnished: 


TABLE OF PHASE SHIFTS FOR N-CHANNEL SYSTEM 


(NV /2 Carrier Frequencies Required) 


Carrier Frequency v v+q v + 2g v + 
Phase {Channel 1 0 0 0 0 
Shift 2 r/N 3a/N 5a/N in/N 
In 3 2a/N 6n/N 107/N l4n NN 
Carrier | 4 3ar/N 9On/N 152/N 21x/N 


TRANSMISSION REQUIREMENTS 


Practical success of a time division multiplex system requires the main- 
tenance of a satisfactory ratio of wanted signal to crosstalk. In order to 
accomplish this, the transmission link must preserve the amplitude and 
phase relations of a group of sidebands. A physical picture of the relations 
involved may be obtained from Fig. 4, which is drawn for the particular 
case of a 5-channel system of the on-and-off switching type. For this 
example the theory previously developed shows that five sidebands of 
equal amplitude are sufficient, namely—the signal itself (which may be 
regarded as a sideband on a carrier of zero frequency), the upper and lower 
sidebands on the switching frequency and on the second harmonic of the 
switching frequency. If we take the phases of the switching fundamental 
and its second harmonic as applied to the first channel as a reference, the 
proper phases of fundamental and second harmonic respectively for the 
other four channels are given by the following table: 


Channel Fundamental Second Harmonic 
Number Phase Phase 

2 72 144 

3 144 288 

4 216 432 

5 288° 576 


In Fig. 4, we have assumed a single-frequency signal component as input 
to the first channel. If the line has distortionless attenuation and phase 
characteristics, the five resulting side frequencies are received in the first 
channel as the five in-phase vectors of equal amplitude shown in (a). Re- 
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ception in the second channel is shown by (0) in which the vector 1 represents 
the directly transmitted signal component (or sideband on d-c). 2 and 


3 «a 


DETECTED COMPONENTS 





ae a al 





(2) RECEIVED IN CHANNEL 


-288° 2 





[1 ‘432° 


(d) RECEIVED IN CHANNEL 4 


(€) RECEIVED IN CHANNEL 


w 


RESULTANT 

crosstack TENS 
25 DB BELOW 
RESULTANT IN (a) 





(f) EFFECT OF AMPLITUDE ANC 
PHASE ERRORS ON RESULTANT 
OF COMPONENTS RECEIVEC 
IN CHANNEL 2 


Fig. 4—Graphical representation of operation of time division multiplex system. Signa 
transmitted in channel 1 of 5-channel 5-sideband system 


represent the detected components from the upper and lower side frequencies 
associated with the fundamental switching frequency, and 4 and 5 th 
components resulting from the upper and lower side frequencies of the 
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TIME DIVISION MULTIPLEX SYSTEMS 
second harmonic. Components 2 and 3 are shifted —72° and +72° respec- 
tively and 4 and 5 are shifted — 144° and +144° in relation to the phase of 
component 1. As shown in (0), the five vectors combine in the form of a 
closed polygon giving a resultant of zero amplitude. Similar vector dia 
grams for reception in the third, fourth, and fifth channels are shown in 
c), (d), and (e). The appropriate diagrams for transmission in channels 2, 
3, 4, and 5 and receiving in any channel can be obtained from (a) — (¢ 
by cyclic permutation of the channel numbers, i.e., transmission in 1 and 
reception in 2 corresponds to transmission in 2 and reception in 3, etc. 

Production of crosstalk by phase and amplitude distortion in the trans- 
mission medium is illustrated by (f), Fig. (4), which shows the resultant 
component received in channel 2 when signal is transmitted in channel 1 
and an imperfect line is used to connect the transmitting and receiving 
terminals. The vector 1 is taken as the reference amplitude and phase. 
The gain characteristic of the line is assumed to be one dé low at the side 
frequency producing vector 2, one db too high for vector 3, 0.5 db low for 
vector 4, and with no error for vector 5. The phase curve is assumed to 
depart from a straight line by —1°, —1°, —4°, +3° at the side frequencies 
from which components 2, 3, 4, 5 respectively are derived. The vector 
polygon fails to close and the resultant represents an unwanted signal 
received in channel 2 at a level 25 db below the wanted signal received in 
channel 1. 

We may make an estimate of the accuracy of the equalization required 
in the general case by writing the transfer impedance Z(iw) in the form: 


Ziiw) = p(w)Zoe aiid 39 


where Bw) represents the departure of the phase shift from a straight line 
and the variation from flat gain is given by 


g(w) = 20 logio p(w) (40) 


The expression (39) may be rewritten as: 
Z(iw) = [1 + 2(iw)|Zoe a (41) 
where 
am p.. \5 18(w) ) 
z(w) = plwe — 1 (42) 


If we assume that the switching function is of the general form (34), we 
calculate from (14) the general relation: 
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\ 
| a+—-—l 
IN+ >> (sli(v+ mg + w;)] + 2*[i(v + mq —w,)]), j = &| 
Vj, = K \n+5-1 L (43 
. (sli(v + mg + w,)]e seciicsieniaiiinidaeie | 


| m=n 


+ s*[i(v + mg — w;)IJe"” sedeachiadiaaines! i # k| 


where A is a complex constant of proportionality. The case of 7 = k whic! 
gives transmission within the channel contains a variation with signa! 
frequency caused by the summation of the departures from ideal transmis- 
sion at the V sideband frequencies. This term presumably will be unimpor- 
tant if the transmission characteristic is sufficiently good to meet crosstalk 
requirements; hence we may neglect the z and <* terms in the case of j 
and write the ratio of interference to desired signal as: 


n+" 1 
: Vix 1 - : b) (Qmn—tatle/N 
Ca="=— 2, (hip togt+ae erry 
Y;; N m=n 


+ z*[i(v + mq — we POP IN) (44 


The crosstalk ratio will in general vary with the signal frequency. The 
requirement would logically be based on the total interference power 
weighted in accordance with the interfering effect at individual frequencies 
Thus we might set 


®b 
Xjx = S; / W jx(w;) | Cin |” de; (45 
@a 
where X ;, represents the weighted interference power received in the ' 
channel when a reference signal wave of mean power 5S; is applied to the ;' 
channel. The limits of integration w, and w, are the lowest and highest 
signal frequencies used. The function Wj(w,;) represents the proper 
weighting with frequency of the interference and takes into account the 
distribution of the interfering signal and the relative importance of the 
different interfering frequencies. 

Equation (45) is sufficient for computation of interchannel interference 
introduced by the line when the transmission characteristics of the line are 
known. A more valuable result, however, would be the expression of the 
required line characteristics in terms of the allowable interference. In 
general this would require some specification of the nature of the departures 
from the ideal characteristic. Except perhaps for systems with very few 
channels, it seems reasonable to assume that the departures are distributed 
fairly uniformly throughout the frequency range transmitted by the line, 
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and hence that for purposes of estimating requirements we may replace 
C,, in (45) by its average value over the band. We may then write (45) 
in the form: 


Cix = Xin = Uj (46) 


wh 
S; / W jx(w;) dur; 


[he value of the right-hand member, which we have designated by the 
symbol U jx , is either known or can be determined for the particular type ot 
signal. Hence our problem is reduced to finding the allowable departures 
in transmission which keep the mean square absolute value of C, from 
exceeding a prescribed maximum value. 

We note that Cj, is the sum of V complex quantities, each of which is 
restricted to a range of values determined by the transfer impedance of the 
line in an individual band of frequencies. A convenient simplification may 
be made by regarding the V complex quantities as .V independent chance 
variables. This is tantamount to assuming that departures in one band 
do not affect departures in any other band; the assumption is not strictly 
true, but should lead to no important error. We may then make use of the 
following theorem*: If 


$= Bits + date + +++ + bate, ” 


where 21, 22, +++ Z, are m independent complex chance variables and 
b,, bg, --- b, are complex constants, 


9 9 9 9 


¢|; = bil} | aip + eee + ns Sa | 48) 


n+ 
ar 1 : 3 : 2 
Cx | = ny? >> (| c[i(v + mg + w;)] ~ + | 2*[i(v + mg — @;)] |) 
= | 2 P/N, (49) 


if the average square of the absolute value of the departure is the same in all 
bands and is equal to | z |", which we shall define as the average squared 
absolute value of the departure for the entire line band used. 

From (42) and (40), 


| s(iw) |" = 1 — 2p(w) cos Ow) + p (w) 


1 — 2.10° cos 6(w) + 10°!” 50) 


5R.S. Hoyt, B. S. T. J., Vol. XI, No. 1, Jan. 1933, p. 64. 
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Since it seems certain that g and @ must remain small to make the syster 
operative, we investigate the nature of (50) when expanded in powers 
gand@. The leading terms are: 
, 9 
2 (log. 10) «,.. 2 

(iw) | = ——— gw) + Ow) + --- 5] 

400 a 
Hence for g and @ small, we have independent of any correlation which ma 

a ; y 

exist between g and @, 


——, log. 10 ) » 7. 
3( ») . = ‘ a ) A"( 52 
1W) ( 0 g Ww + w) Je 
Let 
o} = gw), 03 = Ow) dS 


Then from (46), (49), (52) 


, 


y log, 10\° ? 
— ( A ) ia a] / N 54 


In (54) o; is the r.m.s. departure of the gain in db from a constant and o. 
is the r.m.s. departure of the phase shift in radians from a straight lin 
If o2 is expressed in degress instead of radians, (54) becomes 


Uj, = 10° (13.25 of + .3046 03)/N 55 


The total interference received in any one channel is the sum of the 
individual contributions from the other NV — 1 channels. The additio: 
factor required to express the total in terms of the interference from on 
channel depends on the nature of the individual loads. Thus if the proba 
bility that any one channel is transmitting a signal wave is 7, the averag 
total interference power received in one channel is 


X = r(N — 1)Xu, = US; [ W ix(w;) dw; , (56 

where 
U = (N — 1)rU jn, = ad — 10 *(13.250; + .304603) (57 
For large values of N, the ratio (V — 1)/N approaches unity, and the 


average interference becomes independent of the number of channels 
The average interference may not be the most significant quantity, however 


For example, if there is a considerable probability that all channels are 


* This method of avoiding any assumption concerning correlation of attenuation and 
phase was suggested by Dr. T. C. Fry. 
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carrving energy simultaneously, as would be the case if the channels were 
subdivisions of a common signal band, the peak value of interference would 
probably be of more significance than the average value. 

It is convenient to let 


Sj vilin 

H = 10 logio ¥ (58) 
wp 

F = —10 logi | W jx(w;) dw; (59) 


@a 


H is the ratio expressed in db of mean signal power in one channel to the 
total interference power received in one channel, and F is the weighting 


factor expressed in db. From (56), 


U=10"”™ (60) 
Equation (57) may be written in the form, 
._? 
ate! (61) 
where 
a = 8.69 er. — db 
(N — 1)r 
(62) 
b = 57.3 Vi mk 1) degrees 
ae T ] 


Without the numerical factors, a and } are expressed in nepers and radians 
respectively. 

If we regard o; and 
in which @ and 6 are the semi-axes. 
H — F we may thus represent the requirements on gain and phase variation 
by elliptical boundaries in the oy09-plane. Figure 5 shows such a diagram 
constructed for a large number of channels each active one-fourth of the 
time and with flat weighting. In terms of the symbols above, we have set 
N/(N — 1) equal to unity, 7 = 1/4,and F = 0. Gain and phase variations 
included within a particular ellipse produce average interference power less 
than the amount designated on the boundary in terms of db down on mean 
power in one channel. The requirements imposed on both gain and phase 
variation are considerably more stringent than the corresponding require- 
ments for carrier systems using frequency discrimination and employing 
comparable band widths. 

Requirements on linear transmission of the line are, of course, not the 
only considerations involved in a comparison of time division multiplex 


a2 as variables, (61) determines a family of ellipses 
By assigning values to N, 7, and 
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methods with competing methods of superimposing channels. Other 
aspects to be considered are the synchronization of transmitting and receiv- 
ing switching processes, the effects of non-linearities in the line, and the 
sensitivity of the system to external interference. It is thought, however, 
that the severe restrictions imposed on phase and attenuation characteristics 
when economy of band width is required form the weakest feature of the 
method and will in many cases provide the primary criterion for judging 
its availability in the solution of particular problems. Conversely, if the 
crosstalk requirements of the system are sufficiently mild to enable the 
transmission problem to be solved, the other problems also become relatively 
simplified. 
FURTHER REFERENCES 
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Steady State Delay as Related to Aperiodic Signals 
By R. V. L. HARTLEY 


The concepts of phase and envelope delay, as applied to any linear 
system, rather than only to a medium, are discussed. Criteria are set 
up for the time of occurrence of that part of an aperiodic signal which 
corresponds to a small segment of the spectrum. The original spectrum 
of the signal gives the time of entry and this spectrum as modified by th 
phase characteristic of the system gives the time of exit. 

If the amplitude is constant over the segment, it is shown that when 
the criterion is the time of maximum envelope of the disturbance, th: 
aperiodic delay is identical with the envelope delay. When it is the 
time of maximum absolute value, the delay depends on the signal spe 
trum, the phase shift of the system, and the envelope delay, but not on 
the phase delay. 

If the amplitude vaies rapidly with frequency, the component of an 
aperiodic disturbance which corresponds to a narrow segment of the 
spectrum persists so long that the resulting over-lapping of neighboring 
segments makes their interpretation difficult. 


ie THE earlier applications of steady state theory to transmission prob- 
lems the emphasis was placed on the variation of amplitude with fre- 
quency. The use of long loaded lines made it necessary to take account 
of phase distortion! as well. With the development of telephotography and 
television”, the phase characteristic was found to provide a useful index for 
predicting the overlapping of adjacent picture elements. For these purposes 
it has been found convenient to express the phase characteristic in terms of 
phase or envelope delay. These may be called “steady state delays”’ since 
they are defined and measured in terms of sinusoidal disturbances of ad- 
justable frequency. However, the signals for which they are intended to 
furnish an index are aperiodic in nature. It seemed worthwhile, therefore 
to examine more closely the relations existing between “aperiodic delays,” 
defined in terms of such signals, and steady state delays. 

Let us first review the development of the concepts of steady state delay 
Early in the study of the propagation of sinusoidal waves a distinction was 
made between phase and group velocity. If we fix on a particular distance 
of transmission the ratio of this distance to each of these two velocities 
may be interpreted as a delay associated with the transmission. In the 

1 For discussion and references see ‘Phase Distortion and Phase Distortion Corre: 


tion,” S. P. Mead, B. S. T. J., Vol. VI, p. 195, 1928. 
? Symposium on Television, B. S. T. J., Vol. VI, p. 551. 
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STEADY STATE DELAY AND APERIODIC SIGNALS 223 
communication art, these delays have been called phase and envelope delay, 
respectively. If the medium exhibits dispersion they vary with frequency. 
Let us fix our attention on the conditions throughout the medium at a 
particular instant during the transmission of a sinuscidal disturbance. We 
may determine the total change of phase in passing from the input to the 
output. This may be more than a single cycle. If now we divide this 
phase shift by the frequency, expressed in the same angular units, we get 
the time which will be required for the phase at the input to progress to 
the output, or the phase delay. Also it may readily be shown that the 
derivative of this phase shift with respect to frequency is equal to the 
envelope delay as defined above in terms of the group velocity. The 
simplest treatment of this is based on the consideration of two sinusoidal 
waves of equal amplitude and slightly different frequencies. 

While these delays can be easily interpreted for most media, difficulties 
arise in the case of those substances which exhibit anomalous dispersion. 
Here, in the neighborhood of certain frequencies, the phase shift varies 
rapidly with frequency, and often appears to be discontinuous. Actually 
the apparent discontinuity is a region of very rapid decrease of phase with 
frequency, which leads to a negative value of envelope delay. In the same 
region the transmission varies rapidly with frequency, and selective reflec- 
tion occurs at the entering boundary. This effect can be explained in terms 
of resonance in the elements which make up the fine structure of the 
medium. 

The next step was to dissociate the idea of delay from that of velocity 
in a medium, and associate it with a steady state transfer characteristic 
between any two points of a linear system. This would permit its appli- 
cation to all sorts of complicated networks in which uniform propagation 
cannot be readily visualized. Here two types of characteristic are to 
be distinguished. One, which is associated with what might be called 
“damped” systems, exhibits a reasonably gradual variation of both phase 
shift and attenuation with frequency. This is the analog of a medium 
having normal dispersion. The other, which is associated with ‘‘resonant”’ 
systems, exhibits the phenomena associated with anomalous dispersion. 
In the case of filters and hollow wave guides these resonances give rise to 
regions of high attenuation and reactive impedance, which are the analogs 
of the regions of high absorption and selective reflection at the boundary 
ofa medium. In applying the idea of delay to networks then, we can expect 
the results to agree with our intuitions only so long as we keep away from 
the critical frequencies of resonant systems. 

In computing or measuring the phase shift of a system, at a single fre- 
quency, the result is indeterminate so far as the addition of multiples of 
2x is concerned. This does not affect the envelope delay, which depends 
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only on the derivative, and so this type of delay can be generalized direct\y 
to include the transfer characteristics of arbitrary networks. To give a; 
exact meaning to phase delay some convention would have to be adopted 
for determining what, if any, multiple of 27 is to be added to the computed 
phase for the frequency in question. Apparently no such convention has 
been agreed upon which is of general application. For damped networks 
which transmit frequencies down to zero, it is customary to assume the 
phase shift to be zero at zero frequency, and, for higher frequencies, to ad 
multiples of 27 so that the phase shift varies continuously with frequen 
If, then, B is the computed phase shift, between —z and 7, we may repre- 
sent the continuously varying phase shift by B + 2mr, where m is the 
number of discontinuities in B which have been eliminated in passing from 
zero to the frequency in question. The phase delay may then be defined as 


_ B+ mr 


w 


D, (1 
Any similar convention for resonant systems would be less simple, and 
since, as will appear below, phase delay has little bearing on aperiodi 
signals, it seems unwise to attempt to formulate such a convention here. 
In contrast with steady state delay, let us now examine the delay of a: 
aperiodic signal. If the signal is transmitted without distortion the con- 
cept of delay of the signal as a whole is simple. If, because of distortion, 
the sent and received signals are different we may still agree upon some 
recognizable feature of each as determining its time of occurrence. If the 
distortion is considerable the delay may vary greatly with the distinguishing 
characteristic chosen. For example, if it depends on the behavior of com- 
ponents of high frequency the delay may be quite different from what it 
is if it depends on those of low. In the first case the result would be little 
affected if, before transmission, the signal were sent through a high-pass 
filter and, in the second, if it went through a low-pass filter. In each cas¢ 
we measure a delay associated with a disturbance which comprises only 
those Fourier components of the signal which occupy a particular limited 
range of frequency. We may carry this idea farther and make use of a 
very narrow band-pass filter. By varying the mid-frequency of this band 
we obtain a delay which is a function of frequency. Its value, at any 
frequency, is the delay, as defined by our convention, of a disturbance 
which corresponds to that part of the spectrum of the signal which is in 
the immediate neighborhood of the frequency in question. Our problem 
then is to find recognizable features of a disturbance of this kind such that, 
when they are used as criteria of delay, the result can be related directly 
to the phase or envelope delay as defined in terms of periodic disturbances. 
Compared with the pair of equal sinusoids used in the derivation of 
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envelope delay, this disturbance differs in that, in any finite range of fre- 
quency, there are an infinity of sinusoids, the amplitudes of which need 
not all be the same. For simplicity, we assume the actual filter to be 





replaced by an idealized one in which there is no distortion within the band 


/ and no transmission outside it. If the signal as a whole be represented by 
a Fourier integral, we may obtain the desired disturbance, for an angular 


frequency, «, by integrating from w,; — 6 tow; + 6. The disturbance may 
be represented by 


w +6 


f(t) = real part of M exp [—a + i(wt — 6)] dw, (2) 


w;—8 


where M is a constant dependent on the magnitude of the signal and a@ 
and @ are functions of frequency and position which describe the spectrum 
of the signal at various points in the system. 

The first step is to perform the indicated integration and express the 
resulting function of time in a convenient form. For this we let 


€=w-— a. 
Since we are interested only in small values of € we may replace a by 


/ 
a= a, + aye, 


, 0a — 
where a; and a; are the values of a and — at w,. Similarly, 


OW 
6 = 0+ be. 
We define an instant, 7, , by 
T.=6, (3 
and a time, 7, by 
gra2ti— 7,. (4) 


Substituting these in (2) and performing the integration, we get 


f(t) = real part of 2M exp[—a + i(wit — ((: — w163)] a = us — 
(—a, a tT) 
If we introduce the angles, 


— a 
8 = arc tan —s 


and 


tanh (—da;) 
tanér 


y = arc tan 
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and take the real part, we get 
‘ 2 ‘ , 9) 
[(cosh a1 sin 67)” + (sinh da; cos 67)" |’ 
(ay? + 7°)! 


, 
cos(w;7 — (6; — w10,) + B— yy). 5 


f(t) = 2M exp (—a) 


Let us consider first the extreme case where the spectrum of the signal js 
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Fig. 1—Elementary disturbance corresponding to a narrow segment of the spectrum 


- . . . . ° , . ry 
uniform in amplitude in the neighborhood of w;, so that a; is zero. Then 


sin 6r 


, 
fl) = 26M exp (—a}) cos (wT — (6, — w16;)). (6) 
Here the amplitude includes a constant factor which is proportional to 
the bandwidth, 26, and to the magnitude, M exp(—ay;), at the frequency, 
w;, and a function of time, a plot of which is shown in Fig. 1. This function 
consists of a sinusoidal wave of frequency, w; , the amplitude of which varies 


’ 


with time, the envelope being symmetrical about the instant, 7, = 6, 
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at which it isa maximum. 7,, the time of maximum envelope, is then a 
unique instant which is suitable for defining the time at which the dis- 
turbance occurs. It is determined solely by the slope of the phase frequency 
curve for the spectrum. 

The instant, 7,, may be interpreted, in accordance with the principle 
mal is of stationary phase, as the one at which the sinusoidal components of (2) 
are most nearly in the same phase, and so have the least destructive inter- 


27 r 








a 
PHASE IN RADIANS 
€6-——. ... = 


FREQUENCY (w) 











Fig. 2—Graphical representation of the phase of an elementary disturbance 


ference. This condition will hold when the instantaneous phase angle is 
changing least rapidly with frequency, that is, when 


hen 
0 
—- (wt — 6) = 0, 
(6) Ow 
from which 
to 
, 
Cy, i= 6;. 
T n Wh P .e ° e ° 
. he angle, 6; — w,4;, in (6), gives the phase of the wave at the instant, 
1es _ ° ° e _— ° P e ° e 
; I, , when its envelope isa maximum. The interpretation of this angle will 
6), 


be aided by the geometrical construction of Fig. 2 which is similar to that 
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employed for phase and group velocity’. The abscissae are values of , 
and the ordinates are values of phase in radians. A portion of the function, 
9, in the neighborhood of w: is shown. The distance, OB, is 6,. The slop. 
of the tangent, CA, to the curve at A is 6}. The distance, CB, is w6, 
Consequently, OC, or the intercept of this tangent on the phase axis, js 
6; — w;0;. If, as shown in the figure, the absolute value of this intercept 
is greater than 7, we may transform (6) to a form in which the angle js 
less than 7, by the substitution 


/ 
¢ = 6, — widi + 2nz, 7 
where 7 is an integer and 
le| <7. 


In Fig. 2, is 3, and ¢ is the distance DC. (6) then becomes 
f(t) = 26M exp (—a) = cos (wit — ¢), 
T 


and ¢ is the ordinary phase lag of the sinusoid, relative to an origin of time 
given by the instant of maximum envelope. 

We may choose as the instant at which the disturbance occurs, not 7, 
at which the envelope is a maximum, but 7,, at which the instantaneous 
value of the function has its maximum absolute value. Since 4 is small 
compared with w; , this will occur very nearly at the smallest absolute value 
of rt for which cos(wit — ¢) is +1. This will occur for 


Yy T us 
t = , when —.~- <ogo< 
w Sate 
and for 
xT is T 
T= when ~—x <y< —. or <9 < *#. 
@) 2 2 


From (4), (3) and (7), 
Ze 
where & is an integer such that 


TT 


7 < Y= 1 — ond + be <5. 


The significance of this can be seen from Fig. 3. Here, in addition to the 4 
curve of Fig. 2, there are plotted a series of curves whose ordinates differ 


3’ Lamb, “‘Hydrodynamics,”’ Cambridge U. Press 1916, p. 371. 
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from it by multiples of 7. In so far as any one purely sinusoidal component 
of the disturbance is concerned, values of phase determined by those curves 
which differ by an even multiple of x would be indistinguishable. Those 
differing by an odd multiple would represent a reversal of sign. Let us 


6Tr 
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Fig. 3—Graphical representation of the time of maximum absolute value 
now select that curve for which the tangent at w; intersects the phase axis 
nearest the origin, and call it @+ kx. Since, for the case drawn, 


T 


DC | < =, 
' | 2 


2n. 
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If it were greater we should have 
k = 2n +1. 


It is then obvious that the time of maximum absclute value, 7, , is give 
by the slepe of the line OE. It differs from 7, by the difference in slope oj 
the lines OF and GE. 

We have then deduced from the spectrum of the disturbance its time oj 
occurrence in terms of two definitions of the latter. The next step is t 
compare these times for the input and output and determine the corre- 
sponding delays. Let us consider first the case where the attenuation js 
independent of frequency, so that a; is zero in the output signal also. We 
may then confine our attention to the phase, @. Let us represent its valu: 
at the input by 8, and the phase shift of the system by B. Then at the 
output @ will be equal to 6 + B. If we take the time of occurrence as 
determined by the maximum envelope, these times at the input and output 
are 

T., = bi, 

Te, = +B. 
The delay is then 
1» — Te = Bi, 
which is by definition the envelope delay of the system. 

If we take the time of occurrence based on the maximum absolute value, 
we have, at the input, 


— b + loll 


a 
Ww 


where 
-" < Wo = by a wi “+ kor < 5" 


At the output, 


ts b, + Bi t+ (Ro + ks)4 


? 


(iP 1 


wh 


where 


— 3 < Ys = bi + Bi — woi(bt + Bi) + (ho + bs) < 5. 
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The delay, 
B, + k3r 


’ 


wi 


D, = Ta, — Ta = 


While there is a superficial similarity between this and the phase delay (1), 
it is of little real significance; m, in (1), is determined by the aggregate 
increase in phase shift with frequency, while &, is determined mainly by the 
rate of increase at w,;. An example of a situation in which the two delays 
are very different, is furnished by a wave guide when the frequency only 
just exceeds the cutoff. The phase delay is then almost zero while the rate 
of change of phase shift with frequency is very large. 

Thus the delay based on maximum absolute value depends on both the 
envelope delay and the phase shift of the system, but not on the phase 
delay. There remains to examine this dependence in more detail. The 
value of ks depends on the spectrum of the signal as well as the characteristic 
of the system. It is of interest to see if it can be replaced by a quantity 
derived from the system characteristic alone. The most obvious thing to 
try is a delay which is derived from the phase shift of the system in the 
same way that the time of absolute maximum is derived from that of the 
signal spectrum. This would be 


oi Bi + kor 


@) 


Ds 
where 


—"<W=B-—wBithr <". 


to 


The difference between this and the aperiodic delay based on absolute 
value is 


Ds — D, = ed (ke — ks), 


wy 


— (Wo + Vs — W;). 


@) 


Since ke — ks is either zero or an integer and | W; | is less than > 9 if 


—" Wt < 


T 
7? 


Ds — Da = 9. 
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If 
—r<Wth<—", 
Ds — Da = —™ 
Ww) 
If 
T 
2 < Vo + he < o, 
Dg ~ Da => = 
@! 


Thus the delay as derived from the system characteristic alone may ly 
identical with the aperiodic delay based on maximum absolute value or jt 


may differ from it by +—, that is by half a period. Which condition 
w1 


holds depends on the interrelation of the phase functions which characterize 
the signal spectrum at the input and the transmission of the system, and 
not on either of these functions alone. 

If the attenuation is not uniform, a; cannot be neglected and the expres- 
sion for the output signal becomes more complicated. Both the amplitude 
and phase in (5) vary with time in a manner which depends on the value 
chosen for 6. The expression becomes fairly simple, however, for the case 
where a; is very large, as in anomalous dispersion and in highly resonant 
systems. Then, even when 6 is small, we may assume that 


cosh (Sa) = exp (+ 5a), 
sinh (6a;) = + exp (+ 5a), 


, 
according as a, 2 0. 
The amplitude factor in (5) then becomes 


M exp (—ai + éa;) 
(a? + 7) 


Here the exponent is equal to the value of a at that edge of the segment of 
the spectrum where the amplitude is greatest. The amplitude is sym- 
metrical about + = 0, that is, about ¢ = 6; , at which point it has its maxi- 
mum value. Hence the instant of maximum envelope is still given by the 
slope of the phase, frequency curve, as when a; is small. However, the 
maximum is now extremely flat and its sharpness no longer depends direct}; 
on 6. Over the range of values of r for which 7? << a?, the amplitude is 
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1 
sensibly constant. When 7 = +az, it is reduced to 9 times its maximum. 
) V 


Ve « . ° 
For 72 >> ay’, it varies inversely as | 7 | . 
To investigate the oscillating factor of (5) we note that now 


y= tir +5, 


7 _ 
where the sign of 67 depends on that of ai and that of 3 does not. The 
oscillating factor then is 
— , 
cos [(w: + 5)r — (0; — wi61) — al, 
where 
, 
a) T > 
n = arctan — + -. (8) 
T 2 
The frequency, (w: + 4), is that of the edge of the segment of the spectrum 
where the amplitude is relatively very large. The phase differs from that 
, . . . . . ° 
for small values of a; by a quantity 7 which is an ambiguous function of the 
time tr. This ambiguity may be removed if we assume that the phase 
varies continuously and that, for very small values of 7, the amplitude has 
the same sign as the spectrum component corresponding to an infinitesimal 
/ 


° ay ° ° 
value of 6. As 7 increases through zero, arctan — changes discontinuously 
T 
. + # T ' es a :, oe ; ch 
from + 5 to +5 according as a; S$ 0. To avoid a similar discontinuity, 


in n we say that the sign of 5 in (8) is to be taken opposite for positive and 


negative values of r. If we make it + for 7 < 0, and + forr > 0, according 
asa; 2 0, then 7 is zero in the neighborhood of r = 0. Since the amplitude 
factor is always positive, this corresponds to a spectral component of positive 


amplitude. If we make the sign of 5 + for r < 0, and + for 7 > 0, 9 


becomes + 7, which is the equivalent of a negative amplitude. Hence a 
knowledge of the spectral component of frequency w; enables us to determine 


the sign in (8). For large values of (7), » reduces to aes. 


Here we have assumed the amplitude of the input signal to be independent 
of frequency. If this is not the case the same conditions hold at the input 
as have just been discussed for the output of a resonant system. 

The main conclusion to be drawn from the foregoing is that when the 
amplitude is changing rapidly with frequency, the component of an aperiodic 
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disturbance which corresponds to a narrow segment of the spectrum persists 
for a considerable period so that there is much overlapping of the contriby- 
tions of neighboring segments. It is therefore difficult to deduce the nature 
of the disturbance at any particular time from any narrow region of jts 
spectrum. For the same reason it is difficult to associate the delay ex- 
perienced by an aperiodic signal with the steady state characteristic of a net- 
work when the attenuation of the latter is changing rapidly with freque: Cy 

The net result of our study then is that steady state phase delay has no 
direct relation to the particular types of delay of an aperiodic signal whic! 
we have chosen to investigate. When the amplitude does not change 
rapidly with frequency, envelope delay is identical with the delay produced 
in the maximum value of the envelope of a disturbance corresponding to 
that part of the signal spectrum which is in the immediate neighborhood 
of the frequency in question. The envelope delay, together with the phase 
shift, determines the delay in the maximum absolute value of this dis- 
turbance, subject to an uncertainty of half a period. This uncertainty 
depends on the particular combination of signal spectrum and system 
characteristic. When the amplitude does change rapidly with frequency, 
the envelope delay still gives the delay in the maximum value of the envelope. 
However, this maximum is so flat that the interpretation of the results is 
very difficult. 
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Engineering Requirements for Program Transmission 
Circuits* 


By F. A. COWAN, R. G. McCURDY and I. E. LATTIMER 


Present-day program networks are reviewed from the standpoints of 
engineering, design, and operation as developed to meet the needs of 
the broadcasters. The factors requiring consideration in the further 
development of program networks in anticipation of future needs are 
also discussed. The presentation of the paper is supplemented by a 
demonstration of the quality obtainable by transmission over various 
types of telephone facilities. 


INTRODUCTION 


HE growth of radio broadcasting to the magnitude of a major national 

industry within the last twenty years has been accompanied by the 
development of a nation-wide system of wire-line networks interconnecting 
hundreds of broadcasting stations. Papers have been presented before 
this Institute from time to time!** describing the types of plant used for 
these networks and discussing important features of their design and opera- 
tion. With these twenty years of experience as a background, it should 
now be of interest to review how the various requirements of broadcasting 
have influenced the development of the networks and to consider some of 
the factors which have determined the point to which transmission and 
operating features have so far been carried. 

Simply stated, broadcasting is a means by which sounds originated at one 
place are reproduced simultaneously to large numbers of listeners distributed 
over wide areas. The simplest possible radio broadcasting system would 
consist of a microphone, a radio broadcast transmitter and some radio 
receiving sets. Such a system could serve only the listeners within the 
comparatively limited service area of the transmitter. To serve the whole 
nation many transmitters must be established about the country. Further- 
more the most desirable sources of program are not usually in the neighbor- 
hood of the transmitter to which a particular listener can tune, since talent 
tends to be concentrated in certain parts of the country, and special events 
of interest may occur anywhere. To give a true country-wide service so 
that every listener can hear the programs he enjoys wherever they may 

* Presented at A.I.E.E. Winter Convention, Philadelphia, Pa., January 27-31, 1941. 


Published in Electrical Engineering, Transactions section, April 1941. 
! For all numbered references, see list at end of paper. 
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originate, a Supplementary transmission system must be provided, inter- 
connecting the many studios and broadcasting stations. ‘The wire networks 
that perform this function comprise the subject matter of this paper. 

The present extent of the wire-line facilities which are associated with the 
major portion of these networks is indicated in Fig. 1. The width of the 
lines on this chart has been made proportional to the numbers of circuits in 
3 the various sections. The total length of these circuits is in excess of 110,000 
miles, and it is not unusual for a program originating at some point on a 
network to traverse more than 7,000 miles of circuit before being broadcast 
by the most remote station. 

The requirements which the program networks must meet are in the final 
analysis determined largely by the needs of the broadcasters. The objective 
of a program network service is to meet these needs in as complete and 
| prompt a manner as possible consistent with reasonable cost. With this 
f* objective in mind, it is necessary in planning the plant to consider not only 
-f the day-to-day needs, but the possible future needs as well. The importance 
| of this may be appreciated when it is considered that plant provided today 

for program transmission service will need to be adaptable to the service 
-f requirements ten or twenty years hence. As a result of such planning, 
7 cables and equipment installed five, ten, and fifteen years ago meet present- 
day requirements, and, with some rearrangements, will take care of those 
likely to develop tomorrow. 

The detailed planning of program transmission circuits requires considera- 
tion of: 





1. The numbers of circuits likely to be required, section by section, over 
each route; 









2. The provisions for reliability, flexibility, operation, and supervision 
essential to a high-grade network service; 
3. The transmission requirements, or electrical characteristics, necessary 






to achieve a natural reproduction of the program. 






These three general classes of requirements will be considered in order. 






NUMBER OF CIRCUITS REQUIRED 






The circuits which have been established on a full-time basis for continu- 
ing use form the backbone of the program networks. Even for these circuits, 






however, permanence is relative since frequent extensions and rearrange- 
ments are made to meet changing requirements of the broadcasters. Aside 






from these fulltime circuits there are intermittent requirements occasioned 
by special events and other short-period needs of the broadcasters, some of 
which involve networks almost as extensive as the full-time networks. 






In addition reliability of service requires provision for rerouting the net- 
works in the event of trouble. Figure 2 shows the year-by-year growth in 
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the operated mileage of program circuits for the period 1926 to 1949 
Of the more than 110,000 miles of circuits shown for 1940, about 45.00 
miles have been provided for the short-period services and as stand-by 
facilities for protection. In addition to these, there are still other circuits. 
normally assigned to other services, which are arranged to be readily adapt- 
able to program service to supplement the reserve facilities maintained on a 
full-time basis. 

The time interval necessarily accompanying any extensive constructio 
project makes it necessary to engineer plant considerably in advance of 
actual service requirements to meet, not only the expected growth, but also 
the changes in network routing. Figure 3 shows for two typical sections 
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Fig. 2—Growth in mileage of major network circuits 


along major routes the variations in requirerents for full-time network 
circuits resulting from growth and rearrangerrents required by the broad- 
casters. While, in planning to meet these rapid variations in circuit 
requirements, advantage can be tal en of sorre latitude which exists in the 
choice of routes for occasional services and protection facilities, the task of 
balancing the provision of circuits against requirer ents is an entertaining 
and at times difficult one for the circuit engineer. 


OPERATING REQUIREMENTS 


Considering for a moment the variety of programs originating at many 
different points that can be heard on any horre radio set in the course of an 
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evening without once changing the tuning, it will be apparent that minute- 
to-minute rearrangements of an established interconnecting network must 


be possible. 


For example, studios have to be changed from receiving to 


originating, sections of the network have to be made to transmit first in 
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Fig. 3—Variations in full-time program circuits 


one direction and then in the other, and branches have to be connected and 


disconnected. 


These changes in the network have to be made in the few 


seconds elapsing between the close of one program and the start of the next 
on the receipt of selected cue words or sounds. Even during the course of a 
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single program, switches or reversals may have to be made to change the 
originating point temporarily. To provide for these rapid changes, specig 


operation and special switching and reversing equipment are required a 
many points along the network. Much of this equipment is under remot, 
control from selected points. 

The greater portion of the switching of program circuits is done at about 
25 points throughout the country on the major networks. On the average 
more than 25,000 switching operations per month are performed at these 
25 points. During the busy hours of any typical evening there may be 
something over 500 men on duty at all of the offices about the networks. 

At points where switching requirements are simple, the switching equip- 
ment consists merely of a few keys. At the larger points where the switching 
requirements are complex, the switching equipment consists of elaborate 
relay and control arrangements. These are so designed that it is possible 
to set up in advance the circuit combinations required for the ensuing 
program period without disturbing the programs in progress. The actual 
switching operation takes place at the instant the monitoring attendants 
signal the the receipt of the last of selected cues, and not before then. This 
type of arrange:rent affords a maximum of protection against error, as it is 
possible to check the presetting for the next switch or make a last minute 
change if necessary any tire before the switch has been made. 





Figure 4 shows a picture of such a switching arrangement in use at 
Omaha, Nebraska for one broadcasting company. At this point 13 circuits 
used in various trunk and branch sections of two networks are connected 
to the switching equipment. These are grouped in various combinations to 





take care of as many as five simultaneous programs. A maximum of five 
cues might, therefore, be involved in a switch at this point. 

The operation and maintenance of the networks are carried out by a 
special orginization under centralized authority and trained in the applica- 
tion of uniform methods and procedures found by experience to be produc- 
tive of best results. Transmission is monitored continuously at strategic 
points about the networks. In order to facilitate the activities of this group 
many thousands of miles of intercommunicating telephone and telegraph 
circuits are provided full time for their use. 

A picture of a monitoring position in the program transmission office at 
Washington, D. C. is shown in Fig. 5. It will be noted that the monitoring 
attendant is using an individual headset. This is of a special high fidelity 
type and is used to avoid the confusion that would result from attempting 
to monitor a number of different programs simultaneously with loud- 
speakers. Loud-speakers are available, however, for supplementary checks 
of quality whenever required. 

Accurate transmission measuring equipment is necessary at the various 
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erating points about the networks to insure satisfactory transmission 


specia maintenance results. 
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Fig. 4—Switching panel at Omaha, Nebraska 


TRANSMISSION REQUIREMENTS 
The general transmission requirement for a broadcasting system is that 
the program material be transmitted with a high degree of naturalness. 
Although the exact determination of the transmission characteristics which 
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would accomplish this involves many considerations it will be assumed j; 
this discussion that an ideal transmittion system is one in which the sound 
waves impressed on the listeners’ ears in the home are an exact replica of th 
sound waves striking the microphone in the distant studio. Limitations 
inherent in the human ear, in the program material to be transmitted, and 
in the usual listening conditions, however, make such ideal transmissio 
unnecessary. In expressing the requirements for satisfactory transmission, 
frequency range, attenuation distortion, delay distortion, nonlinearity and 
noise are used as indices of quality. 





Fig. 5—Monitoring position at Washington, D. C. 


Before taking up the transmission requirements of a program circuit, it is 
important to consider further the fundamental factors that are involved in 
fixing the characteristics considered desirable for the entire system. Ac- 
cording to Harvey Fletcher,‘ the zone of audibility of the average normal! 
human ear for pure or single frequency sounds is the area within the curve 
of Fig. 6. The abscissas represent frequency and the ordinates show the 
range of intensity recognizable as sound, between the lower limit or threshold 
of audibility and the upper limit where the sensation of pain is felt. It is 
seen that the extreme frequency range shown on the chart is from about 2() 
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to 20,000 cycles per second. This range is for young people. It is consid- 
erably less for middle-aged and elderly people, and varies with individuals. 
In addition to the limitation of the ear there is the fact that there is little 
energy present in most program material in the extremes of this range, 
particularly in the upper frequencies. The energy versus frequency spectra 
of music and other forms of program have been published elsewhere.° 
Figure 7 shows the frequency range which must be transmitted for a number 
of instruments, speech, and certain noises, so that competent observers 
cannot detect any impairment.’ For whole orchestras, experiment has 
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Fig. 6—Limits of audible sound 


shown that the elimination of frequencies below 40 and above 15,000 cycles 
per second is undetectable.‘ It the upper limit of the transmitted frequen- 
cies is lowered from 15,000 cycles, the impairment is at first barely detectable 
but increases at an accelerating rate. When the limit is materially lower 
than 8,000 cycles, the loss is readily apparent to many people. 

Another important consideration is volume range—that is the difference 
between the maximum and minimum levels of the program. The ordinates 
of Fig. 6 show that for part of the frequency range, the ear can respond toa 
range of intensities of more than 120 decibels, with perhaps 100 decibels as a 
mean. However, the following considerations show that the volume range 
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which the transmission system needs to accommodate is considerably nar | 
rower than the intensity range to which the ear can respond. on 
In the first place, the range of program volumes to which the ear ca: phe 
respond is much less than the range of single-frequency intensities shown by 
the curve. Program waves are in general very irregular in shape, andeven §& hos 
at constant volume contain large and small peaks differing in amplitude by 

many decibels. The range between the volure at which the highst peaks 

reach the maximum instantaneous intensity which the ear can tolerate and 
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Fig. 7—Audible frequency ranges of music, speech, and noise 
Observers voted for entire band by ratio of 60 to 40 over band shown by extremities 
of lines, and by ratio of 80 to 20 over band indicated by circles. Broken lines show range 
of noise accompanying music. 


the volume at which the smallest peaks are just above the threshold of audi- 
bility is therefore less by a number of decibels than the intensity range of 
the ear as measured by single frequencies. 

In the second place, the volume range of the usual program material 
has definite limits. Measurements have shown that a large symphony 
orchestra produces a maximum volume range of about 70 decibels.‘ The 
volume range of most other types of program is considerably less than this, 
for example, being only about 25 to 30 decibles for dance music and as little 
as about 15 decibels for much of the dialogue of actors in radio drama. 
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In the third place, the usual listening conditions impose a definite limit 
on the useful volume range. The loudest passages in the music of a sym- 
phony orchestra correspond toa sound level of about +95 decibels at a point, 
say one-third the way back in an auditorium, but most people in their 
homes prefer a level which is lower than this by 5 to 10 decibels. Figure 8 
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Fig. 8—Residence room noise with radio sets silent 
Average 43 db. Standard deviation 5.5 db. 


shows the results of an extensive survey’ of acoustic room noise in homes. 
It will be noted that the average noise level is +43 decibels on the sound- 
level scale, and that few homes are quieter than +30 decibels even in the 
suburbs. The signal-to-noise range inherent in the listening conditions, 
and allowing nothing for the noise contributed by the transmission system 


or the room noise where the program is being produced, is therefore 
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seen to be somewhere between 45 and 65 decibels. There is, therefore. yo 
advantage to the listener in providing a permissible volume range materia||y 
wider than this in the transmission system. . 

The above discussion applies primarily to the transmission of symphoni 
and similar high-grade program material. Much program material is less 
exacting in its requirements, but on the other hand some sound effects 
such as the tearing of paper require the reproduction of higher frequencies 
for complete naturalness. 

With these broad considerations in mind, the requirements of high-quality 
program circuits may now be taken up. As has been noted, the progran 
network is but one part of the over-all broadcasting system which, in a 


dition, includes microphones and studio equipment, radio transmitters, 
and the home receivers with their loud-speakers. It may be taken as a 
goal for the program networks that their transmission be nearly enough 
distortionless so that the over-all performance in regard to naturalness oi 
reproduction will not be limited by them. 

To meet such a requirement for short program circuits having only one 
or two sections is not difficult technically and does not in general require 
costly types of plant. However, the vast country-wide program networks 
are made up of many sections of circuits in tandem, which as mentioned be- 
fore may total in some cases as much as 7,000 miles. This makes it necessary 
to design and operate the individual circuits to very close limits so that the 
cumulative discrepancies in the whole network will not exceed tolerable 
values; and to consider carefully the types of plant employed lest by virtue 
of sheer numbers of units involved, the total cost be out of line with the 
over-all grade of service being given the listner. These two conflicting 
factors are important ones in the consideration of transmission requirements 
for networks. The determination of the practical working characteristics 
of program networks involves a consideration not only of the physical and 
cost factors discussed above but also of such other factors as cost of studios, 
broadcast transmitters and receivers, and the limitations of the frequency 
allocations of broadcast stations. 

From the standpoint of frequency band the consideration of all factors has 
resulted in the major present-day program networks being set up to transmit 
a frequency band with an upper limit of about 5,000 cycles. All program 
facilities installed in the last ten years or so, however, have been designed 
to be adaptable to the future transmission of frequencies up to 8,000 cycles. 
Operation on an 8,000-cycle basis, however, requires the release of additional 
frequency space now occupied by other services in much of the plant and a 
general readjustment of the program-circuit characteristics. In 1933, 
experimental wire circuits were set up between Philadelphia and Washington 
to transmit frequency bands up to 15,000 cycles. These were employed in a 
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demonstration of stereophonic transmission and reproduction of music.‘ 
Studio-transmitter loops have been provided to transmit wider frequency 
bands than the 5,000 cycles currently provided on the nation-wide net- 
works. At the present time, many of the studio-transmitter loops are being 
set up to transmit bands up to 15,000 cycles. A demonstration will be 
given at the close of the paper of the transmission of programs over cable 
circuits about 1,200 miles in length with frequency bands extending to 
15.000, 8,000, and 5,000 cycles. The 5,000-cycle circuit is of the type in 
present commercial use. The 8,000-cycle circuit is of a type to which 
much of the present program plant can readily be modified. The 15,000- 
cycle circuit consists of a standard carrier system to which has been added 
program terminal equipment now under development. 

In the consideration of transmission requirements for program circuits 
other than nominal frequency band, the variation in performance with length 
and type of circuit is important, since the factors tending to impair trans- 
mission are in the nature of small amounts of distortion or noise which 
accumulate over the length of the circuit. If these effects varied in some 
definite manner with length, transmission requirements could be fixed on 
that basis. However, good engineering practice frequently requires choos- 
ing for the various sections of a long circuit, different types of facilities 
whose contributions to the total effects are not in proportion to their length. 
Even the determination of the maximum permissible distortion and noise 
on a circuit is influenced by outside factors such as are involved in the broad- 
casters weighing operating flexibility and cost against the frequency of 
occurrence of unfavorable network routings and the number of stations 
affected. For example, in order to secure operating flexibility with a mini- 
mum of total network mileage, most of the networks employ the so-called 
“round robin” principle for a part of the network. In this arrangement the 
circuit follows a route from station to station forming a continuous loop 
which returns to its starting point. This naturally results in increased 
circuit mileage between the program source and the more distant listeners 
with an attendant increase in undesired transmission effects. For these 
reasons no exact or specific transmission requirements can be stated for 
even the over-all performance of program transmission service. 


VOLUME RANGE 


The permissible volume range for a program circuit is determined by 
the maximum volume which can be transmitted as limited by nonlinear 
distortion or crosstalk, and the minimum volume which can be transmitted 
without impairment from the noise present on the circuit. 

In connection with their design the various types of program circuit are 
subjected to listening tests in which the transmission of program over a 
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long loop of the circuit is compared with transmission over a local distor. 
tionless circuit. Each type of circuit thus is rated as to the maximum 
volume it can transmit without noticeable distortion. The highest volume 
which can be permitted without excessive crosstalk into other program or 


message telephone circuits is also investigated, and whichever limit is the 
lower determines the maximum allowable working volume for service. ‘The 
range between the maximum permissible volume and the noise level on very 
long lengths of the present program circuits is about 45 or 50 decibels, 
except under some conditions on certain open-wire sections. On the indi- 
vidual links making up the long circuit, the range is 10 or 20 decibels greater 
than this. 


ATTENUATION AND DELAY DISTORTION 


Another important consideration is the amount of attenuation and delay 
(or phase) distortion to be permitted within the transmitted frequency band. 
It is the practice to equip program circuits with adjustable attentuation 
equalizers. By means of these once the desired frequency band has been 
chosen the deviation in attenuation at any frequency within that band, 
compared with that at 1,000 cycles, can be adjusted within close limits. On 
very long circuits, however, experience has shown that even with automatic 
regulating features and careful operation residual variations which may 
amount to several decibels may develop as a result of changing temperature 
and other conditions. These variations are kept within tolerable limits 
by readjustment of the equalizers from time to time. 

Associated with the attenuation distortion is another effect detrimental 
to program quality, namely, differences in time of transmission for different 
frequency components of the signal. In practice, circuits tend to have a 
lower velocity of transmission near the edges of the frequency band than 
in the middle portions. This results in frequency components near the 
edges being delayed as compared to the middle portions of the band. This 
difference in time of transmission is called delay distortion of the circuit. 
Careful listening tests have shown that it becomes noticeable if, at the 
highest transmitted frequency, the delay is more than eight milliseconds 
greater than at 1,000 cycles, and if, at 100 cycles, it is more than about 
15 milliseconds greater than at 1,000 cycles. It is controlled by careful 
attention to the design of loading systems, amplifiers, repeating coils, and 
all other elements of the circuit. Since such small amounts of over-all 
delay distortion are detectable and since networks frequently have 100 or 
more amplifiers in tandem between an originating point and the broad- 
casting stations on the more distant portions of the networks, it is necessary 
that the delay distortion of all individual components of a network be held 


to exceedingly close limits. Accumulations of residual delay distortion 
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PROGRAM TRANSMISSION CIRCUITS 
which cannot be entirely eliminated in design are reduced by the use of 


delay equalizers along the circuits when they are set up. 


CONCLUSION 


From this discussion it is seen that the program networks are comprised 
of many parts, each of which must meet exacting requirements in order 
that over-all results will be satisfactory. It is seen that equally important 
with transmission are the requirements for plant flexibility, adequate 
reserves, uniform practices, and centralized supervision of the networks. 

The features discussed have been those found desirable for present-day 
network service. As indicated earlier, consideration of the needs of the 
future as well as those of the present is an essential feature of the design 
and engineering of the plant for program-network service. As a result 
of having done this it will be possible to provide with present plant, and 
with new plant currently being installed, adequate network facilities as 
the broadcasting art develops toward higher standards of performance. 
With the past experience as a guide, it appears that there should be no 
fundamental difficulty in meeting all reasonable requirements, always 
remembering that in the long run, requirements and costs bear definite 
relations to each other. 
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Abstracts of Technical Articles by Bell System Authors 


Notes on the Time Relation between Solar Emission and Terrestrial Dis. 
turbances.'. CLirFORD N. ANDERSON. Although the correlation betwee 
general solar activity and terrestrial disturbances is quite evident, the 
association of individual storms with specific sunspot groups has never 
been very satisfactory. Disturbances sometimes have occurred when no 
sunspots were visible and at other times large sunspots have been unaccom- 
panied by any abnormal disturbances. A possible explanation of such 
anomalies may lie in longer transit times for the disturbing solar emission 
than is usually assumed. Some indication is given in this paper that these 
transit times may range from periods as short as only one or two days to 
as much as three months. The corresponding velocities for the above 
transit times are of the order of 2000 and 20 kilometers per second. 

Curves show the approximate relation between the angle of emission, 
velocity, day of emission, and the days intervening between the passage of a 
spot through the central meridian of the sun and the corpuscular encounter 
with the earth. 


The Effect of the Earth's Curvature on Ground-Wave Propagation.’ 
CHARLES R. Burrows and Marion C. Gray. Curves are presented for 
the rapid calculation of the ground wave for radio propagation over a 
spherical earth of arbitrary ground constants, antenna heights, and polari- 
zation. 

Based on the pioneering work of G. N. Watson, a rigorous theory of the 
propagation of electromagnetic waves round a spherical earth has been 
developed in the past twenty years. Watson developed his method in 
detail only in the limiting case of an earth of infinite conductivity, but his 
work has since been extended by various authors to cover other values of 
the earth’s conductivity. Theoretically, therefore, solutions are available 
for any values of the earth’s constants (dielectric constant and conductivity 
and for either vertically polarized or horizontally polarized waves. In 
practice, unfortunately, the computations required are lengthy and in- 
volved, and for the most part the recent theoretical papers have confined 
their calculations to a few specific values of the earth’s constants. The 
present paper attempts to summarize the results so far obtained in a manner 


1 Proc. 1.R.E., November 1940. 
? Proc. I.R.E., January 1941. 
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that will make them more easily available to the practical engineer, and to 
fill the gaps in these results by developing a series of curves from which the 
field for any values of the earth’s constants may be read, with all the ac- 
curacy that could be expected in engineering practice. 


Electrical Breakdown of Anodically Oxidized Coatings on Aluminum: A 
Means of Checking Thickness of Anodized Finishes.’ K. G. Compton and 
A. MenpizzA. The existing methods for determining the thickness of 
anodically produced oxide coatings on aluminum are relatively few and are 
almost entirely of a destructive nature. It is a fairly well established fact 
that, within the thickness limits normally encountered in practice, the 
voltage breakdown is a linear function of thickness of oxide film. The 
authors have endeavored to utilize this fact in developing a test method 
for determining the thickness of coatings produced under known and con- 
trolled conditions with practically no injury to the finish. Data are given 
which show the relationship between breakdown resistance, anodizing time, 
thickness of coating, current density and sealing of anodically oxidized 
polished commercially pure aluminum. Statistical data for the values ob- 
tained are also given, indicating the good reproducibility of the breakdown 
values. By calibrating a particular anodic process, satisfactory results 
may be obtained in a relatively short time and often without destroying or 
marring the article. Since the oxide coating is not entirely homogeneous 
it is necessary to obtain a fairly large number of readings for every test 
condition. The authors have found that approximately twenty-five read- 
ings are usually sufficient and can be made in a relatively short time. Al- 
though only one of the many anodizing possibilities has been investigated, 
the applicability of this method of evaluating the thickness of oxide coatings 
may be extended to all commercial treatments. 


Ultrasonic Absorption and Velocity Measurements in Numerous Liquids.’ 
GERALD W. WILLARD. By means of ultrasonic light-diffraction phenomena 
the velocity and absorption of sound in some forty transparent liquids were 
measured in the frequency range of 6to 30 Mc. Among the list of materials 
studied are mixtures of liquids in varying proportions, several solutions of 
solids in liquids, and a non-liquid jell. A novel-construction glass-to-metal- 
to-quartz cell made possible the study of highly solvent liquids. Velocity 
values were obtained from measurements of the diffraction spectra spacing. 
Absorption values were obtained by measurement of the sound radiator 
voltages required to produce certain color transmission effects at measured 
distances from the sound radiator. The use of a mercury arc light-source 


*A.S.T.M. Proc., Vol. 40, 1940. 
‘Jour. Acous. Soc. Amer., January 1941, 
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enhanced the necessary color effects. The relation between sound heam 
width (in the optical direction) and light transmission was studied. |; 
general, the values of velocity obtained were found to be independent oj 
frequency, and the absorption to be proportional to frequency squared and 
unrelated to calculated viscous and thermal losses. A simple calculatio 


is proposed for estimating absorption errors caused by sound beam diffra: 


tion and spreading. These apply as well to absorpiion measured in ot] 


methods than here used. 
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